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SECTION  I 

NTRODUCTION  AND  SUMMARY 


SCOPE  OF  STUDY 

This  is  the  final  technical  report  for  the  FLIR  Stabilization  Study  Program 
(Contract  No.  DAAK70-8I-C-0255).  This  document  describes  the  result'  of  a  5- 
month  technical  effort  to  study,  analyze  and  conduct  preliminary  design  of  a 
precision  pipeline  correlation  for  FLIR  stabilization  applications.  The  period 
covered  by  this  report  is  15  October  1981  to  15  March  1982.  Mr.  Peter  Durenec  is 
the  contract  monitor. 

The  overall  objective  of  this  program  is  to  define  the  appropriate 
architecture  and  design  of  a  pipeline  correlator  for  FLIR  stabilization,  analyze 
expected  performance  and  investigate  other  potential  applications.  The  pipeline 
correlator  is  an  implementation  of  a  precision  correlation  algorithm,  proprietary 
to  CAI,  which  measures  image  shifts  to  accuracies  beyond  the  resolution  of  the 
sensor/detector. 


SUMMARY  OF  RESULTS 

Correlation  stabilization  is  essentially  different  from  inertial  stabilization, 
in  that  with  the  former,  the  imaged  scene  itself  is  tracked,  as  opposed  to  sensor 
inertial  rotations,  as  with  the  latter.  This  is  shown  to  provide  significantly 
improved  performance  at  lower  cost. 

The  salient  results  of  this  program  are  as  follows: 

e  A  modularized  approach  has  been  defined  which  allows  various 

configurations  dependent  upon  application,  and  is  capable  of  expansion 
from  60  to  120  to  180  elements  per  line. 

e  A  basic  stabilization  module  has  been  defined  and  designed  (preliminary) 
which: 


a)  Provides  scan  and  cross  scan  direction  interframe  offsets  on  a 
pixel-by-pixel  and  line-by-line  basis,  respectively 


b)  Can  be  used  with  fixed  or  servo-controlled  sensors 

c)  Can  replace  up  to  three  gyro  transducers  and/or  a  scan 
position  sensor 

d)  Provides  significantly  better  performance  than  gyro  transducers 

•  higher  bandwidth 

•  better  accuracy 

•  no  limit  to  positional  excursion 

•  wider  range  of  angular  rates 

•  no  drift 

•  lower  cost 

e)  Provides  edge  video  output 

f)  Provides  a  measure  of  focus 

#  The  basic  stabilization  module  contains  two  cards  for  pipeline 
correlation,  one  frame  store  controller  card,  and  one  and  one-half  cards 
per  60  elements  of  detector  array  length. 

•  The  addition  of  one  moving  target  indication  (MTI)/range  card  and  one 
control  receiver  card  to  the  stabilization  module  and  a  remote  control 
panel  provides  added  capabilities.  This  configuration  will  provide  video 
format  outputs  for  display  or  use  by  another  processor  of  tlie  following: 

a)  MTI 

b)  Either  relative  motion  mapping  or  (given  detector  size,  focal 
length  and  range  inputs  from  the  control  panel)  absolute  velocity 
mapping 

Given  two  boresighted  sensors,  but  not  requiring  the  frame  storage 
cards,  this  same  configuration  will  outpuf  passive: 

c)  "Within  range"  indication 

d)  Either  relative  range  mapping  or  (given  detector  size,  focal  length 
and  baseline  separation  inputs  from  the  control  panel)  absolute 
range  mapping 


A  special  function  frame  store  module  has  been  defined  and  designed 
(preliminary)  which: 

a)  Provides  video  putput  in  any  combination  of: 
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•  Passively  stabilized  for  translation,  dilation,  rotation  and 
skew  (update  occurs  only  for  those  portions  of  the  stabilized 
scene  currently  within  the  FOV  of  the  sensor;  other  portions 
are  frozen  with  last  valid  video) 

•  Integrated  in  time  to  improve  signal-to-noise  ratio 

•  Scaled  1:1,  2:1  or  4:1 

•  Dynamically  sampled  to  improve  resolution  up  to  3.5x 

b)  Provides  a  window  superimposed  on  sensor  video  to  define  special 
video  location  and  coverage 

c)  Is  used  with  a  plug-in  addition  to  the  remote  control  panel  and  the 
control  receiver  card  option  in  the  stabilization  module 

d)  Consists  of  one  controller  card  plus  one  card  per  60  elements  of 
detector  array  length. 

OPTIONAL  CONFIGURATIONS 

Figures  l.l  through  1.5  illustrate  the  several  configurations  possible  with 
these  modules  and  optional  boards. 

GUIDE  TO  THE  REPORT 

A  detailed  description  cf  the  stabilization  module  and  a  comparison  of  its 
performance  to  that  of  gyro  transducers  are  given  in  section  II.  This  section  also 
presents  some  correlation  concepts  pertinent  to  this  report  for  the  reader 
unfamiliar  with  correlation.  Also  included  in  section  II  is  a  description  of  other 
outputs  provided  with  the  basic  module  and  an  explanation  of  the  need  for  a 
dedicated  frame  store. 

One  optional  configuration  which  provides  motion  mapping  of  passive  ranging 
is  given  in  section  III,  together  with  a  short  explanation  of  passive  ranging. 

Section  IV  is  a  brief  discussion  on  the  relative  merits  of  active  end  passive 
stabilization. 


BASIC  STABILIZATION 
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STABILIZATION  PLUS  MTI 


FIGURE  1.2 


PASSIVE  RANGING 


T  TO  SCAN  MIRROR  SYNCHRONIZER 
nLINE  BY  LINE  RANGE 


FIGURE  1.4 


STABILIZATION  WITH  SPECIAL  FUNCTION  AND 


FIGURE  1.5 


A  description  of  the  features  offered  by  the  optional  special  function  module 
is  given  in  section  V. 

In  section  VI,  the  design  consideration  and  the  conclusions  of  this  preliminary 
design  effort  are  presented.  A  more  detailed  discussion  of  the  recommended 
brassboard  demonstration  system,  together  with  the  specific  card  design,  is  in 
section  VII,  while  production  design  considerations  are  discussed  in  section  VIII. 

Conclusions  and  recommendations  are  given  in  section  IX. 


SECTION  II 

STABILIZATION  MODULE  DESCRIPTION 


BENEFITS 

The  essential  feature  of  correlation  stabilization  is  that  the  stabilization 
results  are  derived  from  the  scene,  rather  than  an  inertial  frame  of  reference. 
This  is  the  source  of  its  many  advantages  over  gyro  transducers.  These 
advantages  are  summarized  in  table  II- 1 . 

For  example,  since  gyro  transducers  measure  only  sensor  rotations,  they  do 
not  account  for  sensor  translation.  Therfore,  on  a  moving  platform,  they  are  not 
effective  in  tracking  a  particular  portion  of  the  environment.  Correlation 
stabilization  conversely  tracks  the  scene  itself,  and  therefore  would  command  the 
sensor  to  rotate  and  maintain  that  scene  in  the  field  of  view  as  the  platform 
translates  relative  to  the  scene.  Note  that  if  the  scene  is  o  moving  target, 
correlation  stabilization  will  track  the  target. 

Because  it  operates  on  the  image,  correlation  stabilization  is  inherently 
tuned  to  the  sensor  and  image  parameters.  A  change  of  focal  length  has  no  effect 
on  performance  relative  to  resolution,  except  for  the  necessity  to  scale  any 
feedback  to  servo  mechanisms.  Accuracy  is  rated  in  terms  of  image  resolution 
and,  therefore,  the  necessary  precision  to  stabilize  a  scene  is  maintained 
regardless  of  sensor/detector  parameters. 

Correlation  provides  more  information  for  stpbilization  than  do  gyro 
transducers.  It  provides  data  on  translation  motions  in  two  axes  and  rotation. 
Hence,  it  actually  takes  the  place  of  three  gyro  transducers.  Furthermore,  it 
tracks  skew  and  dilation  or  scaling  of  the  image  in  two  axes  and  can  passively 
stabilize  the  image  for  these  effects,  which  is  impossible  with  inertial  devices. 


The  B-2  digital  scan  converter  for  FLIR  imagery  includes  a  scan  position 
sensor  (SPjj)  module  to  measure  the  position  of  the  scan  mirror  in  common  module 
FLIR's.  This  module  is  used  to  assure  accurate  location  of  lines  for  reformatting 


PIPELINE  CORRELATION  FOR  FUR  STABILIZATION 

ADVANTAGES 


PROVIDES  OTHER  FUNCTIONS 
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to  video  standards.  Correlation  stabilization  can  replace  this  module.  However, 
it  should  be  understood  that  correlation  precisely  locates  lines  in  relative  terms, 
on  a  frame-to-f’"jme  basis.  It  will  not  correct  any  geometric  distortion  which  may 
have  been  present  in  the  first  frame.  Therefore,  when  such  distortion  caused  by 
the  scan  mirror  is  not  considered  important,  the  correlation  module  will  obviate 
the  need  for  the  SPS.  It  will  "correct"  subsequent  frames  to  the  original 
distortion,  if  any. 

In  every  performance  figure,  correlation  stabilization  is  significantly  better 
than  gyro  transducers.  It  is  more  accurate  and  has  higher  bandwidth  and  greater 
range.  It  can  also  be  located  remotely  from  the  sensor  in  a  protected 
environment,  therby  allowing  smaller,  more  rugged  sensor  designs. 

Furthermore,  the  results  can  be  used  actively  by  feeding  back  to  servo 
mechanisms  on  the  sensor,  or  passively  by  correction  of  pixel  locations  in  frame 
store  or  on  display.  Finally,  as  will  be  detailed  later,  correlation  stabilization 
provides  many  other  important  functions. 


CORRELATION  CONCEPTS 

These  paragraphs  are  intended  to  give  the  reader  unfamiliar  with  correlation 
some  basic  understanding  of  the  process,  so  that  the  pipeline  correlator  may  be 
better  evaluated.  A  more  detailed  treatment  is  given  in  appendixes  A  and  B. 

Correlation  is  a  measure  of  match  between  two  signals  or  functions.  For  the 
application  considered  here,  the  signals  are  video  waveforms  representing 
images.  Consider  two  images  on  transparencies.  Holding  the  transparencies,  one 
on  top  of  the  other  up  to  some  light,  the  best  match  occurs  when,  overall,  the 
most  light  is  transmitted  by  the  pair  of  images.  For  example,  if  the  images  are  of 
the  same  scene,  but  there  is  a  relative  shift  or  translation  between  the  two,  then 
by  sliding  one  transparency  over  the  other  and  measuring  the  total  light 
transmitted,  the  relative  shift  can  be  measured  by  noting  where  the  peak  in 
transmitted  light  occurs. 
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I 
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I 
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Correlation  is  merely  the  mathematical  implementation  of  this  process.  It  is 
the  multiplication  of  the  transmittance  (or  intensity)  of  a  point  on  one  image  by 
the  transmittance  of  the  point  on  the  second  image  which  overlays  the  first  point. 

The  multiplications  are  take  point-by-point  and  summed  over  the  whole 
image.  As  one  image  is  shifted  over  the  other,  new  multiplications  and  sums  are 
taken  and  plotted  as  a  funtion  of  shift.  Locating  the  peak  correlation  value 
defines  the  shift  between  the  images  as  shown  in  figure  2.1. 

At  the  peak  of  the  correlation  plot  or  curve,  all  spatial  frequencies  match  up 
and  contribute  to  making  the  peak  what  it  is.  As  the  images  are  shifted  away 
from  this  optimum  shift,  the  highest  spatial  frequencies  no  longer  match  up,  but 
lower  frequencies  do  still  match  fairly  well  and  contribute  to  an  above-average 
correlation  value.  Thus,  the  correlation  curve  as  a  function  of  shift  slopes  more 
or  less  gradually  down  from  the  peak,  depending  on  the  frequency  content  of  the 
image. 

If  the  two  images  are  strictly  copies  of  one  another,  the  correlation  curve 
obviously  peaks  at  zero  shift.  This  is  called  autocorrelation.  Similarly,  if  the  two 
images  are  identical  but  shifted  versions  of  one  another,  the  correlation  curve 
peaks  at  the  relative  shift  position  but  is  identical  in  shape  as  shown  in  figure  2.2. 

Obviously,  correlation  cannot  truly  be  taken  point  by  point  since  there  are  an 
infinite  number  of  points.  The  signals  must  be  sampled  to  obtain  a  finite  number 
of  discrete  points.  However,  as  soon  as  the  signals  are  sampled,  the  full 
correlation  curve  can  no  long  be  produced.  Only  a* sampling  of  the  curve  results. 
For  autocorrelation,  the  zero  shift  sample  is  still  located  at  the  true  peak  of  the 
curve,  but  this  true  peak  on  the  cross  correlation  curve  is  sampled  only  if  the 
relative  shift  is  an  integer  number  of  sample  spacings  (pixel  dimensions).  Since 
the  object  is  to  find  the  shift,  it  cannot  be  known  a  priori  that  the  shift  is  indeed 
an  integer  number  of  pixels. 

Simply  locating  the  peak  sample  only  defines  the  shift  to  an  accuracy  of  *1/2 
pixel.  However,  since  it  is  known  that  the  autocorrelation  curve  samples  and  the 
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FIGURE  2.1 
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FIGURE  2.2 


cross  correlation  curve  samples  are  samples  of  the  same  curve,  it  is  possible  to 
curve  fit  the  samples  to  obtain  a  far  more  precise  estimate  of  peak  position  on  the 
cross  correlation  curve,  that  is,  to  a  smaller  fraction  of  a  pixel  dimension.  Thus, 
it  is  not  high  spatial  frequencies  that  are  matched  in  precision  correlation,  but  the 
lower  spatial  frequencies  of  the  scene.  A  summary  of  the  correlation  process  is 
given  in  table  11-2. 

In  fact,  precision  pipeline  correlation  implements  only  this  curve  fitting 
operation.  Thus,  there  is  an  implicit  assumption  that  the  peak  is  actually  within 
+2  pixels  of  shift  from  where  it  is  expected.  This  is  of  little  concern  for  the 
application  of  stabilization,  since  the  shift  can  be  tracked  from  line  to  line  and  on 
this  basis,  it  is  nearly  impossible  to  shift  more  than  2  pixels  per  line  (i.e.,  within 
35  ps). 

One  final  point  is  that  the  correlation  curve  sharpness  defines  a  measure  of 
focus.  If  the  scene  is  in  focus,  there  are  high  spatial  frequency  components  to 
match  up  and  produce  a  high  zero  shift  autocorrelation  value.  Conversely,  if  the 
zero  shift  value  is  not  much  larger  than  the  neighboring  shift  values,  there  is  no 
high-frequency  information,  or  it  has  been  destroyed  by  loss  of  focus.  As  a 
natural  by-product  of  pipeline  correlation,  the  autocorrelation  width  is  calculated 
and  can  be  made  available  to  a  focus  drive  mechanism  in  the  sensor  as  shown  in 
figure  2.3. 


PIPELINE  CORRELATION  STABILIZATION 

Precision  pipeline  correlation  is  a  hardware  implementation  of  correlation 
matching  wherein  each  operation  is  sequentially  performed  along  a  bus.  Each 
pixel  datum  enters  the  bus  in  sequence,  is  operated  upon,  and  is  passed  to  the  next 
operation.  Each  operation  may  use  constant  operands  and/or  results  from  other 
operations.  Upon  exiting  the  bus,  the  datum  has  effectively  been  transformed  into 
a  measure  of  offset  or  misregistration  relative  to  the  datum  for  that  pixel  from 
the  previous  frame  as  shown  in  figure  2.4. 
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PIPELINE  ARCHITECTURE 
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FIGURE  2.4 


Correlation  is  defined  over  the  bounds  of  some  window  less  than  or  equal  to 
the  field  of  view.  This  window  may  be  arbitrarily  set,  depending  on  the 
application.  In  pipeline  correlation,  the  transformation  of  data  is  broken  at  a 
specific  point  in  the  pipeline  to  be  summed  over  the  desired  window  before 
proceeding  to  a  final  processing  stage  to  produce  an  offset  datum.  If  the 
summation  is  implemented  as  a  running  sum,  always  subtracting  the  partial  results 
calculated  n  cycles  ago,  results  can  be  produced  at  video  rates.  This  is  illustrated 
in  figure  2.5.  In  fact,  each  result  corresponds  to  a  window  centered  about  each 
pixel  in  the  original  video  as  shown  in  figure  2.6. 

The  geometry  of  common  module  FL1R  imaging  is  particularly  suited  to 
pipeline  correlation  stabilization  in  several  respects.  Essentially,  it  consists  of  a 
linear  array  of  up  to  180  detectors  which  are  scanned  back  and  forth  across  the 
scene.  The  scan  mirror  is  tilted  slightly  at  either  end  of  a  one-directional  scan  to 
provide  interlaced  fields.  Each  one-directional  scan  defines  a  field  and  two 
interlaced  fields  define  a  frame. 

Note  that  no  particular  control  is  used  in  the  pipeline  architecture  to 
separate  lines,  fields  and  frames.  The  operation  simply  accepts  a  continuous 
stream  of  video  data.  This,  of  course,  can  produce  confusion  at  the  boundaries. 
The  back  and  forth  scanning  of  common  module  FLIR's  obviates  this  problem  at 
least  in  one  direction.  More  importantly,  this  scanning  method  allows  tracking  of 
offsets  from  line  to  line  across  field  boundaries.  This  is  very  desirable  for  pipeline 
correlation  as  referred  to  previously. 

Because  the  sensing  elements  are  in  a  linear  array  format,  effects  of  motion 

I 

are  common  to  an  entire  line.  This  means  that  cross  scan  offsets  may  be 
calculated  on  a  line-by-line  basis  and,  in  fact,  there  are  sufficient  detectors  in  the 
line  to  permit  the  window  for  cross  scan  correlation  to  be  one  line  wide.  The 
result  is  that,  for  the  cross  scan  direction,  there  is  no  need  for  a  running  sum  to  be 
taken,  and  a  simplified  design  may  be  implemented.  However,  for  the  scan 
direction,  possible  skews  and  rotations  require  a  block  window  calculation  as 
shown  in  figure  2.5.  Thus,  the  architecture  required  for  FLIR  stabilization  is  as 
shown  in  figure  2.7.  Note  that  edge  images  and  measures  of  focus  are  natural  by¬ 
products  of  pipeline  correlation. 


FLIR  STABLIZATION  MODULE 
PIPELINE  CORRELATOR  BLOCK 


FIGURE  2.7 


GYRO  STABILIZATION 


Gyroscopes  are  devices  which  measure  the  change  in  angular  position  or  (not 
both)  the  angular  rate  to  which  they  are  subjected.  Modern  miniature  gyros  are 
small  canisters  which  contain  a  small  spinning  mass  (rotor).  In  the  angular  positon 
measuring  mode,  if  the  device  is  tilted,  the  rotor  tends  to  remain  at  its  original 
orientation  and  the  difference  in  orientation  between  the  rotor  and  the  canister  is 
measured  by  pickoffs.  In  the  angular  rate  mode,  the  rotor  orientation  is  locked 
relative  to  the  canister  and  only  the  torque  of  attempted  deflection  resulting  from 
an  angular  rate  is  measured. 

Gyro  performance  is  related  to  its  cost.  For  example,  a  simple,  miniature, 
rate-integrating  gyro,  which  costs  about  $2000,  typically  has  G-sensitive  drifts  to 
25°  to  50°/hr/G  and  -G-sensitive  drifts  of  30°  to  60°/hr/G.  Angular  rates  of 
response  on  these  gyros  are  in  the  range  of  1 00°  to  200°/s. 

More  sophisticated  gyros  cost  about  $10,000  and  have  non-G-sensitive  drifts 
of  l5°/hr/G  and  G-sensitive  drift  of  about  l4°/hr/G.  A  portion  of  this  high  cost  is 
attributed  to  the  gas  spin-axis  bearings  which  give  a  mechanical  time  constant  of 
a  little  over  I  ms,  that  is,  the  time  in  which  a  gyro  responds  to  an  excitement. 

The  gas  bearings  also  eliminate  much  of  the  gyro  noise.  Such  gyro's  can  handle 
rate  inputs  up  to  I  l5°/s.  The  rate  limiting  factor  in  most  gyros  is  the  torque 
motor  current  limit  due  to  wire  size# 

Non-G-sensitive  drift  is  due  primarily  to  spring  in  the  internal  leads  that 
carry  the  electrical  signals  into  and  out  of  the  gyro.  This  drift  can  ue  in  a  positive 
or  negative  direction  and  varies  from  gyro  to  gyro*  It  can  be  measured  and 
trimmed  to  zero  electrically  or  magnetically. 

G-sensitive  drift  is  due  to  a  mass  unbalance  incurred  in  the  mechanical 
construction  of  the  gyro.  For  example,  if  on  the  floated  rotor  of  a  gyro  there  is  a 
heavy  spot  on  the  side  marked  as  an  "A"  in  figure  2.8,  it  will  cause  the  rotor  to 
rotate  about  its  output  axis  due  to  the  effect  of  gravity  or  acceleration.  G- 
sensitive  drift  can  be  negated  by  mounting  a  gyro  so  its  output  axis  is  parallel  to 
the  forces  of  acceleration.  Another  way  to  cancel  these  forces  would  be  to  use  an 
accelerometer  to  supply  a  nulling  signal.  As  each  gyro  is  different,  the  scaled 
signal  must  be  differentially  tuned  for  each  device. 


Other  causes  of  gyro  drift  are  a  result  of  bearing  frictions,  magnetic 
influences  and  internal  thermal  gradients.  These  influences  are  normally 
considerably  less  than  the  G  and  non-G-sensitive  drifts,  but  still  must  be 
considered  in  the  design  of  an  accurate  system.  Nearly  all  can  be  minimized,  if 
not  eliminated  by  careful  circuit  design  and  gyro  construction,  as  is  done  for  the 
inertial  navigation  systems  on  aircraft  and  missiles  requiring  extreme  accuracy. 

The  floated  rate-integrating  gyro  is  ruggedized  and  usually  can  withstand  an 
impulse  acceleration  of  200-300  G's  before  being  damaged.  When  a  gyro  is 
damaged,  it  generally  involves  the  output  axis  bearings.  These  bearings  are  either 
of  the  needle-jewel  type  or  small  precision  ball  or  roller  bearings.  It  should  be 
realized  that  a  gyro  sitting  on  end  and  falling  over  to  its  side  on  a  hard  table  can 
easily  experience  shocks  of  300-350  G's,  so  that  the  use  of  the  word  "rugged"  is 
relative.  This  is  illustrated  in  figure  2.9. 

The  schematic  of  a  typical  gyro  is  shown  in  figure  2.10.  A  typical  block 
diagram  for  an  angular  rate  sensing  gyro  is  shown  in  figure  2.1 1.  The  gyro  and  its 
circuitry  should  be  visualized  simply  as  one  block  of  a  servoed  system  normally 
entitled  "single  axis  angular  rate  transducer".  Figure  2.12  is  a  block  diagram  of  an 
angular  position  transducer  using  a  gyro.  Figures  2.1 1  and  2.12  are  only 
representative  and,  depending  upon  the  magitudes  of  the  quantities  to  be 
measured,  may  be  implemented  quite  differently.  Table  11-3  lists  some  of  the 
advantages  and  disadvantages  of  gyros. 

COMPARISON  OF  TECHNIQUES 

I 

Inertial  stabilization  is  a  considerably  different  technique  than  scene 
correlation.  Nevertheless,  they  can  be  quite  easily  comDared.  For  the  most  part, 
measures  of  gyro  transducer  performance  are  fixed  by  the  device.  The  only  real 
exception  is  positional  accuracy  and  maximum  excursion  in  the  position  mode. 
While  there  are  physical  limits  within  the  device  on  excursion  (2°  to  3°),  it  is 
possible  to  induce  a  centering  torque  and  effectively  scale  down  the  rotational 
deviation  pf  the  rotor  relative  to  the  canister.  Thus,  greater  excursions  can  be 
allowed  at  the  expense  of  accuracy.  This  also  means  that  the  indicated  position 
will  drift  back  to  zero  while  the  device  actually  remains  at  its  new  orientation. 
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GYRO  ANGULAR  RATE  TRANSDUCER 


FIGURE  2.11 


GYRO  POSITION  TRANSDUCER 


FIGURE  2.12 


TABLE  11-3 
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HAS  BOTH  ELECTRICAL  AND  MECHANICAL  CONSTRUCTION  PROBLEMS 
VERY  EXPENSIVE  (FOR  A  GOOD  ONE) 

REQUIRES  COMPLEX  CIRCUITRY  TO  OPERATE 
MUST  BE  LOCATED  AT  THE  UNIT  BEING  STABILIZED 


Performance  figures  for  scene  correlation,  on  the  other  hand,  are  intimately 
related  to  the  imaging  sensor /detector.  A  direct  comparison  is  acheived  by 
specifying  imaging  parameters  and  must  be  repeated  on  a  case-by-case  basis.  This 
is  something  of  a  moot  point  in  the  stabilization  of  imaging  sensors,  since  the 
required  performance  is  itself  a  function  of  these  parameters. 

The  computer  simulation  of  pipeline  correlation,  discussed  in  appendix  C,  has 
shown  that  the  rms  accuracy  of  this  technique  is  conservatively  on  the  order  1/20 
of  a  detector  dimension.  This  assumes  at  least  4-bit  digitization,  a  signal-to-noise 
ratio  of  no  less  than  10:1  in  the  video  signal,  and  a  window  of  at  least  64 
elements.  Even  a  signal-to-noise  ratio  of  2.5:1  allows  accuracy  to  1/8  pixel.  The 
rms  error  of  1/20  pixel  is  used  as  the  measure  of  performance  for  angles  and 
angular  rate  accuracies.  The  mean  of  errors  were  shown  to  be  virtually  zero  over 
100  calculations,  indicating  that  there  is  no  drift  with  the  correlation  technique. 
Finally,  the  fact  that  pipeline  correlation  can  lock  on  to  offsets  up  to  +2  pixels 
allows  calculation  of  maximum  rate,  since  the  offset  must  deviate  by  more  than  2 
pixels  from  line  to  line  to  prevent  tracking  of  offsets.  Using  these  results,  table 
11-4  presents  a  comparision  of  gyro  and  correlation  stabilization. 

Several  other  considerations  have  been  mentioned,  such  as  the  ruggedness, 
remote  location,  accounting  for  translation  and  skew,  moving  scenes  and  cost.  In 
each  case,  the  arguments  favor  pipeline  correlation.  The  only  penalties  occur  in 
size,  weight  and  power.  These  must  be  offset  against  the  extra  capabilities 
provided  in  addition  to  stabilization,  such  as  motion  mapping,  passive  ranging  and 
dynamic  sampling. 

I 

In  general,  a  pipeline  correlator,  as  compared  to  gyro  transducers,  offers 
better  performance,  lower  cost  and  greatly  expanded  capabilities  at  the  price  of 
size,  weight  and  power. 
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GYRO  TRANSDUCER  MAY  MEASURE  P05IT10N  OR  KAlt 
EFFECTIVE  3.2  kHz  IN  SCAN  DIRECTION  DUE  TO  DELAY 


OTHER  OUTPUTS 


The  pipeline  correlation  process  performs  several  calculations  which  are 
useful  for  applications  other  than  stabilization. 

For  example,  a  measure  of  edge  is  performed  in  both  horizontal  and  vertical 
axes.  These  measures  are  tapped  for  possible  use  as  edge  maps  of  the  video.  If 
two-axis  edge  mapping  is  desired,  the  edge  outputs  may  be  thresholded  and  the 
results  used  as  inputs  to  a  logical  OR  gate.  The  results  are  output  in  the  same 
video  format  as  the  sensor  output. 

A  second  function  calculated  is  a  measure  of  autocorrelation  width  which  is 
equivalently  a  measure  of  focus,  as  noted  in  the  discussion  on  correlation.  There 
are  other  causes  of  broadened  autocorrelation  peaks,  such  as  image  smear. 
However,  the  object  of  this  module  is  to  eliminate  such  motions;  thus,  the  first 
task  is  to  stabilize  the  image  and  then,  if  the  image  is  still  blurred  ,  the  measure 
of  focus  output  will  permit  sharpening  the  focus  by  means  of  a  servo-controlled 
focus  drive  on  the  sensor.  The  measure  of  focus  outputs  are  given  in  both  axes.  A 
mismatch  in  signals  indicates  some  effect  other  than  defocus  to  be  the  cause  of 
blur.  Thus,  only  the  signal  indicating  sharpest  focus  should  be  used  in  correcting 
focus. 

Several  other  features  require  optional  boards  and/or  modules  and  will  be 
discussed  in  later  sections. 


NEED  FOR  FRAME  STORAGE 

Because  correlation  stabilization  is  based  upon  comparing  sequential  frames, 
it  is  necessary  that  the  video  signal  be  stored  for  comparison  with  the  next 
frame.  It  is  also  necessary  that  the  current  video  signal  and  the  stored  previous 
video  signal  be  presented  to  the  correlator  in  the  same  format.  This  fact  leads  to 
the  requirement  for  two  frame  storage  planes.  If  both  inputs  to  the  pipeline 
correlator  are  to  be  in  the  original  FLIR  sensor  format,  then  one  frame  store  to 
preserve  a  frame  for  correlation  and  a  second  frame  store  to  reformat  for  display 
are  needed.  If,  on  the  other  hand,  both  inputs  are  to  be  in  the  final  output  display 
format,  then  two  frame  stores  are  still  required:  one  to  reformat  and  one  to  store 
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the  previous  frame  in  reformatted  format.  Since  two  frame  stores  are  required, 
regardless  of  configuration,  the  better  approach  is  to  use  the  FLIR  sensor 
format.  Not  only  is  the  back  and  forth  scanning  with  vertical  line  format 
preferrable  to  simplify  the  implementation  of  pipeline  correlation,  but  the  delay 
of  reformatting  is  avoided,  allowing  optimal  bandwidth  of  the  correction  signal. 
Frame  store  configurations  are  presented  in  figure  2.13. 

The  question  arises  whether  a  DSC  reformatter  could  be  used  as  a  simple 
frame  store  for  correlation.  Although  attractive,  such  a  solution  would  be, 
unfortunately,  impractical.  The  preferred  input  for  the  stabilization  module  is 
digital;  this,  however,  is  a  minor  point  since  there  is  little  doubt  digital  output 
from  the  reformatter  could  easily  be  supplied.  The  most  serious  defect  in 
attempting  to  use  a  DSC  reformatter  is  the  rigidity  of  its  operation.  While  it  can 
be  reprogrammed  for  specific  video  formats,  it  cannot  easily  handle  changing 
address  offsets  on  a  pixel-by-pixel  basis  as  required  to  track  offsets  between  the 
current  and  previous  frames.  Finally,  the  amount  of  interfacing  between  pipeline 
correlator  and  frame  store,  as  well  as  modification  to  the  latter  if  it  were  a  DSC 
reformatter  would,  in  the  end,  result  in  a  frame  store  dedicated  to  pipeline 
correlation. 

It  is,  then,  easier  and  more  efficient  to  design  a  frame  store  specifically  for 
this  application. 

FRAME  STORE  CONSIDERATIONS 

The  pipeline  correlator  tracks  offsets  from  frame  to  frame.  If  there  are 
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sections  of  imagery  which  do  not  overlap  the  previous  frame,  the  correlator  can 
provide  no  usable  offset  data  and  an  effective  dead  time  occurs.  Naturally,  if  the 
sensor  is  not  actively  stabilized,  such  occasions  are  bound  to  occur.  Some  dead 
time  is  tolerable.  So  long  as  the  shift  which  occurs  duing  the  dead  time  is  within  a 
few  pixels,  the  correlator  will  be  able  to  pick  up  where  it  left  off  as  shown  in 
figure  2.14. 
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FIGURE  2.14 


TABLE  11-5 

2  3 

Q  (degrees  /second  ) 


WFOV 

NFOV 

360  lines 

15,900 

1,770 

480  lines 

21,200 

2,360 

Assuming  sinusoidal  excitation,  the  peak  product  of  acceleration  and 
velocity  occurs  at  45°.  Therefore,  the  maximum  angular  rate  would  be  J2  times 
the  square  root  of  the  product.  Table  11-6  gives  the  maximum  allowable  angular 
rates  of  sinusoidal  excitation,  above  which  offset  lock  may  be  broken. 

TABLE  11-6 

Maximum  Sinusoidal  Rate 

WFOV 

NFOV 

360  lines 

1 38°/s 

59.5°/s 

480  lines 

206°/s 

68.6°/s 

A  similar  condition  occurs  at  initiation  when  there  is  no  previous  frame 
against  which  to  correlate.  This  is  equivalent  to  an  effective  infinite  angular 
rate.  To  assure  initial  offset  lock,  the  system  merely  operates  as  though  the 
second  field  was  a  reversed  copy  of  the  first  field.  With  the  back  and  forth  scan 
of  the  scan  mirror,  there  is  only  a  I  /2-pixel  cross  scan  offset  at  turnaround. 
Offsets  may  easily  be  tracked  from  a  near  null  condition  by  this  method. 
Therefore,  even  if  the  conditions  become  sufficiently  severe  to  break  offset  lock, 
the  system  will  simply  reinitialize  automatically. 

All  of  the  foregoing  analysis  is  based  upon  frame-to-frame  correlation. 
Clearly,  if  this  frame  store  were  allowed  a  stabilized  mode  in  which  the  scene  is 
locked  and  the  sensor  is  not  actively  stabilized,  the  system  would  quickly  lose 
offset  track  and  be  unable  to  recover  without  locking  in  a  fresh  scene.  Scaling  an 
image  in  the  memory  would  create  a  more  severe  situation,  whether  or  not  the 


sensor  was  actively  stabilized  since  a  long  dead  time  is  assured  for  this 
condition.  For  this  reason,  no  options  of  passive  stabilization  or  scaling  are 
allowed  with  this  frame  store.  Because  passive  stabilization  is  not  feasible  for 
this  storage  plane,  neither  is  scene  integration.  The  frame  store  must  be  a  simple 
temporary  storage. 

Manipulations  of  the  imagery  suggested  above  will  be  handled  by  a  separate 
frame  store  which  is  not  used  for  correlation,  but  will  use  correlation  results. 

This  special  function  frame  store  will  be  described  in  section  V. 

Being  unstabilized,  video  data  is  written  into  the  basic  frame  storage  just  as 
it  appears  from  the  sensor,  that  is  without  offset.  However,  to  track  offsets,  the 
memory  is  read  with  offsets  dictated  by  the  correlation  results.  Because  cross 
scan  offsets  are  calculated  line  by  line  and  with  no  delay,  the  bandwidth  of  this 
offset  is  equal  to  the  line  rate  reciprocal.  The  offset  of  each  line  determines  the 
offset  with  which  to  read  subsequent  lines  for  correlation.  In  the  scan  direction, 
however,  the  windows  are  9  lines  wide.  Thus,  even  though  the  data  output  rate  is 
equal  to  the  pixel  rate,  the  results  are  delayed  4.5  lines.  The  higher  rate  of  output 
allows  correction  for  skew  and  rotation,  but  the  delay  results  in  a  response 
bandwidth  of  the  line  rate  reciprocal  divided  by  4.5. 

These  considerations  are  summarized  in  table  11-7, 


MTI  CONCEPTS 


SECTION  III 
MTI /RANGING  OPTION 


Scene  correlation  is  really  far  more  versatile  than  has  been  described  in 
section  il.  Figure  3.1  illustrates  several  possible  applications. 

Correlation  stabilization  may  easily  be  considered  to  be  the  tracking  of  the 
scene  within  the  FOV  as  a  whole.  Due  to  the  fact  that  perturbations  may  occur 
within  a  frame,  the  correlation  stabilization  described  here  actually  maintains  a 
running  track  on  subportions  of  the  field  of  view. 

The  perturbations  are  not  necessarily  due  to  sensor  motion,  but  may  in  fact 
be  due  to  scene  motion.  If  one  or  more  small  portions  of  the  scene  continuously 
exhibit  a  perturbation  considerably  different  from  the  scene  as  a  whole,  they  may 
safely  be  considered  moving  targets.  Thus,  the  correlation  stabilization  results 
may  be  used  for  moving  target  indication  (MTI)  by  subtraction  of  the  average 
frame-to-frame  offset.  The  difference  is  then  compared  to  a  threshold  value  and 
any  point  sufficiently  different  is  illuminated  as  a  movmg  target.  This  concept  is 
shown  in  figure  3.2. 

At  this  time,  it  should  be  noted^that  the  pipeline  correlation  described 
calculates  cross  scan  offset  on  a  line-by-line  basis  where  the  window  is  exactly 
one  line  in  length  and  width.  Clearly,  this  is  of  little  value  for  detecting  cross 
scan  motion.  More  desirable  is  a  square  window.  Thus,  the  partial  results  for  the 
cross  scan  motion  must  be  tapped  off  before  summation  over  the  I -line  window 
and  instead  summed  over  a  square  window.  This  provides  running  results  for  each 
pixel  in  the  field  of  view. 

The  offset  calculated  by  the  pipeline  correlator  is  proportional  to  angular 
rate  of  motion.  There  is,  in  fact,  no  need  to  subtract  the  average  and  threshold 
the  differences  if  the  application  is  more  suited  to  a  relative  motion  map  than 
MTI.  Thej/ideo  format  offset  results  may  be  taken  directly  to  provide  a  video 
image  of  motion  where  brightness  is  proportional  to  velocity.  Alternatively  the 


IMAGE  REGISTRATION  TECHNIQUES 
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SIGNATURE  MATCHING 

•  IDENTIFIES  TARGET  BY  COMPARISON  TO  "LIBRARY"  IMAGE 

f/Gl/fl£  5./ 


physical  velocity  may  be  calculated  by  multiplying  the  offset  by  the  constant  of 
proportionality  (Cv)  which  is  given  as: 


where: 

P 

fl 

R 

F 


P  R 
fl 


x  F 


=  pixel  size 
=  focal  length 
=  range  to  target 
=  number  of  frames/s 


The  absolute  velocities  may  be  used  either  for  another  processor,  or  the 
system  could  actually  be  configured  to  display  the  velocity  of  a  particular  target 
on  a  control  panel  display  (to  be  described  later).  Alternately,  this  velocity  in 
relative  terms  could  be  fed  to  the  sensor  servomechanisms  to  implement  tracking. 

The  possible  ouputs  which  are  made  available  by  this  simple  addition  to  the 
stabilization  module  are  shown  in  figure  3.3. 

The  performance  of  this  option  is  easily  calculated  from  the  angular  rate 
accuracy  quoted  above  for  stabilization  performance.  The  accuracy  may  be 
thought  of  as  noise  equivalent  angular  rate.  Neglecting  factors  of  oblique  viewing 
angles,  this  works  out  to  .51  mi/h  per  1000  ft  of  range  for  the  wide  field  of  view, 
and  .17  mi/h  per  1000  ft  of  range  for  the  narrow  field  of  view.  Similarly,  the 
maximum  velocities  tolerated  are  determined  by  the  limit  of  2  pixels  of  offset. 
These  relationships  are  graphed  in  figure  3.4. 


PASSIVE  RANGING  CONCEPTS 

Thus  far,  only  the  correlation  of  time  sequential  images  which  leads  to 
various  forms  of  motion  analysis  has  been  discussed.  The  process  could  be 
conceived  as  the  correlation  of  two  video  signals  derived  from  two  sensors 
displaced  in  time.  Consider  instead  the  correlation  of  two  video  signals  derived 
from  twoeensors  displaced  in  space;  that  is,  separated  by  a  fixed  distance.  If 
bo.  sighted,  they  will  image  the  same  scene  but  from  slightly  different  angles  as 


MTI  PROCESSING 


FIGURE  3.3 


DETECTION  LIMITS 


range  ctetaw o 
FIGURE  3.4 


shown  in  figure  3.5.  It  is  easily  shown  by  similar  triangles: 

_R_  _  F 

B  '  D 

and:  R  =  -g— 

where: 

R  =  range  to  target 

B  =  separation  between  lenses 

F  =  focal  length 

D  =  apparent  image  shift 

The  absolute  range  to  a  target  may  thus  be  calculated  using  stereo  sensors  and 
maintaining  precise  alignment  between  the  two  sensors.  Since  the  accuracies  of 
offset  measurement  are  approximately  8  prad,  such  precision  is  clearly  beyond  the 
limits  of  mechanical  tolerance.  These  misalignments  may,  nevertheless,  be 
measured  and  corrected  by  presenting  each  sensor  with  an  artificial  infinity 
target.  Such  techniques  are  known  as  autocalibration,  and  are  required  only  if 
absolute  range  calculations  are  desired.  This  concept  is  shown  in  figure  3.6. 

It  is  important  that  the  detectors  have  fixed  geometric  positions  to  the 
baseline  axis,  so  that  image  offsets  may  be  reliably  calculated.  For  this  reason, 
the  baseline  separation  must  be  in  the  cross  scan  direction.  The  scan  direction 
offsets  are  then  used  to  synchronize  scan  mirrors. 

Since  the  baseline  separation  and  focal  length  are  essentially  fixed 
quantities,  range  error  (  AR)  is  determined  by  differentiating  R  =  BF/D: 

RF 

A  R  =  .  A  D  =  _  _R_  A  D 

where: 

AD  =  error  in  measuring  D  the  image  shift 


RANGING  AND  TRACKING  GEOMETRY 


FIGURE  3.5 


FALSE  MISREGISTRATION  RESULTING  FROM 

AN  OPTICAL  WEDGE 


Since  D  has  been  defined  as  offset  (in  pixels): 
D  =  OxP 


where: 


A  D  =  Px  AO 


=  offset 


P  =  pixel  dimension 

AO  =  error  in  offset  calculation 

Thus,  the  scaling  factor  for  the  reciprocal  of  offset  is  = 

n  2 


and  the  ranging  error  is  AR 


TTxT  p  *°  = 


R*  IFOV 


Note  that  the  same  options  are  available  for  ranging  as  for  motion 
detection.  A  threshold  may  be  set  such  that  only  objects  within  a  specific  range 
are  illuminated.  Alternatively,  as  shown  in  figure  3.7,  the  offsets  may  be  taken  as 
is  for  relative  range  mapping  or  may  be  scaled  appropriately  to  obtain  absolute 
range  values.  Figure  3.8  shows  expected  ranging  accuracies  in  such  cases.  The 
accuracy  in  ranging  is  clearly  different  than  active  methods  which  measure  the 
time  between  the  transmission  of  some  signal  and  its  refelected  return.  Active 
methods  typically  have  fixed  accuracies  regardless  of  range,  while  passive 
stereometric  ranging  has  a  range  squared  dependence.  While  ranging  performance 
is  not  exceptional  because  less  than  an  optimal  configuration  is  used,  it  is 
reasonable,  especially  compared  to  recognition  range.  There  is  little  value  in 
knowing  precise  range  if  it  is  not  known  what  is  being  fired  upon.  Since  the  error 
in  ranging  is  statistical  rather  than  systematic,  considerable  improvement  can  be 
obtained  by  averaging  results.  Results  are  produced  for  each  field,  therefore  a 
10:1  accuracy  improvement  is  obtained  by  averaging  100  samples  yielding  a 
ranging  bandwidth  of  0.3  Hz  for  a  given  target.  This  improved  accuracy  is  also 
graphed  in  figure  3.8  by  reading  the  accuracy  scale  on  the  right  side.  Absolute 
range  values  may  be  fed  to  another  processor  or  displayed  on  the  control  panel  by 
windowing  a  select  target. 

Since  both  video  signals  are  supplied  from  separate  sensors  to  the  pipeline 
correlator;  there  is  no  need  for  frame  storage.  This  is  true  only  if  passive  ranging 


RANGE  PROCESSING 


FIGURE  3.7 


is  the  sole  function.  If  ranging  and  some  motion  analysis  are  desired,  the  frame 
store,  a  third  pipeline  correlation  card,  and  perhaps  the  second  MTI/range  card 
would  be  required. 

CONTROL  PANEL 

For  the  basic  stabilization  module,  very  little  control  is  required  and  no 
control  panel  is  needed.  With  the  MTI/ranging  option,  however,  a  remote  control 
panel  becomes  more  of  a  necessity.  Even  for  the  simple  MTI/"within  range" 
modes,  some  threshold  must  be  supplied.  When  absolute  values  are  desired,  the 
appropriate  terms  to  calculate  a  scale  factor  need  to  be  supplied. 

For  these  purposes,  a  simple  control  panel  consisting  of  a  keyboard,  an 
alphanumeric  display,  a  threshold  potentiometer  and  an  on/off  switch  is 
recommended.  This  panel  allows  remote  operation.  The  display  also  allows 
alerting  the  operator  to  possible  error  conditions  in  the  module  functions. 

An  additional  plug-in  control  panel  can  be  used  to  window  a  selected  target 
and  display  on  the  alphanumeric  display,  absolute  velocities  or  ranges.  This  plug¬ 
in  panel  will  also  be  used  in  conjunction  with  the  special  function  module  and  will 
be  described  more  fully  in  section  V.  For  the  purpose  under  discussion  here,  part 
of  its  function  is  to  define  a  window  on  the  video  display  and  position  that  window 
to  a  desired  location.  The  operator  must  locate  the  window  over  the  selected 
target  and  request  velocity  and/or  range  of  the  target  on  the  keyboard.  The 
appropriate  datum  of  the  absolute  value  output  from  the  motion/range  map  will  be 
latched  and  displayed  on  the  alphanumeric  display^  Figure  3.9  illustrates  the 
several  modes  of  operation. 

When  the  control  panel  is  used,  there  must  also  be  a  control  receiver  card  in 
the  stabilization  module  to  interface  the  function. 


SECTION  IV 

ACTIVE  VERSUS  PASSIVE  STABILIZATION 


There  are  two  means  by  which  imagery  may  be  stabilized: 

•  Active  stabilization  controls  the  sensor 

•  Passive  stabilization  manipulates  the  image  after  the  fact 

Active  stabilization  requires  servo  control  of  the  sensor  which  results  in 
larger,  more  complex  sensors  and  higher  probability  of  damage.  The  advantages 
are  that  it  eliminates  blur  and  allows  continuous  video  update  of  the  scene  of 
interest.  In  addition,  active  control  permits  tracking  of  particular  objects  and 
proper  search  of  a  region  larger  than  the  field  of  view. 

Passive  stabilization  also  has  advantages.  It  is  not  always  feasible  or  even 
possible  to  equip  a  sensor  with  active  stabilization.  Passive  stabilization  permits 
remote  location  from  the  sensor.  Furthermore,  in  some  cases  the  image  is  already 
being  manipulated  for  other  reasons,  and  it  may  be  a  trivial  addition  to  stabilize 
the  image  as  well  in  memory  or  on  a  display.  Passive  stabilization  can  also 
stabilize  for  image  motions  outside  the  realm  of  those  normally  handled  by  active 
means,  such  as  rotation,  skew  and  dilation.  Finally,  a  peculiar  advantage  of 
passive  stabilization  is  the  ability  to  stabilize  one  portion  of  the  scene  while 
allowing  the  sensor  to  wander  about  and  simultaneously  display  other  nearby  areas 
of  possible  interest. 

The  latter  point  suggests  that  active  and  passive  stabilization  are  not 
mutually  exclusive  and  may  in  fact  be  advantageously  used  together  in  certain 
situations.  Table  IV- 1  summarizes  the  above  points. 


COMPARISON  OF  STABILIZATION  MODES 
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MAY  BE  USED  TOGETHER 


SECTION  V 

SPECIAL  FUNCTION  MODULE  -  DESCRIPTION 


BENEFITS 

In  many  instances,  imagery  is  of  little  or  no  value  due  to  lack  of 
stabilization,  poor  scale,  poor  signal-to-noise  ratio,  insufficient  resolution  or 
fleeting  collection  time.  The  special  functions  module  provides  improved  imagery 
when  it  is  not  feasible  or  possible  to  address  these  problems  with  sensor 
modifications. 

For  example,  if  the  sensor  is  fixed  rather  than  servo  controlled,  imagery  may 
still  require  stabilization.  The  special  function  module  provides  passively- 
stabilized  imagery.  If  the  area  of  interest  is  just  a  small  area  of  the  field  of  view, 
the  area  may  be  blown  up  in  the  special  function  module.  Low  signal-to-noise 
ratio  can  be  improved  by  longer  dwell  or  integration  time.  However,  there  are  so 
many  constraints  involved  with  this  parameter,  it  is  generally  impossible  to 
change.  The  module  provides  an  artificial  dwell  time  increase  by  integrating 
frames  in  memory  with  no  effects  on  system  operation.  Insufficient  resolution  is 
also  addressed  with  the  special  function  module  as  it  includes  dynamic  sampling 
capability,  the  integration  of  detector  outputs  by  phase.  Finally,  the  problem  of 
fleeting  scenes  is  generally  handled  by  freeze  frame  operation.  This  module 
incorporates  and  improves  upon  the  operation  by  allowing  image  update  so  long  as 
the  sensor  contains  some  or  all  of  the  frozen  frame  within  its  field  of  view. 

In  general,  the  special  function  module  provides  usuable  imagery  in  even  the 
most  unfavorable  circumstances.  Its  advantages  are  summarized  in  table  V-l. 


DESCRIPTION  OF  FUNCTIONS 

The  special  function  module  is  an  optional  module  which  can  be  easily 
connected  to  the  stabilization  module  and  provides  passively-stabilized,  scaled, 
integrated .qnd/or  dynamically-sampled  imagery.  Each  of  these  functions  can 
significantly  enhance  FLIR  operation  under  degraded  conditions. 


SPECIAL  FUNCTION  MODULE  ADVANTAGES 


The  stabilization  function  makes  use  of  the  offset  data  provided  by  the 
stabilization  module  to  write  fresh  video  information  into  the  appropriate  memory 
location,  i.e.,  that  memory  location  which  contained  the  particular  scene  element 
on  the  previous  frame.  For  the  first  frame,  the  information  is  written  in  normal 
order.  Thereafter,  write  addresses  are  modified  by  the  correlation  offset 
results.  Reading  video  data  to  the  reformatter  for  display  is  always  in  normal 
order  without  offset.  Therefore,  as  shown  in  figure  5.1,  the  output  continues  to 
provide  the  orginal  scene  with  regard  to  the  position  of  scene  elements,  and  is 
merely  updated  in  time  for  fresh  intensity  information. 

Because  the  stabilization  module  requires  unstabilized  frame  storage,  offset 
results  relate  only  to  frame-to-frame  differences.  As  was  shown  in  section  II,  if 
the  angular  rate  is  constant,  the  offsets  are  constant  on  unstabilized  imagery. 
Taking  this  one  step  further,  with  constant  angular  rate  applied  to  the  sensor,  the 
correlator  results  will  be  zero.  To  track  the  position  of  the  stabilized  scene  in  the 
special  function  module  memory  relative  to  the  sensor,  this  module  must 
accumulate  correlator  offset  results,  thereby  calculating  the  total  offset  since 
locking  onto  a  scene.  Figure  5.2  is  a  block  diagram  of  this  concept. 

Whenever  the  total  offset  added  to  the  current  pixel  and  line  addresses 
constitute  a  valid  memory  address,  that  address  may  be  updated  with  current 
video.  In  other  words,  updating  of  stpred  scene  information  is  allowed  only  when 
and  where  the  current  scene  from  the  sensor  overlaps  the  scene  stored  in 
memory.  Whenever  such  updating  is  allowed,  the  current  scene  information  is 
brightened  a  few  gray  levels  to  display  the  location  of  the  stabilized  scene  on  the 
unstabilized  scene  as  shown  in  figure  5.3. 

Should  the  offset  tracking  performed  by  the  pipeline  correlator  ever  break 
lock  under  the  conditions  described  in  section  II,  the  correlator  will  reinitialize 
itself,  but  all  reference  between  stabilized  and  unstabilized  scenes  will  have  been 
lost.  In  this  case,  the  special  function  frame  store  merely  maintains  the  stored 
scene  in  a  true  freeze  frame  mode,  and  a  message  indicating  the  loss  of  reference 
will  be  displayed  to  the  operator  on  the  control  panel  alphanumeric  display. 


FIGURE  5.2 


STABILIZATION  WINDOW  DISPLAY 


WITH  WINDOW  01 
STORED  SCENE 


Image  scaling,  by  itself,  is  simply  the  writing  of  each  pixel  datum  into  a 
block  of  memory  locations.  An  alternative  method  would  be  to  re-read  pixels  and 
lines  several  times  before  proceeding  to  the  next  pixel  or  line.  Although  greater 
flexibility  might  be  achieved  by  this  approach,  implementation  is  more  difficult 
and  it  would  conflict  with  the  control  required  for  dynamic  sampling.  As  will  be 
shown,  dynamic  sampling  requires  filling  the  memory  with  the  scene  portion  of 
interest. 

To  scale,  it  is  necessary  to  define  what  portion  of  the  scene  is  of  interest. 
Among  the  many  methods,  the  easiest  is  to  display  a  window  on  the  unsealed  video 
of  appropriate  size  and  allow  the  operator  to  move  the  window  with  a  joystick  or 
similar  control.  The  scale  factor  determines  the  displayed  window  size.  When  the 
window  is  satisfactorily  located,  a  control  button  enables  the  system  to  select 
that  window  to  display  an  enlarged  format.  This  is  shown  in  figure  5.4. 

Because  memory  chips  can  only  store  or  recall  one  datum  at  a  time,  it  is  not 
practical  to  .attempt  writing  a  pixel  datum  to  several  memory  locations  in  one 
chip.  For  this  reason,  the  frame  storage  architecture  is  specifically  designed  as  a 
4x4  array  of  chips  wherein  each  chip  stores  only  every  fourth  pixel  and  every 
fourth  line.  Hence,  a  4X  scaling  is  accomplished  by  writing  a  pixel  datum  to  the 
same  address  in  all  16  memories  simultaneously,  as  shown  in  figure  5.5.  Similarly, 
2X  scaling  is  accomplished  by  writing  to  only  four  memories  simultaneously.  This 
architecture  does  limit  scaling  factors  to  I,  2  or  4x. 

A  similar  constraint  is  that  the  precision  of  stabilization  is  still  limited  to 
+  1/2  real  pixel  despite  the  scaling  of  pixels.  If  fractional  pixel  offsets  were 
allowed,  there  could  be  up  to  four  different  addresses  to  the  memory  and  would 
require  independent  address  lines,  or  the  writing  would  have  to  be  multiplexed. 
There  is  insufficient  time  for  multiplexed  writing  and  independent  address  lines 
would  be  detrimental  for  size  and  complexity  considerations.  Greater  precision  in 
stabilization  is,  nevertheless,  accomplished  with  dynamic  sampling. 


SPECIAL  FUNCTION  MODULE  MEMORY  ARCHITECTURE 


The  reading  of  scaled  imagery  is  performed  in  normal  fashion  just  as  though 
the  scaled  portion  of  the  scene  were  the  entire  scene. 

It  is  well  known  that  the  sampling  of  signal  frequencies  whose  periods  are 
shorter  than  two  sample  spacings  will  result  in  aliasing  -  the  folding  of  high- 
frequency  information  to  lower  frequency  output  patterns.  In  images,  such 
patterns  are  known  as  Moire  patterns.  The  limiting  frequency  is  known  as  the 
Nyquist  frequency.  Because  of  aliasing,  it  is  not  possible  to  reproduce  an  original 
signal  whose  frequency  is  greater  than  the  Nyquist  value  from  its  samples.  This 
does  not  mean  that  frequencies  at  or  below  the  Nyquist  limit  can  in  fact  be 
reproduced,  as  will  be  shown. 

The  sample  spacing  for  solid-state  detection  is  fixed  by  the  physical  size  of 
the  elements.  It  may  then  appear  that  the  resolution  is  essentially  fixed  by  the 
Nyquist  limit.  This  is  not  the  case.  In  fact,  frequencies  up  to  3.5X  the  Nyquist 
limit  may  be  resolved  without  aliasing  merely  by  oversampling.  That  is,  shifting 
the  scene  by  fractions  of  a  pixel  and  constructing  an  image  from  the  various 
phases  allows  improved  resolution.  The  construction  of  one  image  from  the 
independent  phase  versions  may  be  done  in  one  of  two  ways: 

•  By  considering  each  sample  to  be  a  full  pixel  in  width  but 
averaging  overlapping  fractions  of  a  pixel  according  to  phase 

•  By  considering  each  pixel  to  be  representative  of  a  fraction  of  a 
pixel  and  simply  interleaving  the  results  according  to  phase 

Figure  5.6  illustrates  each  method.  Note  that  even  at  Nyquist  frequency,  under 
static  sampling,  there  may  be  no  information  at  a  particular  phase  ($  =  1/2),  or 
there  may  be  a  slight  phase  shift  and  loss  of  contrast  ($  =  1/4  and  3/4).  The  two 
results  show  that  there  is  improved  contrast  with  the  interleaving  method,  but  no 
decision  should  be  made  with  this  single  case. 

Figure  5.7  compares  the  two  methods  for  several  cases,  sinusoid  and  square 
waves  at  2.67X  the  Nyquist  frequency,  a  single  I /4-pixel  wide  bar  and  a  three-bar 
target  at  3X  the  Nyquist  frequency.  Interleaving  generally  shows  better  contrast, 
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but  often  with  a  phase  shift  or  contrast  reversal.  The  last  example  of  a  three-bar 
target  clearly  shows  the  problem  -  spurious  resolution.  Signal  power  is  shifted 
such  that  two  instead  of  three  bars  appear.  Integration  by  phase  does  show 
correct  phase  relationships,  but  the  contrast  is  generally  so  poor  as  to  be  useless. 
In  fact,  the  three-bar  target  illustrates  that  such  loss  of  contrast  can  even  mask 
the  fact  that  there  are  three  bars,  leaving  only  a  general  hump.  It  is  clear  that 
dynamic  sampling  is  not  a  perfect  solution,  but  it  is  a  definite  improvement. 
Consider  the  first  two  periodic  signals.  Static  samples  clearly  show  aliasing  and 
any  sense  of  a  high-frequency  signal  has  been  lost,  while  the  dynamically-sampled 
version  clearly  resolves  the  pattern  even  if  it  is  spurious  resolution  in  the 
interleaving  case.  The  method  of  interleaving  appears  to  be  a  preferred  approach. 

Fortunately,  interleaving  is  the  most  simply-implemented  method  of 
dynamic  sampling.  Each  pixel  and  line  address  plus  offset  is  merely  scaled  by 
some  scaling  factor  when  writing  data  into  the  frame  store.  Reading  is  performed 
in  normal  fashion.  There  are  several  points  to  be  made  regarding  this  method. 
First,  there  ore  no  restrictions  on  scale  factors  since  only  one  writing  of  data  is 
attempted  for  each  pixel.  Second,  stabilization-mode  operation  is  required  and  it 
is  inherently  tuned  to  the  scale  factor  in  terms  of  precision.  This  is  in  contrast  to 
simple  scaling.  Third,  there  may  be  gaps  in  the  information  stored  if  particular 
phases  of  offset  do  not  occur. 

A  fourth  point  must  also  be  considered.  There  must  be  some  image  motion 
for  dynamic  sampling  to  work,  that  is,  if  there  is  active  stabilization  it  must  be 
imperfect.  The  motion  may  easily  be  induced  and  controlled  if  there  is  active 
stabilization.  Without  active  stabilization,  however,  the  system  will  have  to  rely 
on  the  random  motions  of  the  sensor  induced  by  the  environment.  The  gaps  in 
information  stored  may  be  covered  somewhat  by  storing  the  first  frame  in  a 
standard  scaled  fashion  and  thereafter  updating  those  phases  which  occur.  In  this 
manner,  there  will  be  some  reasonable  data  in  all  memory  locations,  although  it 
does  impose  the  scaling  limitations  previously  mentioned. 


Integration  of  imagery  is  the  storing  of  a  weighted  average  of  the  current 
.  frame  and  the  stored  frame.  The  relative  weights  must  add  to  unity  while  the 
ratio  of  weights  is  a  measure  of  persistance.  To  simplify  hardware 
implementation,  only  eight  possible  weighting  ratios  have  been  allowed:  from  1:0 
(no  integration)  to  1:127  (maximum  integration).  Since  the  video  is  in  8-bit  digital 
format,  a  ratio  of  more  than  1:127  would  be  equivalent  to  a  freeze  frame  since 
curent  video  would  contribute  less  than  a  least-significant  bit. 

Integration  is  defined  by  the  equation: 

S'  =  2"n  C  +  (I  -  2“n)  S 

where: 

S'  =  new  intensity 
C  =  current  pixel  intensity 
S  =  old  stored  intensity 
N  =  persistance  measure  (from  0  to  7) 

Figure  5.8  shows  a  block  diagram  of  this  function.  Note  that  the  information 
from  any  one  frame  is  never  dropped  from  the  calculation,  it  merely  decays. 
Nevertheless,  the  improvement  in  signal-to-noise  ratio  is  equivalent  to  a  straight 
average  of  2n  frames.  The  improvement  in  signal  to  noise  ratio  is  then  >/2rt  as 
shown  in  table  V-2. 

Without  stabilized  (either  actively  or  passively)  imagery,  integration  merely 
results  in  smear.  With  stabilization,  integration  can  provide  considerably 
improved  performance  in  low  signal-to-noise  ratio  conditions.  Integration  may  be 
used  with  simple  scaling  or  dynamic  sampling  as  desired. 

CONTROL  PANEL 

The  control  panel  for  the  special  function  module  is  a  plug-in  addition  to  the 
basic  control  panel  previously  described.  It  includes  switches  for  scaling,  (1:1,  2:1, 
or  4:1)  stabilization  (unstabilized,  stabilized,  dynamically  sampled),  and  an  eight- 
position  switch  for  integration  persistance  an  a  window  control  area.  For  window 
position  control,  a  four-button  array  is  used  as  a  joystick.  It  is  both  less  expensive 
and  more  rugged  than  a  joystick  and  is  easily  operated  with  a  single 


FRAME  INTEGRATION 


finger  or  thumb.  An  additional  pushbutton  control  releases  positional  lock.  Upon 
release,  the  positioning  control  becomes  active  and  the  window  may  be  moved  to 
the  desired  location.  When  the  positioning  control  buttons  are  released,  the 
window  position  is  fixed. 

These  controls  allow  complete  control  of  the  special  function  module.  The 
windowing  position  control  may  also  be  used  for  displaying  the  velocity  or  range  of 
a  particular  target  with  the  appropriate  system  configuration.  In  this  mode,  the 
window  positioning  control  is  active.  If  a  particular  window  is  locked  in  the 
special  function  module,  it  will  be  unaffected.  The  window  merely  defines  which 
output  of  the  MTI/range  card  will  be  latched  and  displayed  on  the  alphanumeric 
display  of  the  basic  control  panel.  This  mode  is  made  operational  by  a  control 
button  on  the  keyboard  of  the  basic  control  panel.  If  constants  need  to  be  input  to 
derive  absolute  velocity  in  range,  the  control  panel  will  so  inform  the  operator. 


SECTION  VI 

DESIGN  CONSIDERATIONS 

MODULARITY  AND  VERSATILITY 

The  basic  premise  of  the  common  module  program  is  the  ability  of  the  user 
to  select  the  particular  components  which  meet  a  set  of  needs,  and  to  easily 
combine  them  into  a  system. 

These  modules  attempt  to  preserve  such  versatility.  Each  block  operation  is 
a  mother  board  configuration.  Frame  store  extensions  from  60  to  120  to  180  pixel 
arrays  are  handled  merely  by  adding  the  necessary  number  of  memory  cards  to  the 
mother  board.  The  whole  frame  store  could  be  used  as  an  independent  stand-alone 
frame  store  if  desired  for  some  purpose.  Similarly,  each  optional  configuration 
described  involves  merely  connecting  the  appropriate  mother  boards  or  inserting 
them  in  their  slots.  For  example,  if  only  a  relative  motion  of  range  map 
capability  were  desired,  the  basic  remote  control  panel  module  may  or  may  not  be 
needed,  depending  on  the  application.  The  basic  stabilization  will  work  perfectly 
well  without  this  panel,  but  when  it  is  needed,  no  system  modification  is 
required.  Each  mother  board  may  be  housed  in  one  or  separate  card  racks. 

The  special  function  module,  providing  passive  stabilization,  scaling, 
integration  and  dynamic  sampling,  is  not  required  in  general  for  stabilization 
applications.  Such  functions  should  only  be  needed  in  fairly  hostile 
environments.  Since  it  is  expected  that  these  functions  will  be  used  less  than 
those  provided  with  the  stabilization  module,  it  is  provided  as  an  optional  plug-in 
addition.  While  it  seems  an  unlikely  situation,  this  module  could  be  configured  as 
a  stand-alone  module  given  some  other  stabilization  signal  with  which  to  stabilize 
imagery  and  perform  dynamic  sampling  or  integration. 

Another  form  of  versatility  has  been  achieved  by  including  multiple 
functions  or  outputs  wherever  possible.  In  general,  there  has  been  little  or  no 
additional  complexity  involved. 


By  a  modular  approach  to  card  layout  on  a  mother  board  module,  the  system 
is  designed  for  easy  growth.  Each  card  performs  a  specific  function  and  is  in  a 
sense  a  submodule.  The  addition,  deletion  or  replacement  of  functions  will  thus 
have  minimal  impact  on  the  entire  system. 

The  control  panel  has  also  been  designed  for  versatility  and  growth.  Being 
based  on  a  single  chip  microprocessor  control  and  having  an  alphanumeric  readout, 
a  wide  variety  of  control  functions,  error  messages  and  operational  cues  can  be 
implemented  with  virtually  no  effect  on  system  interface.  Commands  and 
controls  between  the  panel  and  the  operational  hardware  are  transmitted  and 
received  via  a  three-wire  bus  in  a  serial  format.  Interface  with  I553B  bus 
architecture  has  also  been  considered  and  can  be  readily  implemented. 

As  can  be  seen,  the  concept  of  modularity  has  been  a  prime  concern  in  the 
design  of  the  modules  described. 

SPEED,  SIZE  A ND  COST 

In  the  design  of  a  pipeline  processor,  speed  of  operation  becomes  a  serious 
concern.  Each  operation  must  be  completed  within  I  pixel  time.  After  being 
multiplexed  to  a  serial  format,  a  pixel  time  may  be  on  the  order  of  193  ns  when 
the  format  is  a  field  of  180  pixels  by  480  lines. 

The  operations  to  be  performed  are  addressing,  storing  and  retrieving  of 
information  from  memory,  arithmetic  operations  and  various  forms  of  encoding 
and  decoding  of  control  conditions.  For  current  memories,  access  times  of  150  ns 

I 

ore  considered  fast  but  not  unusual  in  today's  devices.  For  arithmetic  operations, 
additions  are  easily  handled  in  this  time,  multiplications  as  well  as 
multiply/accumulate  operations  are  fast  enough  when  performed  by  special 
devices  on  the  market,  but  standard  divisions  take  too  long,  qsproximately  2.2  y  s, 
for  the  fastest  devices.  Logic  arrays  to  handle  control  functions  are  sufficiently 
fast  to  be  used. 

m 

Size  and  cost  often  suffer  in  the  quest  for  speed.  For  example,  the  special 
multipliers  are  fairly  large  40-pin  devices  and  the  multiply/accumulate  packages 
ere  48  to  60-pin  devices,  depending  on  the  number  of  bits.  If  arithmetic  division 


were  performed,  12  dividers  would  have  to  be  multiplexed  resulting  in  a  large  size 
and  cost.  In  general,  optimal  speed,  size  and  cost  tradeoffs  must  be  carefully 
balanced. 

The  approach  taken  in  this  study  was  to  begin  by  reducing  size  and 
component  count  as  far  as  possible  within  the  constraint  of  operating  speed  and 
with  cost  as  a  consideration;  then,  with  this  near  minimal  size,  to  perform  card 
layouts,  and  finally,  on  those  cards  which  permit  expansion,  to  relax  size 
constraints  if  cost  savings  would  result.  Throughout  the  design  process,  the 
availability  and  cost  of  militarized  parts  has  also  been  a  prime  concern  to 
minimize  redesign  for  production. 

The  final  product  is  an  initial  design  which  is  as  small  and  inexpensive  as 
feasible,  yet  will  operate  at  the  fastest  rate  expected  to  be  required. 


Table  VI- 1  summarizes  the  considerations  discussed  in  this  section. 


SECTION  VII 

BRASSBOARD  DESIGN  AND  OPERATION 


For  the  brassboard  design,  a  single  enclosure,  individual  card  approach  has 
been  taken.  However,  card  designs  have  accounted  for  eventual  transition  to  the 
mother  board  concept. 


PIPELINE  CORRELATOR  CARDS 

The  pipeline  correlation  processors,  shown  in  figures  7.1  and  7.2,  have  fairly 
straightforward  operation.  The  required  inputs  consist  of  two  streams  of  video,  a 
pixel  clock  and  a  line  clock.  Data  manipulation  is  handled  in  random  access 
memory  (RAM)  and  operations  are  generally  handled  in  a  look-up  table  fashion 
with  programmable  read  only  memory  (PROM).  These  methods  are  not  only  the 
most  efficient,  but  offer  considerable  flexibility. 

Summations  over  a  window  illustrate  considerable  evolution  from  the  block 
diagram  given  in  section  II.  For  the  cross  scan  direction,  partial  results  are 
accumulated  for  one  line  and  output  at  the  end  of  a  line  for  final  processing. 

Summations  or  accumulations  rpay  be  performed  with  two  4-bit  adders  and  a 
latch  to  loop  the  sum  back  for  the  next  result.  Alternatively,  a 
multiply/accumulator  may  be  used.  Even  though  the  multiplication  function  may 
be  wasted,  these  devices  are  more  efficient  in  terms  of  space  requirements. 
Naturally,  they  are  also  more  expensive.  If,  howeyer,  the  first  method  results  in  . 
the  circuit  growing  too  large  to  be  fit  on  one  card,  the  cost  savings  may  be  lost  as 
shown  in  figure  7.3. 

The  scan  direction  offset  requires  more  complex  processing,  since  running 
summations  over  windows  covering  several  lines  and  being  output  at  a  pixel  rate 
during  a  line  are  needed.  The  shift  register  architecture  suggested  in  the  block 
diagrams  in  section  II  is  not  a  good  hardware  approach.  Shift  registers  are  neither 
available  in  arbitrary  lengths  nor  flexible.  RAM's  are  considerably  better  choices, 
even  if  they  may  provide  for  more  memory  than  required.  Thus,  instead  of 
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shifting  data  through  shift  registers,  the  data  is  stored  in  one  location  and  the 
address  pointer  loops  through  locations  as  shown  in  figure  7.4.  The  top  diagram 
showing  a  running  sum  of  several  pixels  indicates  a  subtraction  and  addition  to  the 
accumulation.  This  is  clearly  not  feasible.  The  proper  approach  is  a  subtraction 
followed  by  accumulation,  as  shown  in  figure  7.5. 

Subtraction  itself  is  not  performed  as  easily  as  indicated  thus  far.  A  typical 
approach  might  be  to  perform  a  2's  complement  on  the  subtrahend  and  then  add  by 
standard  methods.  Not  only  does  the  complement  operation  require  an  addition 
with  the  number  I ,  but  sign  bits  necessitate  the  use  of  more  components  and  the 
entire  design  becomes  quite  involved.  A  somewhat  cleaner  approach  which  handles 
positive  and  negative  numbers  in  2's  complement  form  is  shown  in  figure  7.6. 
Curiously,  it  is  the  minuend  rather  than  the  subtrahend  which  is  inverted.  The 
final  result  is  then  inverted  (unless  the  result  is  to  be  the  minuend  for  a  second 
subtraction),  to  arrive  at  the  correct  answer  in  2's  complement  notation.  As 
illustrated,  full  resulution  is  achieved  with  a  simple  logical  "OR"  gate.  Less  I -bit 
resolution  consists  of  dropping  an  LSB  and  extending  the  sign  bit  into  the  MSB  pins 
on  the  adders.  Also  shown  in  figure  7.6  is  the  simplest  approach,  a  PROM  lookup 
table. 

Finally,  the  summing  over  block  windows,  shown  in  figure  7.4,  still  shows  an 
adding  tree  to  arrive  at  a  final  sum. ’This  is  cumbersome  at  best.  A  preferred 
approach  is  more  like  the  first  accumulation  shown  in  figure  7.3,  except  that  the 
latch  becomes  a  RAM  at  least  I -line  long  to  store  summations  by  pixel.  The  adder 
tree  and  the  accumulation  approaches  are  compared  in  figure  7.7.  The  tree 
approach  requires  n-l  lines  of  memory,  but  a  separate  device  for  each  line,  while 
the  accumulation  technique  require  n+l  lines,  but  only  two  devices  are  needed. 


FRAME  STORE 

The  design  of  the  frame  store  portion  of  the  stabilization  module  is  far  less 
involved  than  for  the  pipeline  correlator.  This  is  shown  in  figures  7.8  and  7.9. 
Modularity  demands  that  the  memory  plane  be  segmented  into  60  pixel  groups, 
while  commonality  dictates  that  this  architecture  be  identical  to  that  used  for  the 
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special  function  module.  It  has  already  been  shown  that  a  4  x  4  array  of  60-pixel 
by  480-line  segments  is  desirable  for  the  latter,  and  there  is  no  reason  to  differ  in 
the  basic  frame  store.  Thjee  fields  of  storage  are  required  for  proper  operation: 
one  for  the  writing  of  current  data,  one  to  hold  the  previous  field  until  it  can  be 
correlated  with  the  next  field,  and  one  to  read  for  correlation  with  the  current 
field.  It  is  not  possible  to  combine  the  first  and  third  memory  fields  into  one 
because  offset  conditions  could  be  such  that  memory  locations  would  be  written 
with  current  data  before  they  are  required  to  be  read  for  correlation.  Figure  7.10 
illustrates  the  sequence  of  operation.  Note  that  alternate  fields  are  written  into 
memory  in  reverse  order.  There  are  several  benefits  to  this  method.  First, 
addressing  control  is  simplified  since  the  counters  need  only  count  in  opposite 
direction  rather  than  be  reloaded  with  zero.  Second,  no  special  control  is  needed 
for  bootstrapping.  Bootstrapping  merely  consists  of  holding  the  field  control  bit 
for  one  field.  Third,  scene  locations  occupy  nearly  the  same  addresses  in  each 
field.  This  is  important  to  simplify  the  window  locating  and  marking  functions. 

In  summary,  while  the  relationship  between  the  actual  hardware  design  and 
block  diagram  form  discussed  in  section  2  may  not  be  immediately  obvious,  they 
are  functionally  equivalent.  The  hardware  design  is  achieved  with  the  fewest 
number  of  available  components,  while  the  block  diagram  was  "designed"  to 
explain  the  operation. 

It  has  been  mentioned  that  this  memory  operates  in  an  unstabilized  mode  so 
that  current  video  is  written  in  normal  order,  but  to  track  offsets,  the  stored 
image  is  retrieved  with  offsets  added  to  the  addresses.  There  are  several 
implications  associated  with  this  procedure.  The  ^irst,  already  discussed,  is  that 
reading  and  writing  may  overlap  so  that  three  memory  planes  are  require.  Second, 
addresses  will  have  to  be  checked  for  validity  after  offsets  have  been  added.  And 
third,  when  an  image  is  undergoing  either  skew  or  rotation,  the  scan  direction 
offset  may  be  sawtooth  waveform.  Because  of  this  possiblity,  the  scon  direction 
offset  of  the  first  window  in  each  line  must  be  latched  and  reloaded  at  the  end  of 
a  line.  This  allows  the  system  to  track  skew  along  a  line  without  losing  reference 
from  the  fnd  of  one  line  to  the  begirmng  of  the  next  as  shown  in  figure  7.11. 


BASIC  FRAME  STORE  -  MEMORY  OPERATION 


FIGURE  7.10 
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The  most  obvious  method  to  provide  address  with  offset  is  to  add  or  subtract 
a  count  from  the  base  address.  As  previously  pointed  out,  subtraction  is  not  as 
easy  a  procedure  as  it  might  seem.  Therefore,  a  bias  offset  would  seem 
appropriate  to  avoid  adding  negative  numbers.  However,  this  will  complicate 
address  decoding.  Furthermore,  this  method  is  somewhat  inflexible  in  allowing  a 
change  in  offset  of  only  one  location  per  update.  An  improved  method  is  achieved 
by  decoupling  the  nonoffset  addresses  and  the  offset  addresses.  Being  decoupled, 
care  must  be  taken  to  assure  synchronization  between  the  two  address  generators. 

The  alternate  method,  shown  in  figure  7.12,  uses  a  multiply/accumulator.  At 
one  input,  a  constant  "I"  is  presented  to  the  device  while  the  other  input  is  an 
integer  which  is  determined  by  the  offset  data.  For  no  offset  change,  a  I  is 
presented  and  the  accumulator  increments  one  address  position  as  would  be 
expected.  To  effect  a  +1  offset  change,  a  2  is  presented  and  the  address  skips  a 
location.  Similarly,  a  -I  offset  change  causes  a  0  to  be  presented  and  the  address 
holds  for  a  count.  This  device  allows  very  sophisticated  and  flexible  control 
functions.  For  example,  address  resetting  at  the  end  of  a  line  may  be  easily 
implemented  by  presenting  a  -60  to  one  input  and  a  1 ,  2  or  3  to  the  other, 
depending  on  the  sensor  line  length.  Alternatively,  a  slight  skip  at  the  end  of  a 
line  would  allow  line  and  pixel  address  control  in  a  single  device. 

The  rest  of  the  design  and  operation  of  the  frame  store  is  relatively 
straightforward,  involving  chip  selection,  address  decoding  and  read/write  control. 

MTI/RANGWG  CARD 

The  MTI/range  function  consists  of  calculating  offsets  for  small  square 
windows  and  providing  direct  output,  binary  output  of  offsets  greater  than  a 
threshold  or  a  scaled  output.  Since  the  scan  direction  offset  is  already  calculated 
over  the  desired  window  by  the  pipeline  correlator,  these  results  can  be  used 
directly.  For  the  cross  scan  direction,  however,  the  window  of  results  from  the 
pipeline  correlator  is  a  single  line.  Therefore,  the  partial  results  must  be  tapped 
and  summed  over  a  square  window  before  final  processing  to  derive  an  offset. 

For  motion  outputs,  the  average  motion  should  be  calculated  and  subtracted 
from  eoch  result,  so  that  the  results  do  not  include  sensor  motion.  Further,  it  may 
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FIGURE  7.12 


be  desirable  to  calculate  the  magnitude  of  motion  from  the  two-axis  results.  For 
ranging,  on  the  other  hand,  normalization  is  not  desired  nor  are  the  axis  results 
combined.  Range  is  determined  solely  by  the  cross  scan  offsets,  while  the  scan 
offsets  are  used  to  synchronize  the  scan  mirrors  of  the  two  sensors. 

For  range  or  motion,  there  are  three  options:  thresholding  for  "within  range" 
or  MTI  outputs,  unsealed  outputs  for  relative  maps  or  scaled  output  for  absolute 
results.  Figure  7.13  illustrates  the  card  design  in  block  diagram.  Except  for  the 
motion  magnitude,  which  is  performed  as  a  look-up  table  in  PROM,  the  design 
considerations  and  operation  are  virtually  identical  with  the  basic  pipeline 
correlation  function. 


SPECIAL  FUNCTION  MODULE 

The  special  function  module  is  very  similar  to  the  correlation  frame  store. 
The  memory  architecture  is  identical  except  that  only  two  planes  are  required, 
one  for  each  field  since  there  is  no  possiblity  of  contention  between  reading  and 
writing.  While  one  field  is  being  stored  in  one  plane,  the  other  field  is  being  read 
from  the  second  plane.  In  this  memory,  when  the  image  is  to  be  passively 
stabilized,  the  current  video  is  written  to  addresses  with  offset,  while  reading  is 
performed  in  normal  order  without  offset.  This  is  the  reverse  of  the  procedure  for 
the  correlation  frame  store.  . 

It  is  in  the  control  where  the  differences  between  this  special  function 
module  and  the  correlation  frame  store  lie.  The  handling  of  offsets  is  the  same  as 
described  for  the  latter  -  that  is,  using  multiply/accumulators,  but  here  the  power 
of  this  method  is  used  more  fully.  For  example,  when  the  image  is  being 
dynamically  sampled,  the  address  counting,  the  offsets  and  the  window  position 
must  be  scaled  by  the  scaling  ratio.  This  is  easily  implemented  by  including  scale 
ratio  as  an  input  to  the  FPLA,  so  that  the  multiplicand  reflects  this 
consideration.  There  has  been  some  rearrangement  of  how  input  terms  to  the 
multiply/accumulator  are  handled,  but  the  general  operation  is  the  same  for  both 
frame  stores. 


As  previously  discussed,  it  is  necessary  to  accumulate  offset  results  to 
maintain  reference  between  the  unstabilized  correlator  frame  store  and  the 
stabilized  special  function  frame  store.  The  accumulation  is  quite  similar  to  the 
window  partial  results  accumulations  used  on  the  correlation  cards.  A  RAM  is 
used  to  store  pixel  results  along  a  line,  or  line  results  along  a  frame,  and  fresh 
offset  data  is  summed  to  the  sum  of  prior  data.  When  initializing  a  fresh  scene 
lock-on,  the  memory  is  effectively  cleared  by  adding  in  zeroes  for  the  first  frame. 

The  controller  design  is  shown  in  figure  7.14. 

Another  function  which  is  found  in  this  module  is  scene  integration.  There 
are  several  possible  ways  to  implement  this  function.  This  integration  can  be 
written  in  two  forms: 

S'  =  X  x  C  +  (I  -  X)  S 
or:  S'  -  S  +  X  (C  -  S) 

where: 

C  =  current  video 

S  *  stored  video 

S'  =  result  of  integration  to  be  stored 

Figure  7.15  illustrates  the  implementations  which  were  considered.  Method 
I  uses  several  components,  including  two  high-speed  multipliers,  which  use 
disproportionate  amount  of  real  estate  for  the  task.  Method  2  is  somewhat  less 
expensive  in  cost  and  real  estate,  but  still  uses  a  high-speed  multiplier  and  the 
subtraction  function,  as  has  been  previously  stated,  is  not  a  simple  operation. 

Both  of  the  first  two  methods  allow  virtually  any  value  of  persistence.  The  third 
method  limits  persistence  values  to  eight  possible  values.  Prior  discussion  has 
shown  that  eight  integration  ratios  are  quite  sufficient.  Now,  a  simple  look-up 
table  suffices  to  perform  the  (l-X)S  operation,  while  the  X  x  C  operation  is  easily 
coded  in  FPLA.  The  third  approach  results  in  the  least  expensive,  fastest  and 
smallest  configuration,  a  rare  occurrence. 
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CONTROL  PANEL 


For  the  proposed  brassboard  design,  a  single  control  panel  will  be  used 
(shown  separately  as  figures  7.16  and  7.17).  The  control  panel  allows  for  entering 
pixel  size,  focal  length  and  baseline  or  range  data  via  a  pushbutton  keyboard. 
Information  such  as  threshold,  absolute  velocity  or  range,  and  error  conditions  are 
displayed  on  an  alphanumeric  display.  The  window  position  is  controlled  by  a  set 
of  four  pushbutton  switches  arranged  as  an  effective  joystick.  The  switches  are 
used  instead  of  a  joystick  for  their  lower  cost  and  more  rugged  construction.  They 
also  eliminate  the  danger  of  having  a  protruding  lever  on  the  control  panel.  Other 
controls  include  scale,  (lx,  2x,  4x),  stability  and  integration  persistance. 

The  heart  of  the  control  panel  is  designed  with  an  8251  single-chip 
microprocessor.  The  microprocessor  provides  on-chip  ROM,  RAM,  a  serial  port 
and  parallel  ports.  The  version  of  this  microprocessor  that  will  be  used  for  the 
brassboard  design  is  the  8031 ,  which  includes  all  the  above  features  except  on-chip 
ROM,  which  is  supported  by  a  2716  EPROM.  This  redues  the  cost  for  the  initial 
design.  Communications  to  and  from  the  main  electronic  box  are  accomplished  by 
an  EIA  Standard  RS-232C  interface,  which  requires  only  three  wires.  This  feature 
reduces  the  cable  size  requirements.  The  MIL-STD- 1 553B  interface  could  be 
substituted,  but  the  additional  cost  due  to  the  increased  parts  count  and  software 
design  time  is  not  considered  cost-effective  for  a  brassboard  system.  The 
communication  interface  in  the  main  electronics  box  is  also  an  8251  single-chip 
microprocessor  that  supports  ancillary  functions.  These  functions  include  non- 
critical  mathematical  computation,  system  configuation  and  system 
initialization.  Also,  errors  can  be  sensed  and  displayed  on  the  control  panel,  as 
shown  in  figure  7.18. 

The  Intel  827?  keyboard/display  interface  is  used  for  controlling  the 
keyboard  and  display.  This  eliminates  the  need  for  the  microprocessor  to  "poll" 
the  keyboard  (look  for  a  key  pressed)  and  multiplex  the  displays.  It  also  allows 
greater  flexibility  in  the  software  design. 


MTI/RANGE  CONTROL 


FIGURE  7.16 


STABILIZATION  CONTROL 


FIGURE  7. 17 
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SECTION  VIII 

PRODUCTION  DESIGN  CONSIDERATIONS 


A  modular  design  approach  will  be  used  for  production  models.  Four  mother 
boards  will  be  desigied  to  accomdate  the  various  options  in  a  reconfigurable 
manner. 

When  fully  configured,  the  system  supports  the  pipeline  correlator,  frame 
store,  MTl/range  and  communication  interface  and  spacial  functions.  Internal  DIP 
switches  allow  for  selecting  line  width  and  pixel  count.  Each  electronics 
enclosure  can  have  its  own  power  supply  or  a  common  supply  can  be  used. 

Two  control  panels  are  used  for  the  full  system.  This  design  allows  for 
separate  control  panels  for  special  function  and  MTl/range  control.  These  panels 
are  designed  to  meet  MIL-STD-25212.  The  production  units  will  use  the  8251 
single-chip  microprocessor  to  save  board  space  and  reduce  total  package  count.  A 
switch  located  inside  the  main  electronics  box  allows  the  system  power  to  be 
controlled  at  the  control  panel. 

Each  option  is  configured  with  a  frame  store  (memory)  for  60  pixels. 

Additional  frame  storage  is  available  in  60-pixel  increments,  with  a  maximum 

* 

number  of  1 80  pixels. 

Option  I  is  a  basic  stabilization  system  which  uses  only  the  main  electronics 
enclosure  with  the  pipeline  correction,  memory  and  control  boards  installed. 

I 

Option  2  adds  the  MTI  function  to  the  basic  system.  With  this  optign,  the 
MTl/range  control  panel  is  included  along  with  the  MTl/range  board  in  the  main 
electronics  enclosure.  Option  3  provides  the  stabilization  and  ranging  functions. 
This  option  uses  th‘,  main  electronics  enclosure  with  the  pipeline  corelator, 
MTl/range  and  master  control  boards,  however  the  frame  storage  is  not 
required.  The  MTl/range  control  panel  is  also  used.  Option  4  features 
stabilization  and  special  functions.  Both  control  panels  are  used  for  this  option. 
Option  S  contains  all  the  features  (MTl/range,  stabilization,  special  functions)  and 
is  considered  the  full  system  configuration. 


The  only  differences  between  the  brassboard  system  and  production  design 
are  the  operation  of  functions  onto  dedicated  mother  boards,  the  use  of  printed 
circuit  boards,  and  the  use  of  only  military-qualified  versions  of  components. 


SECTION  IX 

CONCLUSIONS  AND  RECOMMENDATIONS 


This  study  and  preliminary  design  has  shown  that  correlation  processing  can 
be  used  for  FLIR  sensor  stabilization  with  better  performance  than  gyro 
transducers,  while  providing  many  other  useful  and  powerful  functions.  In 
particular,  precision  correlation  has  been  shown,  by  a  computer  simulation 
program,  to  measure  image  displacements  with  errors  under  1/20  pixel.  This 
accuracy,  combined  with  the  high-speed  pipeline  implementation,  allows  new 
types  of  information  to  be  displayed  as  images:  motion  maps,  range  maps 
obtained  without  active  sources,  and  dynamically-sampled  imagery  with  resolution 
up  to  3.5  times  the  Nyquist  limit. 

The  preliminary  design  effort  has  developed  a  modular  approach  which 
allows  users  to  select  the  particular  functions  required  for  their  application  and  to 
add  or  delete  functions  as  desired  without  modification  to  other  portions  of  the 
system.  Considerable  attention  has  also  been  paid  to  size,  power  and  cost  in  the 
preliminary  design,  as  well  as  future  transition  to  a  fully  militarized  production 
system.  The  result  is  a  versatile,  affordable  menu  of  options. 

It  is  recommended  that  a  brassboard  system  be  fabricated  to  demonstrate 
the  capabilities  described.  The  system  will  be  used  in  conjunction  with  the 
common  module  digital  scan  converter.  Video  input  to  the  system  will  be  taken 
from  the  digitized,  multiplexed  output  and  results  will  be  displayed  via  the 
reformatter.  It  is  recommended  that  the  brassboard  system  be  configured  for 
120-element  array  sensors.  The  system  will  demonstrate  passive  stabilization, 
motion  mapping,  scaling,  frame  integration  and  dynamic  sampling.  With  two 
boresighted  synchronized  FLIR  sensors,  the  system  can  also  demonstrate  passive 
ranging  and  display  of  a  range  map. 


APPENDIX  A 

PRECISION  CORRELATION 


INTRODUCTION 

Correlation  can  be  used  for  two  purposes: 

•  To  determine  degree  of  match  by  determining  the  height  of  a 
correlation  peak 

•  To  measure  the  relative  alignment  between  two  signals  by  measuring 
the  position  of  the  correlation  peak 

It  is  the  latter  purpose  which  CA1  has  addressed  with  precision  correlation.  There 
are  several  applications  for  precision  correlation.  These  include  stabilization, 
tracking,  MTI,  passive  ranging,  navigation  and  terrain  following. 


CORRELATION  THEORY 

The  applications  of  precision  correlation  depend  upon  the  ability  to 
determine  the  relative  shift  between  two  images.  Succinctly,  it  is  required  to 
measure  how  far  one  image  must  be  displaced  to  match  the  second  image.  The 
first  task  is  then  to  define  match  in  mathematical  terms. 

Since  brightness  is  the  only  information  with  which  to  work,  match  must  be 
defined  as  occurring  when  each  sample  in  one  image  plane  has  the  same  amplitude 
as  its  corresponding  sample  in  the  second  image  plane.  Thus, 

P|  (i»j)  =  P2  G  +  A  i>  j  +  Aj)  for  all  i  and  j  where  Ai  and  Aj 
define  the  shift  required  to  make  the  match  using  Wne  coordinate  system.  In 
reality,  noise,  lens  mismatch,  specular  reflections  and  the  like  prevent  the 
difference  from  ever  being  precisely  zero  over  all  i  and  j;  hence,  the  need  to  find 
the  best  match  or  fit.  One  common  measure  of  fit,  given  numerous  samples,  is  to 
take  the  least  squares  approach  in  which  parameters  are  varied  until  the  sum  of 
the  squares  of  the  differences  shows  a  minimum.  In  this  case,  the  parameters  to 
be  varied  are  A  i  and  A  j .  There  are,  of  course,  other  measures  of  fit  such  as 
minimizing  the  sum  of  the  absolute  values  of  the  differences.  In  this  vein,  it  can 


be  generalized  to  minimize  j  [p,  (i,j)  -  P2  (i  +  Ai»  j  -►  Aj)J  ^  j  ^ 
where,  for  example,  Q  =  I  and  1/2,  respectively  for  the  two  cases  mentioned.  It 
can  be  shown  that  the  least  squares  approach  is  mathematically  equivalent  to  the 
statement  that  the  mean  of  statistical  sample  is  the  best  estimator  of  that 
sample,  while  the  least  absolute  value  is  equivalent  to  the  statement  that  the 
median  is  the  best  estimator.  Much  early  work  in  correlation  matching  was 
devoted  to  determining  which  statement  is  the  more  valid  for  imagery,  with  the 
conclusion  that  the  least  squares  approach  yields  a  less  noisy  correlation  curve. 

Deciding  a  priori  that  average  brightness  and  gain  differences  between  the 
two  images  are  to  be  ignored  in  determining  match,  the  sample  sets  are 
normalized  by  subtracting  the  mean  and  dividing  by  the  standard  deviation.  For 
the  least  squares  approach,  a  match  is  defined  by  locating  the  minimum  of: 


P,  (i,j)-P. 


P2  (i+Ai,  j+Aj)-P2 


Expanding: 


Min  ( A i , Aj ) :  £ 


P9 


P?  (i+Ai,  j+Aj)-P2 

+  £  -£ - - - ~ 
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L  2 


P|  (i,j)-P, 

■2  Op 


P2  (i+Ai,  j+Aj)-P2 


Simplifying: 

Min  (Ai,Aj) 


'2  '  2  11 
E  P,  +  £  P2  -2  £  P,  P2 


The  first  two  summations  both  equal  I  by  definition  of  mean  and  standard 
deviatic/i.  What  remains  is  finding  the  minimum  of  2  £|-e  Pj  P2J  ♦  This  is 
equivalent  to  finding  the  maximum  of  £  P  |  P2  which  is  recognized  as  the  peak 
correlation  coefficient.  Thus,  correlation  is  nothing  more  than  a  statement  that 


match  is  defined  by  minimizing  the  least  squares  difference  between  two  signals 
as  a  function  of  relative  shift,  with  the  implication  that  the  mean  of  the  sample- 
by-sample  differences  is  the  best  estimator  of  the  degree  of  match. 

By  plotting  the  correlation  product  as  a  function  of  image  shift,  a 
correlation  curve  is  generated.  Because  the  images  are,  in  this  case,  focused  on 
discrete  element  sensors,  the  images  are  sampled  in  a  discrete  regular  pattern. 
With  the  images  sampled  discretely,  the  correlation  curve  is  only  defined  at 
discrete  intervals.  Nevertheless,  these  correlation  samples  do  in  fact  represent  a 
continuous  curve.  Furthermore,  this  curve  shifts  precisely  with  the  relative  shift 
of  the  two  images. 


PRECISION  CORRELATION 

Clearly,  the  image  shift  can  be  determined  in  a  gross  sense  merely  by 
locating  the  peak  correlation  sample.  However,  to  measure  a  shift  smaller  than 
the  sample  spacing  (i.e.,  pixel  spacing),  the  curve  must  be  interpolated.  In 
essence,  this  means  assuming  a  mathematical  shape  for  the  curve  and  calculating 
the  necessary  constants  based  upon  the  samples  obtained  by  autocorrelation  (the 
correlation  of  one  image  with  itself)  which  is  known  to  be  centered  at  the  origin. 
Once  the  constants  are  evaluated,  'he  fractional  pixel  shift  can  be  determined  by 
applying  these  values  to  the  samples  of  the  correlation  curve  obtained  by  cross 
correlation  of  the  two  different  images.  This  procedure  assumes  that  the  curve 
maintains  its  shape  as  it  shifts  with  the  relative  image  shift.  In  practice,  the 
mathematics  can  be  arranged  so  that  the  constants  never  need  to  be  evaluated 
explicity,  but  the  cross  correlation  sample  values  and  autocorrelation  sample 

I 

values  are  combined  in  one  equation  to  directly  yield  the  fractional  pixel 
displacement. 

The  actual  shape  assumed  for  the  curve  can  be  gaussian,  exponential, 
polynominal,  spline,  trigonometric  or  any  one  of  a  wide  variety  with  various  levels 
of  generalization.  The  performance  of  any  algorithm  making  an  assumption  of 
curve  shape  must  be  assessed  both  in  terms  of  how  well  the  assumed  shape  fits 
actual  shapes  generated  with  real  imagery,  ease  of  implementation  and  processing 
speed.  A  mismatch  between  the  assumed  shape  and  actual  shapes  of  correlation 
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curves  obviously  produces  errors  in  the  calculated  interpolation.  Furthermore, 
care  must  be  taken  in  deciding  which  samples  of  the  curve  to  use  since  the  curve 
is  quickly  dominated  by  noise  effects  immediately  around  the  peak,  yet  the  peak 
itself  on  the  autocorrelation  curve  has  correlated  noise  included  which  will  not  be 
found  in  cross  correlation.  Finally,  filtering  of  the  scene  can  help  or  hurt  the 
correlation  accuracy  by  shaping  the  correlation  curve  to  some  extent. 

Even  if  real  correlation  curves  could  be  well  described  in  practice  by  some 
mathematical  shape,  two  error  sources  remain.  The  first  is  noise,  which  is  always 
present  to  degrade  results.  The  second  is  fine  scene  structure.  Spatial 
frequencies  about  and  greater  than  Nyquist  frequency  are  strongly  phase 
dependent  in  their  effects  on  sample  amplitudes.  This  phase  dependence  can 
produce  fluctuations  in  the  correlation  curve  which  are  not  detectable  solely  by 
autocorrelation  samples.  Filtering  can  be  a  solution  this  problem  only  so  long  as 
the  filter  does  not  greatly  reduce  the  slightly  lower  spatial  frequencies  which  help 
to  sharpen  the  correlation  curve  and  allow  the  accurate  measurement  of  small 
displacements. 

CAI  has  expended  considerable  effort  in  determining  the  best  assumed  shape 
consistent  with  ease  of  implementation  and  speed  of  operation.  The  findings  have 
shown  that  without  filtering,  RMS  correlation  error  of  I /50th  of  a  pixel  can  be 
achieved.  Further  improvement  maybe  obtained  with  noise  and/or  scene 
filtering. 

EXPERIMENTAL  RESULTS 

The  accuracy  of  the  algorithm  is  dependent  on  the  signal-to-noise  ratio 
inherent  in  the  signal  and  on  the  number  of  signal  samples  used  in  the 
calculation.  Figures  A.  I  and  A.2  show  computer  and  analytic  results  (which  have 
been  verified  on  a  physical  demonstration  unit)  of  rms  error  versus  signal-to-noise 
ratio  and  number  of  samples.  Further  accuracy  is  obtained  by  filtering  the 
correlation  results  with  a  low-pass  filter  as  shown  in  figure  A.3. 


RMS  ERROR 


Generally,  a  basic  correlation  accuracy  of  I /50th  sample  spacing  is  used  for 
calculating  expected  performance.  In  many  instances,  however,  as  will  be  shown, 
the  application  may  lend  itself  to  filtering  to  improve  the  basic  results. 


USES  AND  LIMITATIONS 

The  most  obvious  immediate  applications  for  precision  image  registration 
are  tracking  and  its  related  function,  stabilization.  The  high  precision  afforded  is 
critical  for  high-resolution  imaging,  which  is  important  in  tracking  for  fire  conrol 
and  useful  in  general  tracking.  A  less  obvious  application  is  passive  ranging. 
Ranging  is  identical  to  tracking,  except  that  the  images  displaced  in  time  used  for 
tracking  are  replaced  by  images  displaced  in  space  to  obtain  stereo  images  as 
shown  in  figure  A.4. 

The  formula  for  range,  R,  given  two  parallel  optical  systems  separated  by  a 
distance,  B,  is  (noting  similar  triangles): 


where: 

F  =  focal  length  of  he  optics 

D  =  relative  displacement  between  the  object  images 

Assuming  B  and  F  are  constant,  the  error  is  determining  range  is  limited  by  the 
error  in  measuring  D: 

RF  r2 

dR  =  -  dD  =  -  —p~  dD 

Since  the  sample  spacing  on  the  correlation  curve  is  equivalent  to  the  sample 
spacing  of  the  signal,  that  is  a  resolution  element  or  pixel  spacing  at  the  image 
plane: 

R2 

dR  =  BF  Pe 

where: 

P  =  spacing  between  pixels 
e=  algorithm  accuracy 
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In  contrast  to  active  rangefinding  with  fairly  constant  range  error,  this 
method  has  a  range-squared  dependence.  This  is  a  disadvantage  in  the  sense  that 
to  accommodate  greater  ranges,  the  physical  package  must  expand  for  this  passive 
technique;  while  an  active  technique  requires  only  increased  power  to  overcome 
weakening  due  to  beam  spread.  Currently,  power  seems  to  be  a  less  expensive 
commodity.  Not  only  is  increased  package  size  undersirable  in  itself,  but 
maintaining  rigid  alignment  becomes  impossible  and  other  complexities  must  be 
incorporated  to  compensate.  On  the  other  hand,  there  is  no  range  gating  so  no 
ambiguity  of  ranges  con  occur.  It  is  also  generally  true  that  the  nearer  the  object, 
the  higher  the  accuracy  desired,  or  at  least  the  more  critical  is  that  accuracy. 
Stereoscopic  ranging  does  not  waste  data  bits  producing  highly-accurate  results 
for  far  objects  to  achieve  high  accuracy  on  near  objects.  There  is  also  a  practical 
near-range  limit  due  to  excessive  image  shift.  Even  if  the  taget  image  has  not 
moved  off  the  sensor  format  at  very  close  ranges,  the  processing  cannot  afford  to 
generate  a  correlation  curve  several  hundreds  of  samples  long  to  locate  a  match. 

Applications  of  passive  rangefinding  include:  fire  control,  altimeters, 
TERCOM,  terrrain  avoidance,  target  detection  and  possibly  target  recognition. 

For  all  but  the  first  two  applications,  it  is  desirable  to  map  correlation  results 
over  the  image  format.  However,  using  standard  techniques  for  calculating  a 
correlation  curve,  locating  a  peak  sample  and  interpolating  the  position  of  the 
curve  peak,  only  a  few  results  per  frame  could  be  computed.  CAI  has  found  that 
its  precision  correlation  algorithm  lands  itself  to  pipeline  implementation, 
allowing  the  kind  of  high-speed  mapping  desired. 


APPENDIX  B 

PIPELINE  PRECISION  CORRELATION 


OPERATION 

Pipeline  precision  correlation  (PPC)  is  a  specialized  form  of  CAI's  precision 
correlation  algorithm  wherein  each  arithmetic  operation  is  separately 
implemented  in  hardware.  Two  serial  streams  of  video  data  are  input  and 
processed  to  produce  two  serial  streams  of  partial  results  on  a  pixel-by-pixel 
basis.  These  partial  results  are  summed  over  a  desired  correlation  "window."  The 
sums  are  then  processed  in  a  final  stage  to  produce,  on  a  pixel  shift-by-pixel  shift 
basis  if  desired,  the  calculated  misregistration  for  each  window.  If  suitably 
designed,  the  dimensions  of  the  window  can  be  dynamically  adjusted. 

Conceptually,  incoming  video  information  is  simply  transformed  into 
misregistration  results.  The  signal  may  be  interlaced  or  noninterlaced  and  the 
results  will  be  similar  in  format.  The  non  interlaced  format  is  nonetheless 
preferred  because  the  ability  to  correlate  on  small  objects  is  improved. 


-VIDEO  LIKE"  VERSUS  BLOCK  RESULTS 

As  shown  in  figure  B.I,  results  Ore  calculated  as  running  results  for  all 
overlapping  windows  mXn.  Results  could  also  be  calculated  for  adjacent  windows 
only,  but  this  would  be  only  minimally  easier  to  accomplish.  Running  "video  like" 
results  is  to  an  extent,  but  not  entirely,  oversampling.  Not  only  is  aliasing 
reduced  by  oversampling,  but  the  larger  number  of  results  allow  smoothing,  which 
improves  the  overall  accuracy. 

A  particularly  attractive  implementation  occurs  if  either  the  desired  window 
width  equals  the  scan  line  or  the  window  height  is  a  single  scan  line.  In  either 
case,  the  need  for  an  array  of  shift  registers  or  their  equivalent  is  unnecessary. 


RELATIVE  VERSUS  ABSOLUTE  RESULTS 


With  time  sequential  applications  such  as  tracking,  absolute  motion  can  be 
determined  only  insofar  as  the  motion  of  the  sensor  is  known.*  It  is  rarely 
necessary  that  absolute  motion  of  some  object  be  determined  in  dynamic 
conditions.  In  ranging,  however,  absolute  values  are  often  desired  even  for 
dynamic  conditions.  The  implication  of  absolute  range  with  a  steroscopic  sensor  is 
that  the  base'ine  separation  and  angular  alignment  be  precisely  known.  Typically, 
ranging  appnations  requiring  accuracies  of  I  /50th  sample  spacing  are  equivalent 
to  determining  parallax  to  a  few  microradians.  Therefore,  the  relative  angular 
alignment  must  be  held  or  measured  to  such  tolerances. 

In  terms  of  mechanical  stability,  this  is  a  difficult  proposition  at  best.  An 
alternative  to  exceedingly  rigid  mechanical  tolerarices  is  incorporation  of  some 
form  of  autocalibration  in  the  optomechancial  system.  This  may  be  some  means 
of  allowing  one  sensor  to  "see"  the  other  via  most  of  the  optical  path  within  the 
physical  limits  of  the  system,  or  for  each  to  "see"  a  reticle  or  pair  of  reticles,  as 
is  used  in  highly-precise  visual  rangefinders  shown  in  figure  B.2.  The  reticle 
pattern  is  separated  from  the  scene  both  spectrally  and  in  field  position,  so  that 
each  may  be  correlated  simultaneously  to  obtain  range  regardless  of  any  optical 
wedge  within  the  system.  Because  the  reticles  can  be  an  optimum  pattern,  well 
illuminated  with  good  contrast  and  of  large  extent,  the  correlation  precision 
achieved  is  expected  to  well  exceed  The  typical  I /50th  sample  spacing  accuracy  of 
non-ideal  objects. 


Because  of  the  increased  physical  and  processing  complexity  required  to 
ensure  accurate  absolute  ranging,  it  is  appropriate*  to  examine  which  applications 
require  it.  Certainly  fire  control  appliations  do  require  absolute  range  to  solve 
the  gun  directing  equation.  A  TERCOM  application  conversely  can  use  relative 
data,  since  it  is  only  the  terrain  contour  which  is  of  interest.  Care  must  be  taken 
when  matching  to  normalize  not  only  for  offsets  in  absolute  range  due  to 
uncertainty  of  altitude,  but  for  offsets  in  reciprocal  range  due  to  misalignment  of 
optical  channels. 


*  With  a  stationary  background,  the  sensor  motion  is  easily  calculated  as  the 
reverse  of  the  apparent  background  motion. 
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For  terrain  avoidance,  it  may  be  initially  assumed  that  absolute  ranging  is 
needed.  If  it  is  required,  there  are  several  options.  First,  the  necessary 
autocalibration  could  be  included  in  the  system  design.  Second,  a  separate,  active 
rangefinder  could  provide  the  absolute  range  to  a  point  in  the  scene  which  would 
allow  scaling  the  entire  relative  map.  In  this  way,  use  of  an  active  sensor  can  be 
kept  to  a  minimum  with  an  occasional  random  flash.  Third,  the  absolute  range 
could  be  roughly  calculated  by  measuring  the  angular  rate  of  some  off-axis  portion 
of  the  scene.  The  formula  for  range  parallel  to  the  line  of  flight  (R|_qF^  's  9iven 
by: 
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where: 
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V  =  platform  velocity 

Pj  =  image  location  of  an  object  at  time  T 

T  =  current  time 

T'  =  time  at  which  the  image  location  of  on  object  was  previously  noted 


There  are  several  potential  error  sources  in  the  latter  option  of  range 
calculation,  but  over  a  period  of  time  and  over  several  image  points,  the  accuracy 
should  be  acceptable.  In  fact,  if  a  human  is  the  sole  user  of  the  relative  range 
map,  it  is  possible  that  observation  of  the  angular  rates  on  display  would  be 
sufficient  to  estimate  absolute  range  to  objects  and  take  appropriate  action. 
These  three  options  are  not  mutually  exclusive  and  may  be  used  in  combination  to 
obtain  greater  accuracy  than  any  one  alone.  With  the  cues  given  by  angular  rate, 
it  is  not  clear  that  the  explicit  calculation  of  absolute  range  is  definitely  required 
for  terrian  avoidance. 


LIMITATIONS 

Because  PPC  implements  only  the  correlation  curve  interpolation,  the 
correlation  curve  itself  is  never  calculated.  This  fact  has  several  implications. 
First,  there  is  no  checking  for  degree  of  match.  It  must  then  be  assumed  that  the 
two  input  video  signals  represent  a  matched  pair.  Similarly,  the  peak  sample  is 
not  located  explicitly.  It  is  assumed  that  the  peak  sample  is  within  one  or  two 


sample  spacings  of  some  nominal  value.  PPC  calculates  only  the  interpolated 
peak  position.  The  results  are  valid  strictly  if  the  actual  misregistration  is  within 
+  1  sample  spacing  of  nominal.  In  fact,  the  algorithm  will  perform  nearly  as  well  if 
the  misregistration  is  within  +2  sample  spacings  of  nominal,  but  larger 
misregistrations  up  to  about  3  or  4  sample  spacings  will  be  interpreted  as  I  or  0. 
Larger  misregistration  will  produce  results  which  vary  about  the  nominal. 

However,  if  the  misregistration  can  be  expected  to  vary  smoothly  as  in  a  tracking, 
a  down-looking  or  a  range-tracking  application,  as  opposed  to  discontinuously  as  in 
a  forward-looking  application,  any  degree  of  misregistration  can  be 
accommodated  by  delaying  one  of  the  video  signals  relative  to  the  other. 

Assuming  that  the  appropriate  relative  delay  is  initially  set,  whenever  a 
misregistration  equal  to  I  sample  spacing  is  calculated,  a  delay  is  introduced  so 
that  the  calculated  result  is  zero.  Thus,  so  long  as  the  misregistration  varies 
somewhat  smoothly,  continued  high-accurate  correlation  is  possible.  Similarly,  if 
the  alignment  between  two  ranging  sensors  is  not  well  maintained,  delays  can  be 
incororated  to  account  for  expected  variations. 

To  reiterate,  the  limitations  of  PPC  described  so  far  are  that  the  two  signals 
must  be  known  a  priori  to  match  and  that  expected  misregistrations  are  either 
limited  or  smoothly  varying.  If  a  sufficiently  large  and  complex  system  were 
allowable,  even  the  latter  limitation  could  be  avoided.  Generally,  a  fixed 
rectangular  window  is  chosen,  but  the  only  real  requirement  is  for  the  window  to 
have  flat  horizontal  edges  and  that  the  horizontal  width  be  constant  although  the 
vertical  edges  may  curve  and  the  window  may  even  be  a  composite  of  several 
windows.  One  final  limitation  is  that  normalization  is  not  performed  so  the  two 
input  signals  should  be  well  matched  in  dc  level  and  gain. 

CALCULATION  RATE 

With  a  pipeline  implementation,  the  rates  at  which  correlation  results  may 
be  produced  is  limited  by  the  slowest  single  operation.  For  this  algorithm,  the 
slowest  operation  sets  a  limit  to  8  million  correlations  per  second  with 
commercial  chips  and  chip  sets.  This  corresponds  to  maximum  input  video 
bandwidths  of  4  MHz.  Alternatively,  CAI  has  developed  an  algorithm  to  perform 
normally  slow  arithmetic  functions  at  rates  up  to  100  million  operations  per 
second.  This  is  another  pipeline  implementation  and,  a  such,  is  most  useful  when  a 


steady  steam  of  high-rate  calculations  are  required.  With  such  special 
electronics,  PPC  can  be  performed  with  input  and  output  bandwidths  of  50  MHz. 


APPENDIX  C 

COMPUTER  SIMULATION  PROGRAM 


INTRODUCTION 

A  considerable  portion  of  the  effort  expended  under  this  program  went  to 
the  development  and  execution  of  a  computer  program  to  simulate  pipeline 
precision  correlation.  This  appendix  describes  the  program,  presents  listings  of 
the  program  and  results,  and  analyzes  the  results. 

It  is  important  to  note  that  simulation  programs  are  inherently  suspect 
because  it  is  impossible  to  model  all  possible  aspects  of  a  system  without  actually 
building  the  system.  The  programmer  must  make  some  sort  of  evaluation  of  what 
factors  are  important,  what  details  of  those  factors  are  important  and  how  to 
most  effectively  model  them,  considering  both  accuracy  and  program  size  and 
speed. 

This  program,  i r  addition  to  the  general  problem  of  simulation,  involves  two 
other  difficulties.  First,  processing  architecture  of  a  computer  is  not  well  suited 
to  simulate  a  pipeline  processing  architecture.  Second,  the  simulation  involves 
performing  calculations  on  scenes.  Just  defining  adequate  descripters  of  scenes  is 
the  subject  of  a  vast  amount  of  research.  Actually,  generating  a  reasonable 
facsimile  of  a  scene  in  general  is,  then,  a  premature  undertaking.  There  is, 
however,  no  alternative  but  to  make  the  attempt  and  decide  sub,  tively  if  the 
results  are  "scene  like".  The  results  of  any  process  dependent  on  such  artificial 
scenes  can  only  be  considered  a  general  indication,  since  there  is  no  way  to  prove 
that  the  scenes  could  be  real  scenes,  let  alone  real  scenes  likely  to  be 
encountered. 

The  sole  justification  of  results  from  this  simulation  program  is  that  they 
basically  agree  with  results  obtained  from  a  hardware  demonstration  unit  which 
implemented  a  more  standard  "one  target  window"  correlator,  rather  than  a 
pipeline  correlator. 


This  simulation  program  was  also  a  design  process  to  determine  the  values 
which  several  variable  parameters,  such  as  window  size  and  digitization  level 
(number  of  bits),  shold  take  in  the  preliminary  hardware  design.  The  factors 
against  which  precision  pipeline  correlation  were  tested  included  actual  image 
shift  and  signal-to-noise  ratio.  Other  factors  such  as  lens  MTF,  high-frequency 
vibration  and  other  image  degradations  were  not  tested.  These  factors  can  only 
be  adequately  modeled  with  transform  manipulations  and  were  not  considered  cost 
effective  for  the  purposes  of  this  study. 

The  program  was  initially  written  with  very  broad  calculations  in  the  sense 
that  many  combinations  of  design  parameters  were  calculated.  At  each  stage, 
choices  were  made  that  eliminated  several  combinations.  As  the  design  choices 
were  narrowed,  the  testing  against  image  shift  and  noise  was  broadened.  This 
approach  does  not  guarantee  that  final  choices  for  design  parameters  were 
optimal  for  all  cases,  but,  as  will  be  seen,  the  choices  were  either  so  clear-cut  for 
limited  cases  that  there  is  little  likelihood  of  a  bad  choice  or,  towards  the  end 
with  broader  tests,  so  close  as  to  make  little  difference. 

DESCRIPTION 

The  program  consists  of  two  sections.  The  first  section  generates  imagery 
and  the  second  operates  on  the  imagery  to  simulate  pipeline  correlation. 

CAI  has  in  its  library  a  program  to  generate  large  images  which  includes  the 
capability  to  superimpose  targets  on  background  with  horizons,  define  ground  and 
sky  clutter  statistics,  set  independent  motions  for  .target,  ground  and  sky,  apply 
MTF's  and  other  filters  and  so  forth.  Because  of  its  versatility,  it  is  a  large 
program  and  not  well  suited  to  the  requirements  of  this  study.  Pipeline 
correlation  by  its  nature  requires  continuous  streams  of  video  data.  Hence,  rather 
than  generate  scenery  on  a  frame-by-frame  basis,  the  approach  taken  was  to 
generate  scenery  "on  the  fly".  This  approach  requires  a  different  set  of  image 
descriptors  than  he  generation  of  a  full  frame.  The  primary  descriptor  is 
probability  of  edge,  that  is,  the  likelihood  that  a  pixel  intensity  will  be  unlike  any 
of  its  previously-calculated  neighbors.  Neglecting  initiation  and  edge  effects, 
each  pixel  intensity  being  generated  will  have  had  four  of  its  neighbors  previously 
generated.  The  generation  of  scenes  proceeds  in  the  following  manner.  First  ,  a 


uniformly-distributed  random  number  is  used  to  decide  if  the  current  pixel  will  be 
an  edge.  If  so,  the  average  of  its  four  previously  generated  neighbors  is  taken  and 
the  new  pixel  is  uniformly  allowed  to  be  any  value  outside  a  tolerance  of  this 
average,  but  between  0  and  I .  Other  distributions  such  as  an  agaussian  (an  upside- 
down  gaussian  curve)  led  to  intractable  mathematics.  If,  on  the  other  hand,  the 
current  pixel  is  not  to  be  an  edge  and  therefore  like  one  or  more  of  its  neighbors, 
the  next  step  is  to  decide  what  combination  of  neighbors  it  is  to  be  like.  Again,  a 
uniform  random  variable  detrmines  how  many  neighbors  it  will  be  similar  to:  one, 
two,  three  or  four.  If  four,  no  further  decisions  are  required.  If  one  or  three,  yet 
another  uniform  random  number  determines  which  of  the  four  is  used  or  excluded, 
respectively.  Finally,  if  two,  a  uniform  random  number  determines  which  pair. 
Once  the  combination  of  pixels  to  be  used  has  been  determined,  the  average  of 
their  intensities  is  calculated.  Using  this  average  as  the  mean  and  a  linear 
function  of  the  number  of  pixels  as  the  standard  deviation,  a  random  number  of  a 
gaussian  distribution  is  calculated.  The  random  number  is  used  as  the  intensity  of 
the  current  pixel.  This  procedure  is  reasonably  efficient  in  calculation.  Its 
purpose  is  to  allow  small  gray  shade  differences,  as  well  as  edges,  to  propagate 
randomly  in  any  direction. 

As  one  scene  is  being  generated,  a  shifted  version  of  it  must  also  be 
generated.  The  second  version  is  produced  by  calculating  a  weighted  (depending 
on  shift)  average  of  two  neighboring  pixels  in  the  fist  version.  High-frequency 
effects  are  thought  to  be  adequately  considered  by  this  method.  It  may  appear 
that  the  averaging  method  avoids  high-frequency  effects.  If  it  is  examined  in 
reverse,  however,  as  though  the  first  were  derived  from  the  second  versions,  it 
would  seem  that  high-frequency  effects  are  continuously  generated;  the  version 
"one"  pixels  rarely,  if  ever,  being  an  average  of  version  "two"  pixels. 

Typical  images  generated  by  this  procedure  as  shown  in  fight  C.l  and  C.2. 
The  latter  shows  several  copies  of  the  original  at  varying  shifts. 

While  the  images  are  being  generated,  their  statistical  standard  deviations 
are  calculated.  These  values  are  used  to  calculate  a  noise  term  added  to  each 
pixel  and  yielding  the  desired  signal-to-noise  ratio. 
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FIGURE  C.l  SAMPLE  IMAGE 
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FIGURE  C.2  SAMPLE  IMAGE  WITH  VARIOUS  SHIFTS 


The  second  section  of  the  program  is  the  actual  simulation  of  pipeline 
correlation.  This  is  a  fairly  straightforward  mathematical  procedure.  To  simulate 
a  pipeline  architecture,  numerous  partial  results  must  be  stored  in  files.  To 
compare  design  parameters,  the  most  efficient  calculation  method  is  to  execute 
the  program  a  single  time  and  only  calculate  alternative  partial  or  final  results  at 
the  appropriate  point  in  the  operation.  For  example,  all  results  use  the  same 
scene  data,  so  the  scene  is  generated  once.  The  first  differentiation  occurs 
immediately  after  the  scene  generation  stage  where  loops  are  set  up  to  test 
different  terms  in  the  precision  correlation  algorihm.  Within  this  loop  is  a  loop  of 
different  digitization  levels  of  the  video  to  be  tested.  There  is  then  a  loop  to  test 
different  window  sizes.  Finally,  after  the  results  for  all  these  combinations  have 
been  calculated,  is  a  loop  to  test  for  different  lengths  of  a  running  average  to 
smooth  the  results  somewhat  and  improve  the  accuracy  of  the  final  output.  As 
can  be  seen  from  the  many  combinations  possible,  files  storing  partial  results  can 
be  quite  large.  All  of  these  loops  are  executed  for  each  new  line  of  video 
generated.  Over  many  lines,  the  statistics  of  final  results  are  calculated  and 
these  statistics  are  printed  for  each  combination  at  the  end  of  program 
execution.  During  the  study,  as  decisions  were  made,  possible  combinations  were 
eliminated.  Both  the  program  and  files  became  more  tractible,  allowing  more 
extensive  analysis  of  the  effects  of  noise  and  image  shift.  At  each  decision  point, 
the  program  was  modified  to  more  efficiently  calculate  the  remaining  options. 

m 

As  with  any  development,  during  modification  to  optimize  the  program,  new 
and  more  efficient  calculation  procedures  were  conceived  and  implemented  so  the 
modifications  were  not  always  limited  to  deleting  existing  code.  Also,  errors  were 
found  which  for  the  most  part  had  minimal  effect.'  Those  results  which  were,  or 
may  have  been  affected  by  the  errors  not  detected  until  later  ore  marked  as 
suspect.  In  general,  it  was  not  worthwhile  backtracking  and  rerunning  such 
executions  as  the  errors  did  not  affect  the  conclusions. 

Listings  of  the  major  versions  of  this  program  will  be  found  at  the  end  of  this 
appendix.  The  first  listing  includes  a  few  utitlity  subroutines  such  as  a  gaussian 
distribution  random  number  generator  and  an  image  printing  routine. 


RESULTS 


Figures  C.3  through  C.5  are  graphs  of  some  of  the  later  results  whe  most  of 
the  design  parameters  had  been  decided.  Figure  C.3  shows  RMS  error  in  fractions 
of  a  pixel,  versus  window  size.  Not  surprisingly,  larger  window  sizes  provided 
better  performance,  although  the  ratio  of  improvement  is  not  as  predictable  as 
may  have  been  expected.  The  crossovers  and  general  nonlineraity  of  improvement 
is  probably  due  to  the  high-frequency  effects  mentioned  previously,  where  large 
windows  may  include  some  scene  information  which  causes  loss  of  accuracy  not 
seen  by  a  smaller  window.  Figure  C.4  shows  the  effects  of  digitization  level  on 
rms  error.  Clearly,  at  4  bits,  there  is  a  severe  breakpoint.  Binary  l-bit 
correlation  is  even  worse,  as  will  be  illustrated  below.  Because  4  bits  is  a 
breakpoint,  and  thus  some  what  dangerous  to  use,  the  design  effort  used  8  bits  of 
digitization.  Some  fairly  strong  performance  is  shown  in  figure  C.5,  where  error 
versus  signal-to-noise  ratio  is  graphed.  Note  that  even  at  signal-to-noise  ratio  of 
2.5:1  the  performance  is  still  adequate  for  stabilization. 

There  were  essentially  five  stages  of  analaysis  in  the  simulation  study. 
However,  before  presenting  and  analyzing  the  results,  some  clarifications  should 
be  made.  Initially,  three  different  sets  of  terms  were  used  in  the  precision 
correlation  algorithm  to  attempt  to  determine  the  optimum  set.  However,  the 
sets  chosen  exhibited  a  preference  for  shift  in  one  direction.  Therefore,  results 
for  positive  and  negative  shifts  were  calculated  for  each  set  of  terms.  The  two 
cases  are  referred  to  as  "Type  I"  and  "Type  2"  on  the  results  printouts.  Because  a 
relatively  wide  diffrence  existed,  during  the  first  two  stages  of  program  results, 
the  third  and  fourth  stages  calculated  results  both  ways,  and  determined  which 
shift  direction  was  the  more  likely,  based  on  partial  results.  In  the  fifth  stage,  a 
set  of  terms  that  had  been  worked  our  to  be  independent  of  shift  direction  were 
used  and  shown  to  exhibit  quite  good  performance. 

The  results  printouts  present  four  statistics  for  each  combination  of 
parameters.  The  statistics  are,  in  order:  the  average  error,  the  rms  error,  the 
maximum  error  during  the  trials,  and  the  standard  deviation  of  error.  For  the 
first  two  stages,  however,  the  average  is  not  divided  by  the  number  of  trials  and 
thus  is  the  sum  of  errors,  and  the  second  statistic  is  neither  divided  by  the  number 
of  trials  nor  is  the  square  root  taken  and,  hence,  this  statistic  is  the  sum  of  the 


FIGURE  C.4  RMS  ERROR  VS  SHIFT  AND  DIGITIZATION  LEVEL 


squares  of  errors.  The  purpose  of  leaving  the  statistics  as  sums  is  to  allow  their 
combined  statistics  to  be  more  easily  calculated. 

Obviously  bad  results,  those  indicating  shifts  greater  than  5  pixels,  were 
artificially  set  to  five.  The  integer  at  the  end  of  each  line  is  the  number  of  trials 
considered  valid.  A  minor  program  error  actually  added  7  or  8  to  this  number,  as 
indicated. 

At  the  first  stage  of  analysis  (figure  C.6),  the  second  and  third  sets  of 
algorithm  terms  were  eliminated  as  showing  poorer  performance  relative  to  the 
first  set.  In  terms  of  digitization,  the  low  counts  of  valid  trials  for  I -bit 
digitization,  leads  to  the  conclusion  the  binary  correlation  is  a  poor  choice.  It  is 
also  noted  that  windows  larger  than  8x8  provide  only  minimal  improvements  in 
performance.  Finally,  it  is  clear  from  these  first  stage  results  that  smoothing  by 
taking  a  running  average  of  results  offers  reasonable  improvements  in 
performance.  All  of  this  analysis  is  based  upon  images  with  0.1  pixels  of  relative 
shift  between  them  and  no  noise  imposed  on  the  video  signal. 

The  second  stage  maintains  the  zero  noise  condition.  Only  the  first  set  of 
algorithm  terms  is  tested  with  a  more  limited  set  of  window  sizes  and  digitization 
levels,  and  the  smoothing  running  average  length  is  fixed  at  eight.  Conversely, 
the  number  of  shift  positions  is  expanded  to  include  relative  shifts  of  0.1,  0.2,  0.3, 
0.7  and  0.9.  Note  that  the  first  and  last  statistics,  sum  of  errors  and  standard 
deviation,  for  the  6x6  windows  in  the  first  three  blocks  are  suspect  because  of  an 
error  detected  in  a  later  stage.  The  conclusion  fom  this  stage  of  analysis  is  that 
the  preference  of  the  set  of  terms  for  one  shift  direction  is  intolerable  and  other 
terms  will  have  to  be  defined  which  are  shift-direction  independent.  Further,  2- 
bit  digitization  is  unacceptable  with  windows  less  than  8x8,  and  4-bit  digitization 
is  only  marginally  acceptable. 

At  the  third  stage,  the  technique  of  calculating  both  cases  of  shift  direction 
and  choosing  the  most  likely  is  tested  with  the  same  parameters  as  for  the  second 
stage.  Nate  that  the  first  two  statistics  are  now  proper  mean  error  and  rms 
error.  Results  are  considerably  improved.  Although  even  a  6  x  6  window  with  2- 
bit  digitization  produces  acceptable  results  for  stabilization,  noise  has  still  not 
been  considered.  Therefore,  larger  windows  and  greater  resolution  is  still 
preferred. 
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0.390  -0.069 

0.015  ! 

!  -5.783 

0.331  -0.061 

0.014  ! 

!  108 

10X10 

-0.025 

0.001 

0.009 

0.002  I 

-5.077 

0.254  -0.034 

0.013  ! 

-4.290 

0.392  -0.067 

0.015  ! 

-5.801 

0.333  -0.061 

0.014  ! 

108 

PIPELINE  CORRELATION  RESULTS! 
DIGITIZING  TO  10  SITS 
ALOORITHH  2  TYPE  1 


RUNNING  AVE  LENGTH  1248 
HINOOH 


4X  4 

-1.013 

7.174 

1.232 

0.238 

-4.000 

4,444 

0.470 

0.104 

-7.1B1 

2.583 

0.382 

0.140 

-7,517 

1.144  -0.208 

0.077 

107 

7*  7 

4.202 

17.540 

1.470 

0.401 

-0.474 

0.185 

1.130 

0.202 

-3.218 

5.147 

0.743 

0.214 

-1.047 

2.342  0.448 

0.147 

106 

BX  B 

-0.424 

11.318 

-2.087 

0.324 

-5.444 

4.202 

-1.113 

0.234 

-8.210 

3.441 

-0.424 

0.142 

-4.701 

1.430  -0.340 

0.107 

108 

OX  0 

3.837 

7.078 

1.514 

0.234 

-1.223 

3.542 

0.741 

0.181 

-1.038 

1.724 

0.344 

0.123 

-0.043 

0.733  -0.140 

0.082 

108 

10X10 

-0,384 

7.040 

-1.528 

0.234 

-5.444 

3.001 

-0.010 

0.184 

-3.837 

2.312 

-0.530 

0.134 

-3.503 

1.347  -0.331 

0.100 

106 

PIPELINE  CORRELATION  RESULTS! 

BI6ITIZING  TO  IQ  BITS 
ALOORITHH  2  TYPE  2 

RUNNING  AVE  LENGTH  1  .2  4  B 


BINBOH 
OX  0  ! 

0.095 

0.004 

0.05B 

0.008  ! 

-4.058 

0.244 

-0.059 

0.013  1 

-4.040 

0.365  -0.068 

0.015  1 

-3.481 

0.208  -0.042 

0.014  1 

I  108 

7X  7  ! 

0.051 

0.010 

0.044 

0.009  ! 

-5.004 

0.233 

-0*068 

0.014  1 

!  -4.018 

0.361  -0.071 

0.015  ! 

!  -5.30* 

0.280  -0.040 

0.013  1 

!  108 

BX  B  1 

-0.047 

0.011 

-0.043 

0.010  1 

-3.100 

0.244 

-0.074 

0.015  ! 

-4.123 

0.376  -0.074 

0.016  « 

-5.434 

0.310  -0.048 

0.013  1 

108 

OX  0  ! 

-0.050 

0.012 

0.074 

0.010  ! 

I  -5.100 

0.243 

-0.077 

0.014  ! 

!  -4.270 

0.393  -0.077 

0.016  1 

!  -5.783 

0.333  -0.044 

0.015  ! 

1  108 

10X10  ! 

-0.057 

0.133 

0.302 

0.038  1 

-3.108 

0.317 

-0.155 

0.024  ! 

-4.304 

0.420  -0.123 

0.022  1 

!  -3. BIO 

0.347  -0.080 

0.018  1 

I  107 

PIPELINE  CORRELATION  RESULTS! 
BIBITIZING  TO  10  BITS 
ALOORITHH  3  TYPE  1 


i 


RUNNING  AVE  LENGTH  1  2  4  S 

HINDOU 


6X  6  1 

-5.041 

13.872 

-1.330 

0.355 

!  -10.043 

8.828 

-1.013 

0.270 

>-10.083 

5.108  -0.831 

0.104 

!-12.SP4 

2.012 

-8.495 

0.116  1 

!  107 

-7X  7  ! 

1  -2.179 

4.328 

1.345 

0.100 

!  -7.154 

2.474 

0.477 

0.134 

!  -0.232 

1.788  -0.315 

0.004 

1  -8.002 

1.144 

-8.217 

0.071  ! 

1  102 

8X  8  f 

3.266 

6.503 

1.542 

0.244 

!  -1.782 

2.042 

0.820 

0.143 

1  -3.154 

1.377  0.302 

0.100 

!  -2.445 

0.403 

-8.161 

0.864  ( 

107 

9X  9  > 

0.878 

11.177 

1.800 

0.322 

1  -4.107 

5.307 

0.817 

0.220 

!  -4,471 

2.723  0.717 

0.133 

!  -3.072 

1.373 

-8.330 

0.107  ! 

!  1*7 

10X10  1 

-2.273 

6.450 

-1.318 

0.243 

1  -7.312 

3.440 

-0.730 

0.171 

I  -7.813 

2.044  -0.424 

0.118 

1  -7.382 

1.242 

-8.279 

0.883  1 

104 

PIPELINE  CORRELATION  RESULTS! 
BIBITIZING  TO  10  SITS 
ALOORITHH  31YPE  2 


RUNNING  AVE  LENGTH  1 

2 

4 

8 

HINSON 
4X  4  • 

-0.300 

0.044  -0.104 

0.024 

-3.441 

0.320  -0.131 

0.021 

-4.301 

0.420  -0.110 

0.010 

-3.043 

0.3S4  -0.082 

0.014 

188 

7X  7  ! 

-0.477 

0.414  -0.428 

0.042 

-5.734 

0.503  -0.340 

0.043 

-4.711 

0.527  -0.224 

0.032 

-4.223 

0.423  -0.130 

0.024 

188 

OX  8  > 

-0.033 

0.103  0.288 

0.031 

-3.088 

0.308  -0.101 

0.025 

-4.102 

0.404  -0.088 

0.021 

-5.720 

0.334  -0.073 

0.018 

197 

OX  0  > 

0.084 

0.033  0.150 

0.017 

-4.044 

0.242  -0.070 

0.018 

-4.140 

0.384  -0.070 

0.018 

-5.400 

0.323  -0.044 

0.013 

188 

10X10  > 

-0.211 

0.030  -0.112 

0.017 

-3.243 

0.280  -0.108 

0.017 

-4.470 

0.422  -0.002 

0.018 

-3.073 

0.337  -0.074 

0.014 

187 

! 


PIPELINE  CORRELATION  RESULTS! 
DIGITIZING  TO  8  PITS 
ALSOR1THN  1  TYPE  1 


RUNNING  AVE  LEN6TH  1248 
HINDOO 


6X  6  ! 

-1,009 

3.413  -0.905 

0.183  ! 

!  -4.004 

2.128  -0.538 

0.129  ! 

'  -6.482 

1.298  -0.397 

0.092  ! 

!  -6.731 

0.952  -0.241 

0.070  ! 

1  106 

7  X  7  ! 

0.952 

1.845 

0.740 

0.130  ! 

!  -4.009 

0.924  0.341 

0.085  ! 

!  -6.225 

0.742  -0.205 

0.060 

1  -5.454 

0.485  -0.149 

0.042  1 

1  107 

IX  0  < 

3.402 

3.339 

1.045 

0.173  ! 

-1.422 

1.484  0.519 

0.116  ! 

-4.086 

0.734  0.220 

0*073 

-2.247 

0.261  0.095 

0.044 

108 

9X  9  ! 

-1.397 

3.0S8  -1.547 

0.148  ! 

I  -4.441 

1.977  -0.897 

0.121  ! 

!  -7.468 

1.394  -0.518 

0.090 

i  -6.233 

0.810  -0.279 

0.065 

1  108 

10X10  ! 

1.019 

0.347 

0.337 

0.058  1 

1  -4.033 

0.340  0.155 

0.044  • 

-4.405 

0.269  -0.097 

0.029 

-3.711 

0.160  -0.067 

0.017 

108 

PIPELINE  CORRELATION  RESULTS! 
DIGITIZING  TO  8  PITS 
AL80RITHH  1  TYPE  2 


RUNNING  AVE  LENGTH  12  4  8 

H1NDOH 


4X  6  1 

1.884 

3.415 

1.832 

0.177  ! 

!  -3.143 

1.825  0.900 

0.127  ! 

-4.121 

1.031  0.434 

0.090  ! 

!  -3.428 

0.526  0.195 

0.062  i 

108 

7X  7  1 

!  0.032 

0.020 

0.099 

0.014  ! 

1  -5.020 

0.263  -0.082 

0.014  ! 

-5.873 

0.3S1  -0.088 

0.017  * 

!  -4.620 

0,217  -0.062 

0.013  < 

108 

ax  8  < 

!  0.918 

0.474 

0.717 

0.079  ! 

-4.138 

0.596  0.455 

0.064  j 

-4.822 

0.452  0.193 

0.047  ! 

!  -4.300 

0.287  -0.084 

0.033  ! 

108 

9X  9  1 

!  -0.043 

0.014 

0.085 

0.011  ! 

-5.090 

0.263  -0.068 

0.015  ! 

-5.976 

0.357  -0.069 

0.015  i 

!  -5.606 

0.313  -0.062 

0.014  ! 

108 

10X10  ! 

-0.144 

0.004 

0.022 

0.006  ! 

-5.190 

0.269  -0.066 

0.014  1 

-6.289 

0.393  -0.071 

0.016  1 

-5.833 

0.338  -0.063 

0.015  ! 

108 

PIPELINE  CORRELATION  RESULTS! 
DIGITIZING  TO  8  PITS 
AL60RITHN  2  TYPE  1 


RUNNING  AVE  LENGTH  12  4  8 

NINDOW 


4X  4  ! 

-0.581 

0.382 

1.573 

0.279 

-5.570 

5.152 

0.870 

0.212  ! 

-4.477 

2.853 

0.491 

0.151  1 

-4.479 

1.135  -0.211 

0.063  i 

107 

7  X  7  1 

4.407  18.212 

1.789 

0.409 

!  -0.547 

9.728 

1.358 

0.300  ! 

-2.524 

5.427 

0.724 

0.223  ! 

!  -0.457 

2.459  0.400 

0.151  < 

!  104 

BX  8  > 

-0.227 

11.187 

-2.051 

0.322 

-5.270 

4.121 

-1.098 

0.233  ! 

-7.515 

3.342 

-0.607 

0.142  ! 

-5.952 

1.573  -0.339 

0.107  ! 

108 

9X  9  ! 

3.729 

7.558 

1.552 

0.242 

\  -1.324 

3.748 

0.765 

0.186  ! 

i  -1.498 

1.831 

0.382 

0.129  ! 

!  -0.728 

0.794  0.178 

0.084  ! 

1  108 

10X10  ! 

1.440 

4.593 

-1.173 

0.204 

-3.415 

2.449 

-0.770 

0.153  ! 

-3.777 

1.585 

-0.442 

0.114  ! 

-3.580 

0.979  -0.315 

0.089  ! 

105 

PIPELINE  CORRELATION  RESULTS! 
DIGITIZING  TO  8  PIT8 
ALGORITHM  2  TYPE  2 


RUNNING  AVE  LENGYH  1  2  4  B 

HINDOO 


6X  6  1 

0.141 

0.158 

0.284 

0.038  j 

-4.911 

0.3(4  -0.128 

0.031  ! 

-5*827 

0.375  -0.098 

0.024  1 

■  -5.043 

0.269  -0.082 

0.018  1 

107 

7X  7  ! 

1  0.054 

0.815 

0.686 

0.087  i 

I  -5.008 

0.494  0.388 

0.065  ! 

!  -5.538 

0.545  0.167 

0.049  ! 

!  -4.562 

0.325  -0.122 

0.035  ! 

!  108 

•X  8  I 

0.249 

0.139 

0.178 

0*036  ! 

-4.804 

0.317  -0.143 

0.031  ! 

-5.292 

0.335  -0.102 

0.026  ! 

-4.921 

0.269  -0.074 

0.020  1 

108 

9X  9  I 

!  -0.129 

0.090 

-0.195 

0.029  ! 

!  -5.174 

0.287  -0.135 

0.019  1 

1  -6.132 

0.384  -0.101 

0.018  ! 

!  -5.736 

0.331  -0.075 

0.015  ! 

!  108 

10X10  f 

-1.304 

2.337 

-1.512 

0.147  ! 

1  -4.353 

1.414  -0.785 

0.098  1 

-7.499 

1.036  -0.427 

0.069  1 

-7.016 

0.718  -0.241 

0.049  I 

107 

PIPELINE  CORRELATION  RESULTS! 

DIGITIZING  TO  8  BITS  I 

ALGORITHM  1  TYPE  1 


RUNNING  AVE  LENGTH  12  4  8 

HIN80U 


4X  6 

-3.204 

13.498 

1.447 

0.355 

-8.214 

7.744  -0.981 

0.257 

-8.925 

4.552  -0.822 

0.188 

•9.883 

2.342  -0.474 

0.115 

106 

7X  7 

-8.408 

7.323 

1.433 

0.240 

-5.481 

3.734  0.739 

0.179 

-6.947 

2.101  0.376 

0.124 

-3.551 

1.153  -0.217 

0.090 

103 

•X  • 

4.505 

8.334 

1.483 

0.273 

-0.542 

3.471  0.810 

0.184 

-1.404 

1.703  0.390 

0.125 

-0.744 

0.587  -0.143 

0.073 

108 

9X  9 

0.542 

11.372 

1.528 

0.324 

-4.504 

5.349  -0.740 

0.223 

-4.625 

2.570  0.667 

0.148 

-4.074 

1.197  -0.312 

0.098 

107 

10X10 

-3.429 

9.130 

-1.594 

0.289 

-8.448 

5.338  -0.941 

0.208 

-8.811 

2.894  -0.525 

0.142 

-8.239 

1.454  -0.303 

0.097 

107 

PIPELINE  CORRELATION  RESULTS! 

DIGITIZING  TO  8  PITS 
ALGORITHM  J  TYPE  2 

RUNNING  AVE  LENGTH  1  2  4  G 


M1NMN 
4X  4  ! 

-0.245 

1.340  -0.974 

0.111  ! 

!  -5.295 

0.847  -0.548 

0.073  < 

-4.034 

0.438  -0.303 

0.053  1 

!  -3.284 

0.408  -0.172 

0.037 

I  108 

7X  7  f 

-2.093 

3.320  -1.587 

0.174  ! 

!  -7.140 

2.029  -0.848 

0.120  ! 

•7.358 

1.249  -0.501 

0.084 

1  -4.793 

0.115  -0.274 

0.040 

1  108 

•X  •  1 

0.884 

0.410  0.415 

0.075  1 

-4.144 

0.478  0.279 

0.034  I 

-4.901 

0.380  0.101 

0.038 

-4.445 

0.233  -0.085 

0.023 

108 

9X  9  i 

8.425 

0.237  0.313 

0.047  ! 

!  -4.425 

0.333  0.134 

0.034  ! 

-5.718 

0.388  -0.098 

0.028 

1  -5.350 

0.314  -0.078 

0.021 

1  108 

10X10  ' 

-0.273 

0.298  -0.447 

0.052  ! 

i  -5.328 

0.423  -0.275 

0.039  i 

-4.413 

0.473  -0.172 

0.029  ! 

-3.943 

0.382  -0.114 

0.023 

!  108 

PIPELINE  CORRELATION  RESULTS 
DIGITIZING  TO  A  SITS 
ALGORITHM  1  TTPE  1 


RUNNING  AVE  LENGTH  I 
HINSON 

AX  A  I  2.050  P.204  1.SP7  0.2P1  I  -2.P44 

7*  7  I  I.PA4  2.673  9. SOS  0.15 A  !  -2.979 

GX  8  !  2.A12  3.753  -1.177  0.1B5  I  -2.40P 

PX  P  ■  -1.5A1  4.AA5  -l.PAA  0.207  <  -A.A18 

10X10  >  1.408  0.513  0.358  0.0A8  I  -3.A4A 


PIPELINE  CORRELATION  RESULTS! 
8I61TIZIH6  TO  A  SITS 
ALGORITHM  1  TYPE  2 

RUNNING  AVE  LENGTH  1 


WINDOW 
4X  4  t 

0.858 

1.185 

-0.483 

0.104 

-4.148 

7X  7  l 

0.400 

0.187 

0.243 

0.041 

!  -4.447 

•X  9  • 

1.445 

0.574 

0.442 

0.072 

-3.447 

tX  9  ! 

0.201 

0.044 

0.112 

0.025 

!  -4.858 

10X10  ! 

0.053 

0.081 

0.102 

0.027 

-4.999 

PIPELINE  CORRELATION  RESULTS! 
DIGITIZING  TO  A  SITS 
ALGORITHM  2  TYPE  1 


RUNNING  AVE  LENGTH  1 
HINDOH 


4X  4  t 

-1.505 

8.182  l.OPA 

0.275  I 

1  -4.437 

7X  7  t 

1.PP2 

17.243  -1.774 

0.3PP  ! 

1  -2.941 

BX  8  ! 

3.0P4 

7. 029  -1.191 

0.253  ! 

-1.938 

9X  9  ! 

4.904 

P.S37  1.324 

0.294  ! 

!  -0.13P 

10X10  ! 

2.758 

10.774  1.819 

0.313  1 

-2.282 

PIPELINE  CORRELATION  RESULTS! 
DIGITIZING  TO  A  SITS 
ALGORITHM  2  TYPE  2 


RUNNING  AVE  LENGTH  1 
HINSON 


4X  4  1 

0.414 

0.999 

0.500 

0.094  ! 

-4.424 

7X  7  1 

-0.239 

2.854 

-0.979 

0.143  1 

!  -5.344 

ex  e  f 

1.708 

1.070 

0.503 

0.098  ! 

-3.374 

n  9  ! 

0.470 

1.442 

-0.581 

0.115  ! 

1  -4.579 

10X10  1 

5.048 

4.934 

1.419 

0.209  ! 

-0.001 

PIPELINE  CORRELATION  RESULTS! 

DIGITIZING  TO  A  SITS 
ALGORITHM  3  TYPE  1 

RUNNING  AVE  LENGTH  1 
HINSON 

AX  A  !  -7. PAP  8.554  -1.384  0.272  1-12. P2A 

7X  7  I  -2.385  5.343  -1.0P7  0.221  I  -7.385 

SX  8  !  -0.38P  8.537  -1.PP7  0.281  I  -5.428 

PX  P  !  0.A10  7.2AP  1.331  0.25P  !  -4.423 

10X10  •  -1.017  5.240  1.A07  0.220  I  -A. 027 


2  4 

4.085  0.71P  0.1P3  !  -3.4A5  2.008  -0.427 
1.072  -0. 363  0.096  !  -5.273  0.A18  -0.222 
1.7AA  -0.AP4  0.12A  >  -4.843  O.PSO  -0.3A0 
2. SAP  -l.OPA  0.151  !  -A. 337  1.740  -0.A08 
0.43P  0.221  0.054  !  -3.241  0.252  -0.100 


2  4 

0.748  -0.38A  0.074  !  -5.115  0.503  -0.235 
0.33P  -0.135  0.03A  1  -5.280  0.387  -0.178 
0.481  0.383  0.058  <  -2.PP3  0.270  0.1A7 
0.282  -0.111  0.024  I  -4.7A4  0.257  -0.0P2 
0.302  -0.123  0.02A  !  -5.P13  0.37P  -0.103 


2  4 

4.841  -0.572  0.203  !  -7.7P3  2.A3A  -0.3A4 
8.417  -0.P27  0.278  I  -5.234  4.A2S  -0.A8A 
3.551  -0.AB2  0.180  I  -4.550  2.175  -0.47A 
4.1A2  0.73A  0.1PA  !  0.352  1.A24  0.3A1 
5.800  1.103  0.231  I  -2.0P5  2.P73  0.5P7 


2  4 

0.7 JO  0.2P4  0.071  !  -4.P50  0.507  0.132 

1.751  -0.54P  0.117  I  -4.0A1  0.P84  0.321 

0.A71  0.2P3  0.072  !  -1.AS3  0.3P5  0.17A 

0.758  -0.305  0.072  I  -4.2P4  0.407  -0.1P4 

3.4P7  1.305  0.180  I  -0.814  2.007  0.A37 


2  4 

A.04P  -0.827  0.204  (-13.577  4.432  -0.543 
3.234  0.734  0.15P  I  -7.330  1.857  0.3A1 
4. IPO  -1.08P  0.1P0  I  -3.313  1.PP3  -0.587 
3.A74  -0.683  0.180  I  -3.AP2  1.514  -0.353 
2.701  0.A53  0.148  !  -3.PP3  1.535  0.35P 


0.133  <  -3.3P3  1.134  -0.248  0.0P8  '  108 

0.058  !  -3.104  0.21P  -0.134  0.035  >  108 

0.082  !  -2.2P5  0.385  -0.211  0.05A  (  108 

0.113  !  -4. PSP  0.P37  -0.321  0.081  ■  108 

0.038  !  -2.44A  0.114  -O.OA7  0.023  !  108 


8 

0.04P  !  -3.117  0.188  -0.128  0.030  I  108 
0.035  !  -0.7P5  0.074  -0.085  0.025  !  108 
0.042  !  -2.337  0.133  -0.0A7  0.028  (  107 
0.021  !  -4.SAA  0.223  -0.0A7  0.017  I  108 
0.023  !  -3.523  0.318  -0.07P  0.018  >  108 


8 

0.13P  (  -7.453  1.11P  -0.224  0.073  !  107 

0.201  !  -1.442  1.517  0.304  0.118  >  10A 

1.134  !  -2.845  O.PAB  -0.207  0.0P1  !  107 

0.123  !  0.985  0.5PB  0.1A8  0.074  !  108 

0.1A5  I  -2.2BP  1 .743  -0.425  0.12A  I  104 


8 

0.031  I  -2. PSP  0.211  -0.134  0.035  I  107 

0.088  I  -1.8A7  0.41A  0.180  0.040  I  107 

0.05P  !  -1.714  0.231  0.094  0.043  !  108 

0.047  I  -4.301  0.27A  -0.120  0.031  •  108 

0.13A  I  -0.551  0.PA4  0.2P4  0.0P4  I  108 


8 

0.1SP  (-14.715  3.411  -0*348  0.114  I  103 
0.112  !  -5.AP5  1.113  -0.214  0.087  I  102 
0.124  !  -3.0AP  0.872  -0.275  0.077  I  107 
0.113  I  -4.114  0.P18  -0.320  0.084  !  107 
0.104  I  -A.2P5  0.955  -0.207  0.074  •  105 


PIPELINE  CORRELATION  RESULTS! 
DIGITIZING  TO  A  8 ITS 
ALGORITHM  3  TYPE  2 


RUNNING  AVE  LlNSTN  1 

2 

4 

8 

MINION 
4X  4  1 

1.994 

5.808  1.770 

0.231  ! 

-2.998 

2.943  0.842 

0.143  f 

-3.433 

1.503  0.394 

0.114  * 

-2.290 

0.718  -0.229 

0.079  1 

103 

n  7  » 

0. 144 

4.73P  -1 .274 

0.209  ( 

1  -4.888 

2.758  -0.879 

0.153  ! 

-4.003 

1.521  -0.501 

0.113  ' 

-3.488 

0.774  -0.279 

0.078  I 

>  108 

ex  e  i 

0.S2A 

3.484  -0.P78 

o.ieo  i 

-4.505 

1.844  -0.590 

0.125  t 

-4.781 

1.147  -0.349 

0.094  I 

-4.890 

0.744  -0.247 

0.071  1 

108 

9X  9  ! 

-0.881 

2.543  0.P44 

0.154  ! 

!  -3.101 

1.557  0.518 

0.110  ! 

!  -4.191 

1.073  -0.400 

0.082  ! 

-4.424 

0.425  -0.231 

0.045  < 

I  108 

10X10  f 

-2.000 

2.550  -1.189 

0.153  ! 

-7.033 

1.774  -0.448 

0.110  P 

-8.527 

1.401  -0.373 

0.082  P 

-9.199 

0.942  -0.244 

0.055  ! 

108 

PIPELINE  CORRELATION  RESULTS 
B1GITIZIN6  TO  4  BITS 
ALGORITHM  1  TYPE  1 


I 


\ 

fc 


I 


RUNNING  AVE  LENGTH  1 
WINDOW 


AX  A  ! 

-2.288 

10.256 

1.567 

0.307 

-7.352 

7X  7  ! 

-0.657 

4.217 

1.086 

0.197 

!  -5.685 

8X  8  j 

3.205 

3.241 

1.250 

0.171 

-1.828 

»X  »  ! 

1.380 

1.096 

0.604 

0.100 

!  -3.682 

10X10  > 

0.834 

1.257 

0.600 

0.108 

-4.269 

2 

4 

5.321  -0.864 

0,211  ! 

!  -7.360 

2.795  -0.617 

2.754  -0.944 

0.151  ' 

\  -6.851 

1.615  -0.577 

1.934  0.725 

0.133  ! 

'  -3.402 

1.076  0.371 

0.786  0.382 

0.078  i 

!  -3.524 

0.564  0.231 

0.954  -0.293 

0.085  ! 

-3.267 

0.593  0.229 

8 


0.146  I 

-7.151 

1.713 

-0.387 

0.107  j 

105 

0.105  1 

!  -3.366 

0.684 

-0.279 

0.073  ! 

!  106 

0.095  ! 

1.081 

0.456 

0.189 

0.064  ! 

l  108 

0.064  ! 

i  -0.501 

0.264 

0.166 

0.049  ! 

t  107 

0.068  » 

-0.888 

0.287 

0.162 

0.051  ! 

108 

I 

I 

I 


PIPELINE  CORRELATION  RESULTS: 
DIGITIZING  TO  4  PITS 
ALGORITHM  1  TYPE  2 


RUNNING  AVE  LENGTH  1  2 

WINDOW 


AX  A  1 

-1.562 

4.0A1  -1.322 

0.193  1 

!  -6.648 

2.612  -0.738 

0.143 

7X  7  1 

-0.833 

1.SA2  O.AOA 

0.120  ! 

!  -5.915 

1.168  -0.297 

0.088 

DX  8  1 

-1.045 

2. SPA  -0.PS7 

0.155  1 

-6.214 

1.679  -0.536 

0.111 

»X  9  1 

-0.884 

1.054  -O.ASA 

0.098  ! 

'  -6.094 

1.042  -0.389 

0.080 

10X10  1 

-0.723 

0.A2P  -0.283 

0.076  ! 

-5.882 

0.815  -0.309 

0.068 

-5.354 

4 

1.397  -0.394 

0.102  ( 

i  -7.303 

8 

0.977  -0.237 

0.067 

102 

-4.709 

0.659  -0.199 

0.065  1 

!  -4.673 

0.480  -0.183 

0.051 

108 

-2.696 

0.698  -0.268 

0.076  1 

-1.041 

0.295  -0.135 

0.051 

107 

-2.332 

0.510  -0.235 

0.065  1 

!  -1.389 

0.297  -0.136 

0.051 

10A 

-3.404 

0.442  -0.172 

0.056  • 

-2.350 

0.260  -0.123 

0.044 

108 

PIPELINE  CORRELATION  RESULTS! 
DIGITIZING  TO  4  DITS 
ALGORITHM  2  TYPE  1 


RUNNING  AVE  LENGTH  1  2  4  8 

WINDOW 


AX  A  f 

-7.P08 

14.440  -1.473 

0.338 

1-12.927 

8.740  -l.OPA 

0.258 

>-11.902 

5.017  -0.420 

0.185 

1-14.875 

4.479 

-0.470 

0.131  ! 

i  105 

7X  7  1 

-3.844 

20.294  -1.P37 

0.430 

1-10.870 

10.945  -1.22A 

0.302 

>-10.335 

5.844  -0.492 

0.211 

1  -8.110 

2.810 

-0.423 

0.143  ! 

1  104 

BX  8  1 

3.012 

13.083  -1.A33 

0.347 

!  -2.083 

7. 748  -0.90A 

0.247 

-0.970 

4.143  -0.309 

0.196 

!  1.330 

1.985 

0.278 

0.135  1 

!  107 

PX  P  ! 

S.2PA 

11.241  -1.747 

0.31P 

1  0.184 

3. 447  1.173 

0.225 

1.347 

2.473  0.592 

0.151 

I  3.451 

1.241 

0.349 

0.102  i 

!  105 

10X10  ! 

2.784 

A. 822  -1.174 

0.250 

!  -2.313 

3. 453  0.725 

0.183 

-0.388 

2.189  0.497 

0.142 

■  1.143 

1.243 

0.272 

0.107  1 

105 

I 


PIPELINE  CORRELATION  RESULTS! 
DIGITIZING  TO  4  BITS 
ALGORITHM  2  TYPE  2 


RUNNING  AVE  LENGTH  1248 
WINDOW 


AX  4  ! 

-3.633 

16.573  -1.645 

,0.390  ! 

-8.656 

8.891  -1.054 

0.275  ! 

-9.130 

5.244  -0.660 

0.203 

•  -15.064 

4.483  -0.373 

0.149 

i  103 

7X  7  ! 

!  -2.216 

11.454  -1.386 

0.325  ! 

1  -7.258 

5.684  -0.965 

0.219  1 

1  -7.076 

2.922  -0.470 

0.151 

-8.925 

2.023  -0.308 

0.109 

104 

BX  8  1 

!  -0.949 

9.681  -1.254 

0.299  1 

1  -6.127 

6.200  -0*838 

0.233  ! 

-2.091 

3.312  -0.551 

0.174 

-1.813 

1.858  -0.329 

0.130 

107 

9X  9  1 

!  -0.501 

13.017  1.800 

0.347  ! 

1  -5.661 

7.808  1.050 

0.244  ! 

I  -2.010 

4.464  0.814 

0.202 

-1.456 

2.334  0.457 

0.147 

107 

20X10  1 

i  -3.717 

12.054  -1.600 

0.332  I 

i  -0.881 

7.444  -1.167 

0.249  1 

-6.928 

3.780  -0.593 

0.176 

-5.440 

1.790  -0.334 

0.118 

108 

PIPELINE  CORRELATION  RE8ULT8! 
DIGITIZING  TO  4  DITS 
ALGORITHM  J  TYPE  1 


RUNNING  AVE  LENGTH  1  2 

WINDOW 


4X  4  1 

1  -5.170 

19.126  -1.850 

0.418 

f-10.131 

9.063  -1.114 

7X  7  1 

!  0.606 

13.904  -1.571 

0.339 

!  -4.371 

7.189  -0.951 

8X  8  1 

!  -3.865 

4.784  1.054 

0.207 

!  -0.960 

3.498  -0.455 

9X  9  i 

!  -6.277 

0.125  -1.433 

0.248 

1-11.423 

6.494  -1.344 

10X10  1 

!  -6.303 

5.864  -1.130 

0.224 

1 -11.439 

4.153  -0.850 

I 


0.274 

1-11.028 

4 

4.545  -0.677 

0.170 

J -17, 139 

8 

4.606  -0.599 

0.132 

1  102 

0.255 

f  -5.244 

3.825  -0.564 

0.182 

1  -4.902 

1.684  -0.374 

0.114 

!  103 

0.160 

I  -7.661 

2.120  -0.374 

0.121 

1  -6.237 

1.028  -0.250 

0.079 

!  104 

0.221 

1  -9.535 

4.380  -0.945 

0.181 

!  -8.220 

2.865  -0.630 

0.144 

!  104 

0.165 

!-10.065 

2.535  -0.507 

0.122 

!  -7.606 

1.418  -0.327 

0.090 

\  103 

PIPELINE  CORRELATION  RESULTS! 
DIGITIZING  TO  4  BITS 
ALGORITHM  3  TYPE  2 


PUNNING  AVE  LENGTH  I  2 

WINDOW 


6X  6 

2.757 

16.930 

1.800 

0.395 

F  -2.208 

8.461  0.941 

0.279 

7X  7 

-7.247 

18.094  -1.989 

0.404 

f-12.259 

11.126  -1.106 

0.300 

BX  8 

-3.872 

9.014 

1.456 

0.287 

•  -8.983 

5.029  -0.838 

0.199 

4X  9 

-5.669 

8.533  -1.331 

0.276 

■-10.818 

5.052  -0.789 

0.192 

10X10 

-7.294 

11.083 

1.567 

0.313 

>-12.488 

7.454  -1.094 

0.236 

-3.496 

4 

4.201  0.578 

0.195 

J  -7. 872 

8 

2.533  0.402 

0.133  ! 

104 

-13.603 

7.347  -0.706 

0.228 

1-14.610 

4.B49  -0.540 

0.143  ! 

104 

-7.392 

2.596  -0.446 

0.139 

!  7.510 

1.444  -0.300 

0.092  > 

105 

-9.990 

3.052  -0.477 

0.140 

!  -8.833 

1.803  -0.334 

0.100  \ 

105 

-11.597 

5.292  -0.098 

0.194 

!  -8.809 

3.093  -0.550 

0.140  I 

105 

L 


PIPELINE  CORRELATION  RESULTS! 
DIGITIZING  TO  2  (ITS 
ALGORITHM  1  TYPE  1 


RUNNING  AVE  LENGTH  1248 
HINDOU 


4X  4  1 

-3.787 

2.298 

0.547 

0.142  ! 

-8.924 

2.316  -0.436 

0.121 

-10.379 

1.996  -0.306 

0.096  ! 

-6.512 

0.848  -0.223 

0.065 

88 

7X  7  ! 

1  -0.312 

1.918 

-0.529 

0.133  ! 

!  -5.441 

1.510  -0.441 

0.107 

-5.449 

0.994  -0.268 

0.082  1 

i  -1.773 

0.434  0.142 

0.061 

i  104 

8X  8  ! 

1.432 

1.452 

0.400 

0.115  1 

!  -3.457 

1.033  -0.238 

0.092 

-4.913 

0.796  -0.205 

0.073  • 

-2.051 

0.456  -0.231 

0.062 

107 

9X  9  1 

!  0.027 

0.989 

0.400 

0.096  ! 

f  -5.052 

0.929  -0.233 

0.080 

-5.434 

0.693  -0.201 

0.062  ! 

!  -2.859 

0.365  -0.184 

0.052 

1  108 

10X10  ! 

-0.843 

0.762 

0.300 

0.084  ! 

-5.950 

0*867  -0.213 

0.071 

-6.070 

0.639  -0.184 

0.052  ! 

-3.651 

0.320  -0.153 

0.043 

108 

PIPELINE  CORRELATION  RESULTS! 

DIGITIZING  TO  2  (ITS 
ALGORITHM  1  TYPE  2 

RUNNING  AVE  LENGTH  1248 
HINDOU 

AX  A  ■  -4. 245  2. 143  0.471  0.13A  !  -9.383  2.343  -0.4A0  0.119  !  -5.324  1.272  -0.2A4  0.097  !  -3.A50  0.A79  -0.193  0.072  •  91 

7X  7  !  -1.857  1.378  -0.433  0.112  !  -A. 970  1.519  -0.483  0.100  !  -2.912  0.724  -0.217  0.077  !  -1.254  0.343  -0.121  0.055  ■  102 

SX  8  I  -0.389  0.931  -0.2A7  0.093  I  -5.484  1.042  -0.279  0.084  >  -1.038  0.522  -0.170  0.0A9  !  0.223  0.383  -0.211  O.OAO  •  107 

9X  9  I  -0.944  0.710  -0.225  0.081  !  -4. 055  0.913  -0.258  0.073  !  -2.885  0.450  -0.144  0.059  !  -1.352  0.294  -0.148  0.051  !  108 

10X10  !  -1.477  0.470  0.344  0.077  !  -4.790  0.923  -0.221  0.048  !  -3.890  0.429  -0.148  0.052  !  -1.804  0.244  -0.151  0.045  !  108 


PIPELINE  CORRELATION  RESULTS! 
DIGITIZING  TO  2  BITS 
ALGORITHM  2  TYPE  1 


RUNNING  AVE  LENGTH  1248 
HINDOU 


4X  4  ! 

!  -7.422 

2.495  -0.475 

0.142 

1-12.544 

3.091  -0.449 

0.123  ! 

1  -9.519 

1.904  -0.288 

0.099  ! 

!  -5.559 

0.749  -0.196 

0.065  f 

42 

7X  7  1 

!  -4.835 

3.175  -0.400 

0.159 

1-11.943 

3.555  -0.547 

0.144  1 

1  -8.028 

2.200  -0.471 

0.122  j 

!  -5.708 

1.116  -0.335 

0.087  1 

47 

8X  8  > 

-2.440 

2.851  0.547 

0.141 

!  -7.448 

2.501  -0.344 

0.135 

-4.212 

1.404  -0.252 

0.107  ! 

-1.720 

0.770  -0.221 

0.083  ! 

48 

9X  9  ! 

-2.487 

5.959  -1.100 

0.234 

!  -7.447 

3.487  -0.485 

0.171 

!  -5.098 

2.154  -0.370 

0*133  1 

!  -2.780 

1.269  -0.304 

0.105  1 

80 

10X10  1 

-0.990 

4.930  -0.850 

0.213 

!  -4.011 

3.730  0.542 

0.177 

-3.991 

2.342  0.377 

0.143  ! 

-0.871 

1.381  0.331 

0.113  ! 

B5 

PIPELINE  CORRELATION  RESULTS! 
DIGITIZING  TO  2  BITS 
ALGORITHM  2  TYPE  2 


RUNNING  AVE  LENGTH  1248 
HINDOU 


6X  6  ! 

-1.935 

2.413 

0.471 

0.155 

!  -7.119 

2*200 

0.384 

0.127  ! 

-2.580 

1.202  0.267 

0.103  1 

-1.085 

0.638  0.185 

0.076  1 

42 

7X  7  ! 

!  -2.470 

2.395 

-0.600 

0.147 

!  -7.862 

2.050 

-0.454 

0.117 

1  -2.445 

0.811  -0.250 

0.083  1 

!  -1.266 

0.447  0.160 

0.063  ! 

47 

BX  8  I 

-1.348 

2.070 

0.471 

0.138 

!  -6.564 

1.821 

-f .379 

0.115  1 

-0.937 

0.786  -0.229 

0.085  ! 

1.324 

0.648  -0.299 

0.076  ! 

48 

9X  9  ! 

!  -2.184 

3.194 

0.614 

0.171 

f  -7.386 

2.583  -0.525 

0.139  ! 

1  -2.843 

1.339  -09376 

0.108  ! 

1.264 

0.966  -0.301 

0.094  ! 

80 

10X10  ! 

-1.985 

3.193 

0.700 

0.171  ! 

I  -7.255 

2.954  -0.579 

0.151  ! 

-4.121 

1.678  -0.397 

0.119  t 

2.329 

0.956  -0.235 

0,092  ! 

85 

PIPELINE  CORRELATION  RESULTS! 
DIGITIZING  TO  2  BITS 
ALGORITHM  3  TYPE  1  -- 


RUNNING  AVE  LENGTH  1 

HINDOU 

4X  4  !  -4.049  3.417  -0.700 

0.174 

1-11.203 

2 

3.479  -0.450 

0.144  ! 

-7.878 

4 

2.009  0.280 

0.115  1 

-5.940 

8 

1.014  -0.195 

0.080  1 

S3 

7X  7  !  -5.105 

3.366  0.650 

0.170 

f -10.155 

3.041  0.392 

0.139  1 

1  -8.974 

2.298  -0.345 

0.120  ! 

!  -5.288 

1.223  -0.218 

0.094  1 

42 

BX  8  1  -4.012 

4.160  0.900 

0.188 

1-11.162 

3.743  0.445 

0.155  I 

-7.444 

2.253  0.350 

0.127  1 

-3,297 

1.464  -0.318 

0.112  ! 

47 

9X  9  !  -4.535 

5.059  -1.100 

0.208 

t-11.631 

4.581  -0*650 

0.174 

1  -9.324 

3.019  -0.412 

0*143  ! 

I  -4.584 

1.761  -0.347 

0.120  ( 

80 

10X10  !  -4.774 

5.867  0.733 

0.229 

f  -9.889 

4.898  -0.567 

0.192  1 

!  -7.141 

3.161  -0.427 

0.158  1 

'  -1.044 

1.941  -0.286 

0.134  1 

84 

PIPELINE  CORRELATION  RESULTS! 
DIGITIZING  TO  2  BITS 
ALGORITHM  3  TYPE  2 


RUNNING  AVE  LENGTH  1  2  4  8 

HIN80H 


6X  6  1 

-4.798 

3.567  -0.600 

0.170 

•—11 .948 

3.477  -0.550 

0.148  ! 

-8.085 

2.041  -0.313 

0.115  ! 

-6.312 

1.199  -0.237 

0.088  ! 

53 

7X  7  ! 

I  -4.111 

2.995  0.567 

0.142 

f  -9.077 

2.934  -0.443 

0.142  ! 

i  -5.048 

1.650  -0.303 

0.114  i 

i  -3.192 

1.033  0.224 

0.093  t 

57 

BX  8  ' 

-5.220 

3.92G  -0.600 

0.184 

•-10*370 

4.042  -0.400 

0.168  ! 

-4.841 

2.860  -0.434 

0.150  t 

1.003 

1.447  -0.274 

0.124  ! 

49 

9X  9  ! 

-4.090 

4.795  -0.900 

0.203 

1-11.185 

4.821  -0.733 

0.184  ! 

-9.041 

3.557  -0.564 

0.161  ! 

0.500 

2.233  -0.374 

0.144  1 

•1 

10X10  ! 

•5.411 

3.679  0.614 

0.178 

1-10.649 

3.784  -0.547 

0.159  1 

-3.544 

2.217  -0.394 

0.139  1 

4,696 

1.877  -0.245 

0.124  ! 

79 

PIPELINE  CORRELATION  RESULTS 
DIGITIZING  TO  1  SITS 
ALGORITHM  1  TYPE  1 


RUNNING  AWE  LENGTH  1 


UlMDQU 
AX  4  ! 

-5.454 

1.175 

0.300 

0.090 

1-10.804 

?X  7  ! 

-1.710 

0.902 

0.329 

'.090 

!  -4.822 

•X  8  ! 

0.789 

0.0?' 

0.329 

0.087 

!  -4.324 

9X  9  • 

-0.712 

0.52.. 

0.233 

0.049 

!  -5.807 

10X10  ! 

-1.700 

0.407 

0.184 

0.059 

!  -4.794 

PIPELINE  CORRELATION  RESULTS! 
DIGITIZING  TO  1  SITS 
ALGORITHM  1  TYPE  2 


RUNNING  AVE  LENGTH  1 
HINSON 


4X  4 

-5.991 

1.177 

-0.300 

0.088 

(-11.104 

7X  7 

!  -1.943 

0.743 

0.200 

0.081 

‘  -7.054 

fX  8 

0.245 

0.474 

0.18. 

0.079 

!  -4.855 

»X  9 

1  -1.024 

0.444 

0.173 

0.045 

1  -4.121 

10X10 

-2.010 

0.375 

0.1.7 

0.054 

!  -7.102 

PIPELINE  CORRELATION  RESULTS! 
DIGITIZING  TO  1  SITS 
ALGORITHM  2  TYPE  1 


RUNNING  AVE  LENGTH  1 
HINSON 


4X 

4  I 

i  -8.090 

1.113  0.233 

0.0.9 

•-13.140 

7X 

7  1 

1  -4.212 

1.31.  0.300 

0.094 

(-11.242 

BX 

8  ! 

-5.501 

1..U  -0.500 

0.123 

(-10.551 

9X 

9  ! 

!  -4.549 

1 ,.59  0.400 

0.117 

(  -9.599 

10X10  1 

i  -3.703 

2.1.4  -0.700 

0.137 

!  -8.753 

PIPELINE  CORRELATION  RESULTS! 
DZ8ITIZ1N8  TO  1  SITS 
ALGORITHM  2  TYPE  2 

RUNNING  AVE  LENGTH  1 


UINDOU 
.X  .  1 

-..899 

1.217 

0.300 

0.085 

(-11.949 

7X  7  ! 

-5.555 

1.28. 

0.275 

0.09. 

(-10.405 

SX  8  I 

-5.985 

1.501  -0.300 

0.104 

(-11.035 

9X  9  ! 

-..493 

1.703  -0.500 

0.110 

(-11.543 

10X10  ! 

-..173 

1.494  -0.433 

0.103 

(-11*223 

PIPELINE  CORRELATION  RESULTS! 
DIGITIZING  TO  1  SITS 
ALGORITHM  3  TYPE  1 


RUNNING  AVE  LENGTH  1 
HINSON 


4X 

4  1 

-•*057 

1*225  0.300 

0.07. 

1-13.107 

7X 

7  ! 

!  -4*218 

2*083  -0.500 

0.12. 

1-11.2.8 

•X 

0  1 

-4.374 

2.558  -0.400 

0.142 

1-11.42. 

9X 

9  ! 

-7*584 

3.072  -0.400 

0.153 

(-12.434 

10X10  1 

-7.473 

3*138  -0*400 

0.155 

1-12.723 

PIPELINE  CORRELATION  RESULTS! 
DIGITIZING  TO  1  SITS 
ALGORITHM  3  TYPE  2 


RUNNING  AVE  LENGTH  1 
HINSON 


4X 

4  ! 

•  -9.114 

1.137  -0.300 

0.058 

1-14.1.4 

71 

7  1 

!  -5»48> 

2.758  0,500 

0.151 

1-10.733 

•X 

8  1 

i  -4*443 

2,9.7  -0.433 

0.154 

1-11. .93 

9X 

9  1 

1  -4.775 

3.304  -0.529 

0.1.3 

1-11. *25 

10X10  1 

-7.318 

3.1*.  -0.529 

0.157 

1-12.439 

> 


2 

1.804  -0.305 

0.082 

1-12.210 

4 

1.954  -0.240 

1.199 

0.214 

0.084 

!  -7.257 

1.091 

0.179 

0.854 

0.201 

0.079 

!  -4.430 

0.734 

0.142 

0.777  -0.150 

0.044 

(  -5.484 

0.480  -0.148 

0.778  -0.150 

0.057 

1  -4.583 

0.490  -0.147 

2 

4 

1.924  -0.350 

O.GBS  1 

-9.940 

1.520  -0.238 

1.084  -0.150 

0.074 

1  -5.899 

0.814  0.154 

0.809  -0.150 

0.074  1 

-3*490 

0.598  -0.130 

0.749  -0.150 

0.043  ! 

-4.941 

0.570  -0.130 

0.792  -0.150 

0.055  1 

-5.910 

0.588  -0.130 

2 

4 

2.115  0.183 

0.049 

(-11.878 

1.714  0.154 

2.044  0.225 

0.091 

!  -9.999 

1.435  0.237 

2.378  -0.350 

0.112 

!  -9.288 

1.872  -0.230 

2.078  -0.317 

0.107 

!  -8.334 

1.419  -0.221 

2.271  -0.450 

0.120 

!  -7.490 

1.483  -0.280 

2 

4 

2.028  0.183 

0.081 

(-10.484 

1.594  0*194 

1.950  0.225 

0.092 

!  -9.343 

1.529  0.237 

2.082  -0.305 

0.094 

1  -9.804 

1.424  -0.245 

2.374  -0.517 

0.103 

(-10.384 

1.921  -0.404 

2.1ft  -0.440 

0.097 

(-10.144 

1.789  -0.370 

2 

ft 

4 

2.171  0.250 

0.073 

1-11.845 

1*709  0.142 

2*440  -0*350 

0.114 

(-10.005 

2.041  0*343 

3*113  -0*400 

0.133 

1-10.213 

2.437  -0*434 

3*444  -0*414 

0.142 

(-11*434 

3*028  -0*532 

3.773  -0.588 

0.145 

(-11*573 

3*138  -0.550 

2 

4 

2.1*9  -0.250 

0.055 

(-12*902 

1.782  -0.1SS 

3.314  0.450 

0.144 

1  -9.471 

2.704  0.437 

3.477  -0.440 

0.143 

(-10*414 

2.909  -0.395 

3.891  -0.532 

0.155 

(-10.813 

3.274  -0.494 

3.742  -0.531 

0.147 

(  -9.423 

2.875  -0.479 

S 


0.073 

1-10.084 

1.272  -0.178 

0.055  1 

75 

0.075 

!  -5.215 

0.411  0.134 

0.058 

98 

0.072 

!  -2,477 

0.477  0.153 

0.042 

108 

0.059 

!  -3.499 

0.419  -0.139 

0.052 

100 

0.052 

1  -4.338 

0.414  -0.125 

0.045  1 

108 

s 


0.075 

-8.090 

0.933  -0.147 

0*057  ! 

77 

0.047 

1  -4.018 

0.444  0.125 

0.054  ! 

!  97 

0*044 

-1.873 

0.388  -0.142 

0.057  ! 

108 

0.054 

!  -3.043 

0.352  -0.128 

0.050  1 

!  108 

0*049  ( 

-3.991 

0.352  -0.119 

0.043  ! 

108 

8 


0.042 

l  -9.824 

1.143  -0.117 

0.048  ! 

18 

0.081 

1  -7.948 

1.044  0.123 

0.045  ! 

31 

0.100 

'  -7.383 

1.259  -0.147 

0.084  ■ 

47 

0.093 

!  -4.443 

1.114  0.174 

0.082  ■ 

58 

0.104 

!  -5.338 

1.114  -0.192 

0.088  ( 

72 

S 


0.071  ! 

-8.433 

1.008  -0.117 

0.034  ! 

IS 

0.082  ! 

!  -7.291 

0.951  0.143 

0.043  1 

31 

0.083  ! 

-4.909 

0.947  -0.219 

0.048  ! 

47 

0.092  ! 

I  -4.749 

1.144  -0.284 

0.082  1 

58 

0.088  1 

-5.842 

0.947  -0.244 

0.074  1 

72 

■ 


0.042  1 

-9.793 

1.117  -0.117 

0.044  1 

10 

0.102  1 

-7.954 

1.279  -0.213 

0.080  1 

24 

0.117  1 

-4.343 

1.270  -0.280 

0.091  1 

35 

0.130  1 

-7.427 

1.7*7  -0.373 

0.109  1 

47 

0.133  1 

-4.871 

1.829  -0.371 

0.114  1 

42 

s 


0.047 

(**10*850 

1.233  -0.142 

0.034  1 

10 

0.132 

I  -7*419 

1.803  0.333 

0.109  1 

30 

0.131 

(  -3.388 

1.414  0.3*8 

0.110  1 

37 

0.142 

(  -2*857 

1.703  -0.34* 

0.123  1 

49 

0.137 

1  -2*892 

1.732  -0.348 

0.124  ! 

44 

II. 37. OR  >1«S 

520. SS  ARU'Si  .83  CONNECT  HRS 
LOSSES  OFF  AT  It. 37. I I  ON  03MARS2. 


STAGE  II 


(  DIGITIZING  TO*  6  BITS 
ALGORITHM  1  TYPE  1 

I  (  RUNNING  AVE  LENGTH  8 


WINDOW 
AX  A  *- 

1.043 

0.391 

O._0?8  ! 

105 

7X  7  ! 

-57358 

1.42? 

-0.305 

0.103  ! 

103 

8X  8  1 

1.08A 

1.587 

0.269 

0.121  > 

104 

PIPELINE  CORRELATION  RESULTS! 
DIGITIZING  TO  6  BITS 
ALGORITHM  1  TYPE  2 


c 

( 

( 

( 

c 

c 


RUNNING  AVE  LENGTH  8 
WINDOW-^ — ^ 


AX 

A  'L_=0O58' 

0.641 

0.218 

103 

7X 

7  1  -0.A92 

0.394 

-0.174  0.060  ! 

105 

8X 

8  !  1.316 

0.608 

0.172  0.074  ! 

107 

PIPELINE  CORRELATION  RESULTS! 
DIGITIZING  TO  4  BITS 
ALGORITHM  1  TYPE  1 

RUNNING  AVE  LENGTH  8 


WINDOW  - 

_ 

AX  6  k  -0 . 04ff 

1.061  -0.290 

0.094)  ! 

101 

7X  7  1-13.185 

5.108  -0.563 

0.180  ! 

101 

BX  8  !  -2.220 

1.043  -0.311 

0.096  ! 

100 

PIPELINE  CORRELATION  RESULTS! 
DIGITIZING  TO  4  BITS 
ALGORITHM  1  TYPE  2 


RUNNING  AVE  LENGTH  8 
WINDOW 


AX 

6  ! 

-5.279 

3.931  -0.544 

0.184  ! 

98 

7X 

7  ! 

-2.A37 

1.921  -0.295 

0.1«1  ! 

:  104 

8X 

8  ! 

!  -1.207 

1.233  -0.35? 

0.106 

!  103 

( 

PIPELINE  CORRELATION  RESULTS! 
DIGITIZING  TO  2  BITS 
1  ALGORITHM  1  TYPE  1 


RUNNING  AVE  LEN6TH  8 
WINDOW 


AX 

A  ! 

-1.74A 

0.53? 

-0.21? 

0.069  ! 

65 

7X 

7  ! 

9.604 

1.842 

0.400 

0.096  ! 

81 

8X 

8  ! 

9.633 

1.684 

0.35? 

0.087  ! 

105 

( 

( 


L 


PIPELINE  CORRELATION  RESULTS! 

DIGITIZING  TO  2  BITS 
ALGORITHM  1  TYPE  2 

RUNNING  AVE  LENGTH  8 
WINDOW 

AX  A  !  -7.32A  0.7A0  -O.ISO  0.04?  '  57 

7X  7  !  10.303  1.3A8  0.233  O.OAO  !  73 

BX  8  !  -A. 230  0.?05  -0.1?7  0.07!  '  94 

> 

13.4A.5A  >ou»r*  ti»# 

123.70  ARU'S.  .49  CONNECT  HRS 

13.47.02  >run  pIpvc  tntt  iu» 

EXECUTION! 

> 


r 


( 

C 

( 

( 


PIPELINE  CORRELATION  RESULTS! 
DIGITIZING  TO  4  BITS 
ALGORITHH  1  TYPE  I 


RUNNING  AVE  LENGTH  8 

WINDOW  _ _ ^ 

4X  4  i  0.14*  4.345 
7X  7  Pttr023  3.130 
BX  8  ! “11. 080  3.514 


-0.409 

^0.201^1 

10 

-0.520 

0.134  1 

9 

-0.419 

0.148  ! 

10 

C 

c 

c 

c 


PIPELINE  CORRELATION  RESULTS! 
DIGITIZING  TO  4  BITS 
ALGORITHH  1  TYPE  2 


RUNNING  AWE  LENGTH 
WINDOW, 


4X 

4 

r~  -<U01?  ' 

0.299  -0.240 

^01053-'! 

107 

7X 

7 

!  -0.305 

0.349  -0.149 

0.058  ! 

105 

BX 

8 

I  0.139 

0.247  -0.174 

0.048  ! 

107 

PIPELINE  CORRELATION  RESULTS! 

DIGITIZING  TO  4  BITS 
ALGORITHH  1  TYPE  1 

RUNNING  AVE  LENGTH  8 

WINDOW  - — _  T - _ 

4X  4  5.143  -0.483  0.211J  98 

7X  7  !  -4.475  2.298  -0.332  0,132  !  99 

8X  8  !  -4 .458  1.772  -0.411  0.121  !  93 


PIPELINE  CORRELATION  RESULTS! 
DIGITIZING  TO  4  BITS 
ALGORITHH  I  TYPE  2 


RUNNING  AVE  LENGTH  8 
WINDOW 


4X 

4  < 

-7.018 

4.440  -0.414 

0,192 

102 

7X 

7  1 

-0,050 

1.787  -0.421 

0.129 

101 

BX 

8  ! 

4.843 

1.404  0.324 

0.113 

103 

PIPELINE  CORRELATION  RESULTS! 
DIGITIZING  TO  2  BITS 
ALGORITHH  1  TYPE  1 


RUNNING  AVE  LENGTH  8 
WINDOW 


4X 

4  ! 

-4.481 

1.310 

-0.241 

0.092  ! 

45 

7X 

7  ! 

4.445 

2.295 

0.482 

0.140 

81 

8X 

8  ! 

4.757 

1.074 

0.319 

0.078  ! 

105 

PIPELINE  CORRELATION  RESULTS! 
DIGITIZING  TO  2  BITS 
AL80RITHH  1  TYPE  2 


RUNNING  AVE  LENGTH  3 
WINDOW 


4X  4 

1-11.273 

1.841  -0.234 

o.oao  > 

57 

7X  7 

!  7.149 

1.085  0.245 

0.073  ! 

73 

8X  8 

!  -8.805 

1.534  -0.228 

0.087  ! 

94 

PIPELINE  CORRELATION  RESULTS! 

DIGITIZING  TO  4  BITS 
ALGORITHM  1  TYPE  1 

RUNNING  AVE  LENGTH  8 

WINDOW  _ - .  _ _ — - 

AX  A  k  0.02^13. 947  -0.879  0.35!^  89 

7X  7  !-13.S4S  7.17A  -0.A88  0.275  !  91 

8X  8  '-20.701  10.117  -0.719  0.239  1  9A 


PIPELINE  CORRELATION  RESULTS! 
DIGITIZING  TO  A  BITS 
ALGORITHM  1  TYPE  2 


RUNNING  AWE  LENGTH  8 

WINDOW - ^ 

AX  A  0.444  -0 . 2A3  0. OAA>  •  108 

7X  7  I  -2.334  0.784  -0.232  0.082  I  10A 

8X  8  !  -1.34A  0.388  -0.179  0.059  I  107 


PIPELINE  CORRELATION  RESULTS! 
DIGITIZING  TO  4  BITS 
ALGORITHM  1  TYPE  1 

RUNNING  AVE  LENGTH  8 


WINDOW 


Am 

AX 

/UN 

A 

0.04Z  14.349  -0.950 

'0.389>  ! 

84 

7X 

7 

-21.739  9.231  -0.729 

0.212  ! 

92 

8X 

8 

-20.951  8.033  -0.A79 

0.192  ! 

92 

PIPELINE  CORRELATION  RESULTS! 

DIGITIZING  TO  4  BITS 
ALGORITHM  1  TYPE  2 

RUNNING  AVE  LENGTH  8 
WINDOW 

AX  A  !  -5.012  2.093  -0.407  0.131  I  103 

7X  7  !  -7.201  5.533  -0.A51  0.21A  !  101 

BX  8  I  -A. 738  0.938  -0.2A9  0.0A9  I  101 


PIPELINE  CORRELATION  RESULTS! 

DIGITIZING  TO  2  BITS 
ALGORITHM  1  TYPE  1 

RUNNING  AVE  LENGTH  8 
WINDOW 

AX  A  1-21.098  7.285  -0.518  0.171  !  A3 

7X  7  I-11.13A  5.259  -0.53A  0.195  !  79 

8X  8  I  -7.925  1.871  -0.330  0.109  I  104 


PIPELINE  CORRELATION  RESULTS! 
DIGITIZINO  TO  2  BITS 
ALGORITHM  1  TYPE  2 


RUNNING  AVE  LENGTH  8 
WINDOW 


AX 

A  ! -20.SA1 

7.03A  -O.AOA 

0.170  ! 

57 

7X 

7  1-10.513 

3.405  -0.37A 

0.143  ! 

73 

8X 

8  1-22.348 

4.484  -0.557 

0.138  ! 

94 

> 


14.3A.2A  >lo« 

570.78  ARU'S.  1.27  CONNECT  HRS 
LOGGED  OPP  AT  14.3A.30  ON  03HAR82H 


c 


rjruiiu  bvnntkn  t  a  uif  nc.auk.iat 

DIGITIZING  TO  6  IITS 
ALGORITHM  1  TYPE  1 


RUNNING  AVE  LENGTH  G 


6  L-0.ili>  4G.2SG_zlUllS^*uI*Ll _ Zp 

7X  7  !-60. 552  46.101  -1.23?  0.362  I  77 
»X  a  1-51. 8?3  35.P05  -1.201  0.31?  •  82 


PIPELINE  CORRELATION  RESULTS* 
DIGITIZING  TO  6  IITS 
ALGORITHM  1  TYPE  2 


RUNNING  AVE  LENGTH  S 


HJNOOIU-^-^  - 

<6X  6  vS*gQ.12£>- 


JL05^ 


} -0*707— 0*346  *'0.081')!  101 
7X  7  I  -3.565  0.632  -0.248  0.06?  !  106 
BX  8  I  0.421  0.670  -0.230  0.07?  !  106 


PIPELINE  CORRELATION  RESULTS.* 
DIGITIZING  TO  4  BITS 
ALGORITHM  1  TYPE  1 


c 

( 


RUNNING  AVE  LENGTH  6 

WINDOW— - - T-5-^c 

^-6X  58.753  -1J22_^*ZJ£)L-71 

7X  7  1-51.495  32.383  -1.115  0.26?  !  83 
•X  8  !-57.0?8  38.5??  -0.?34  0.27?  !  82 


PIPELINE  CORRELATION  RESULTS! 
DIGITIZING  TO  4  BITS 
ALGORITHM  1  TYPE  2 


c 

( 

c 

c 

( 

( 


RUNNING  AVE  LENGTH  8 
WINDOW  . 

AX  6  ! -11 .836  5.645  -0.727  0.201  • 

7X  7  1-25.187  ?.223  -0.756  0.176  ! 

BX  8  1-11. ?27  3.28?  -0.4?5  0.135  1 


PIPELINE  CORRELATION  RESULTS! 
DIGITIZING  TO  2  BITS 
ALGORITHM  1  TYPE  1 

RUNNING  AVE  LENGTH  8 
WINDOW 


100 

?4 

100 


6X  6 

1-51.472 

31.703  -0.882 

0.258  ! 

57 

7X  7 

1-28.751 

15.008  -0.864 

0.261  ! 

75 

•X  8 

1-26.800 

8.758  -0.630 

0.140  ! 

103 

'  1 

PIPELINE  CORRELATION  RESULTS! 
DIGITIZING  TO  2  BITS 
ALGORITHM  1  TYPE  2 


RUNNING  AVE  LENGTH  8 
WINDOW 


AX 

6 

1-47.542 

26.756  -0.886 

0.232  < 

57 

7X 

7 

1-34.170 

15.732  -0.772 

0.213  1 

73 

IX 

1 

!-I?.4SG 

6.554  -0.644 

0.168  ! 

?4 

> 


STAGE  III 


execution: 

>.05>.lr0.»58> 

,  THE  SCENES  HERE  GENERATED  WITH  0.050  PROBABILITY  OF  EDGE 

V  RELATIVE  SHIFT  IS  0.100  PIXELS 

HEAN  INTENSITY  IS  0.491 

STD  DEV  IS  0.198 
C  1 • /SNR  IS  0.0 


PIPELINE  CORRELATION  RESULTS! 
DIGITIZING  TO  4  BITS 
1  FOR  100  LINES 

RUNNING  AVE  LENGTH  8 
C  UINDOU 


6X 

6  ! 

-0.008 

0.034  -0.103 

0.033  1 

107 

7X 

7  ! 

-0.001 

0.009  -0.024 

0.009  ! 

107 

8X 

8  ( 

!  0.000 

0.007  -0.024 

0.007  1 

1  107 

PIPELINE  CORRELATION  RESULTS: 
DIGITIZING  TO  4  BITS 
C  FOR  100  LINES 

RUNNING  AVE  LENGTH  8 
C  UINDOU 


6X 

6  ! 

-0.002 

0.049 

-0.126 

0.049  1 

107 

7X 

7  1 

0.001 

0.014 

0.039 

0.014  i 

107 

8X 

8  1 

0.004 

0.010 

0.035 

0.009  ! 

!  107 

PIPELINE  CORRELATION  RESULTS! 
DIGITIZING  TO  2  BITS 

<  FOR  100  L'VES 

RUNNING  AVE  LENGTH  8 

<  UINDOU 


6X 

6  1 

-0.093 

0.094  -0.100 

0.013  ! 

97 

7X 

7  ! 

i  -0  088 

0.089  -0.100 

0.016  1 

103 

8X 

8  i 

-0.088 

0.089  -o.y>o 

0.015  I 

107 

(  .05. .3.0. >589> 

THE  SCENES  HERE  GENERATED  UITH  0.050  PROBABILITY  OF  EDGE 
RELATIVE  SHIFT  IS  0.300  PIXELS 
C  HEAN  INTENSITY  IS  0.553 
STD  DEV  IS  0.181 
1 ./SNR  IS  0.0 


V  PIPELINE  CORRELATION  RESULTS! 

DIGITIZING  TO  6  BITS 
FOR  100  LINES 

C 

RUNNING  AVE  LENGTH  8 
UINDOU 

6X  6  !  -0.024  0.042  -0.101  0.035  I  107 

„  7X  7  !  -0.015  0.051  -0.171  0.049  !  107 

8X  8  !  -0.012  0.032  -0.095  0.030  !  100 


l  PIPELINE  CORRELATION  RESULTS! 

DIGITIZING  TO  4  BITS 
FOR  100  LINES 

RUNNING  AVE  LENGTH  8 
UINDOU 

&X  6  i  -0.024  0.044  -0.117  0.037  I  107 

7X  7  •  -0.020  0.0S4  -0.194  0.051  !  100 

•X  8  !  -0.012  0.034  -0.113  0.031  I  106 


C  PIPELINE  CORRELATION  RESULTS! 

DIGITIZING  TO  2  BITS 
FOR  100  LINES 

c 

RUNNING  AVE  LENGTH  6 
WINDOW 

<f  iX  i  !  -0.092  0.140  -0.358  0.105  !  93 

7X  7  f  -0 • 044  0.119  -0.324  0.100  !  99 
8X  8  I  -0.044  0.105  -0.295  0.083  !  100 

(f  > 


e 

•  OS  * ♦ 5  9  0 • ? 43  9 

((  THE  SCENES  WERE  GENERATED  WITH  0.050  PROBABILITY  OF  EDGE 
RELATIVE  SHIFT  IS  0.500  PIXELS 
MEAN  INTENSITY  IS  0.534 
(r  STD  DEV  IS  0.195 
1 • /SNR  IS  0.0 

(< 

PIPELINE  CORRELATION  RESULTS! 

0  DIGITIZING  TO  6  BITS 
FOR  100  LINES 

RUNNING  AVE  LENGTH  8 


WINDOW 


6X 

6  < 

-0.020 

0.048 

0.114 

0.043  1 

!  106 

7X 

7  ! 

0.008 

0.025 

0.085 

0.024  1 

107 

ax 

8  ! 

!  0.005 

0.018 

0.066 

0.017  1 

!  107 

PIPELINE  CORRELATION  RESULTS! 
DIGITIZING  TO  4  BITS 
FOR  100  LINES 

l  RUNNING  AVE  LENGTH  8 


WINDOW 


6X 

6  < 

-0.020 

0.047 

0.113 

0.043  1 

107 

7X 

7  \ 

0.005 

0.027 

0.089 

0.026  1 

I  107 

BX 

8  1 

0.003 

0.019 

0.069 

0.019  1 

i  107 

..  PIPELINE  CORRELATION  RESULTS! 

DIGITIZING  TO  2  BITS 
FOR  100  LINES 

AVE  LENGTH  8 

-0.137  0.181  -0.446  0.117  !  95 

-0.078  0.130  -0.318  0.104  j  106 

-0.073  0.100  -0.255  0.068  I  106 


•  0Sr.7»0.>  753 » 

THE  SCENES  WERE  GENERATED  WITH  0.050  PROBABILITY  OF  EDGE 
(  RELATIVE  SHIFT  IS  0.700  PIXELS 

HEAN  INTENSITY  IS  0.493 

STD  DEV  IS  0.223 
(.  1  ./SNR  IS  0.0 


V  RUNNING 
WINDOW 
.  6X  6  I 

L  7X  7  I 

sx  e  i 

> 


PIPELINE  CORRELATION  RESULTS! 
DIGITIZING  TO  6  BIT8 
C  FOR  100  LINES 

RUNNING  AVE  LENGTH  8 
t  WINDOW 


6X 

6  • 

-0.036 

0.109 

-0.282 

0.103  ! 

107 

7X 

7  1 

-0.002 

0.053 

-0.123 

0.053  ! 

107 

ax 

S  ! 

>  0.008 

0.039 

0.115 

0.038  ! 

!  107 

PtPFITNC  CORRELATION  RESULTS! 


PIPELINE  CORRELATION  RESULTS! 

DIGITIZING  TO  4  BITS 
POR  100  LINES 

RUNNING  AVE  LENGTH  8 
WINDOW 

AX  6  I  -0.023  0.10S  -0.282  0.102  I  107 

7X  7  !  0.007  0.044  -0.126  0.044  !  107 

BX  8  !  0.013  0.036  0.107  0.033  f  107 


PIPELINE  CORRELATION  RESULTS! 

DIGITIZING  TO  2  BITS 
FOR  100  LINES 

RUNNING  AVE  LENGTH  8 
WINDOW 

6X  6  I  -0.135  0.196  -0.504  0.143  !  97 

7X  7  I  -0.103  0.167  -0.481  0.131  !  101 

8X  8  I  -0.088  0.148  -0.517  0.119  !  105 


,05.1.1.0. ,32. 

THE  SCENES  WERE  GENERATED  WITH  0.050  PROBABILITY  OF  EDGE 
RELATIVE  SHIFT  IS  1.100  PIXELS 
HEAN  INTENSITY  IS  0.519 
STD  DEV  IS  0.242 
1 • /SNR  IS  0.0 


PIPELINE  CORRELATION  RESULT^! 

DIGITIZING  TO  6  BITS 
FOR  100  LINES 

RUNNING  AVE  LENGTH  8 
WINDOW 

6X  6  I  -0.059  0.095  -0.203  0.074  I  107 
7X  7  !  -0.025  0.058  -0.185  0.052  !  107 


8X  8  I  -0.019  0.051  -0.190  0.047  I  107. 


PIPELINE  CORRELATION  RESULTS! 

DIGITIZING  TO  4  BITS 
FOR  100  LINES 

RUNNING  AVE  LENGTH  8 
WINDOW 

6X  6  •  -0.065  0.099  -0.236  0.075  !  107 

7X  7  !  -0.022  0.065  -0.191  0.061  !  107 

8X  8  !  -0.016  0.055  -0.198  0.052  !  107 


PIPELINE  CORRELATION  RESULTS! 

DIGITIZING  TO  2  BITS 
FOR  100  LINES 

RUNNING  AVE  LENGTH  8 
WINDOW 

6X  6  ‘  -0.146  0.229  -0.556  0.176  !  91 

7X  7  >  -0.083  0.165  -0.431  0.142  !  103 

SX  8  I  -0.058  0.134  -0.278  0.121  I  107 


.05. .9.0. ,98754» 

THE  SCENES  WERE  GENERATED  WITH  0.050  PROBABILITY  OF  EDGE 

RELATIVE  SHIFT  IS  0.900  PIXELS 

HEAN  INTENSITY  IS  0.483 

STD  DEV  IS  0.211 

l./SNR  IS  0.0 


PIPELINE  CORRELATION  RESULTS: 

DIGITIZING  TO  6  BITS 
FOR  100  LINES 

RUNNING  AVE  LENGTH  8 
WINDOW 

AX  6  I  -0.101  0.159  -0.357  0.123  f  106 

7X  7  I  -0.008  0.033  -0.122  0.032  I  107 

8X  8  I  -0.010  0.042  -0.159  0.041  *  107 


PIPELINE  CORRELATION  RESULTS: 

DIGITIZING  TO  4  BITS 
FOR  100  LINES 

RUNNING  AVE  LENGTH  8 
WINDOW 

AX  6  I  -0.102  0.166  -0.335  0.131  !  105 

7X  7  !  -0.006  0.036  -0.129  0.036  *  107 

8X  8  I  -0.010  0.049  -0.184  0.048  !  107 


PIPELINE  CORRELATION  RESULTS! 
DIGITIZING  TO  2  BITS 
FOR  100  LINES 


RUNNING  AVE  LENGTH  8 
MXNDOr 


6X 

6  ! 

-0.185 

0.241  -0.507 

0.154  1 

!  91 

7X 

7  ! 

-0.030 

0.152  -0.409 

0.149  1 

!  100 

BX 

8  ! 

-0.025 

O.lOO  -0.243 

0.097  I 

!  104 

> 


14.59.11  > 1 oN 

188.39  ARU'S.  .25  CONNECT  HRS 
LOGGED  OFF  AT  14.59.14  ON  05HAR82I 


STAGE  IV 


execution: 

> ■ 05  .  0 »  » *  1 . 3435  > 

THE  SCENES  WERE  GENERATED  WITH  0.050  PROBABILITY  OF  EDGE 

RELATIVE  SHIFT  IS  0.0  PIXELS 

MEAN  INTENSITY  IS  0.409 

STD  DEV  IS  0.257 

l./SNR  IS  0.100 


PIPELINE  CORRELATION  RESULTS: 

DIGITIZING  TO  4  BITS 
FOR  100  LINES 

RUNNING  AVE  LENGTH  8 
WINDOW 

9X  9  !  -0.002  0.029  -0.082  0.029  I  107 

>.05.08 .102. *1.3435. 

THE  SCENES  WERE  GENERATED  WITH  0.050  PROBABILITY  OF  EDGE 

RELATIVE  SHIFT  IS  0.200  PIXELS 

MEAN  INTENSITY  IS  0.409 

STD  DEV  IS  0.251 

l./SNR  IS  0.100 


PIPELINE  CORRELATION  RESULTS: 

DIGITIZING  TO  4  BITS 
FOR  100  LINES 

RUNNING  AVE  LENGTH  8 
WINDOW 

9X  9  !  0.002  0.026  0.065  0.026  !  107 

> 


•05. .4. . 1 >79865. 

THE  SCENES  WERE  GENERATED  WITH  0.050  PROBABILITY  OF  EDGE 

RELATIVE  SHIFT  IS  0.400  PIXELS 

MEAN  INTENSITY  IS  0.542 

STD  DEV  IS  0.193 

l./SNR  IS  0.100 


PIPELINE  CORRELATION  RESULTS: 

DIGITIZING  TO  4  BITS 
FOR  100  LINES 

RUNNING  AVE  LENGTH  8 
WINDOW 

9X  9  !  -0.003  0.057  0.149  0.039  I  107 

>.0S..6> .1>2149> 

THE  SCENES  WERE  GENERATED  WITH  0.050  PROBABILITY  OF  EDGE 

RELATIVE  SHIFT  IS  0.600  PIXELS 

MEAN  INTENSITY  IS  0.518 

STD  DEV  IS  0.215 

l./SNR  IS  0.100 


PIPELINE  CORRELATION  RESULTS: 

DIGITIZING  TO  4  BITS 
FOR  100  LINES 

RUNNING  AVE  LENGTH  8 
WINDOW 

9X  9  !  -0.012  0.046  -0.139  0.037  !  106 


.05' *  8 ' *  1 r  6551 1 

THE  SCENES  MERE  GENERATED  WITH  0.050  PROBABILITY  OF  EDGE 
RELATIVE  SHIFT  IS  0.800  PIXELS 
HEAN  INTENSITY  IS  0.525 
STD  DEV  IS  0.198 
1 » /SNR  IS  0.100 


PIPELINE  CORRELATION  RESULTS) 

DIGITIZING  TO  4  BITS 
FOR  100  LINES 

RUNNING  AVE  LENGTH  8 
WINDOW 

9X  9  I  -0.024  0.047  0.084  0.029  I  107 

> 


•  05  fl.'.l'41S'7f 

THE  SCENES  WERE  GENERATED  UITH  0.050  PROBABILITY  OF  EDGE 
RELATIVE  SHIFT  IS  1.000  PIXELS 
HEAN  INTENSITY  IS  0.484 
STD  DEV  IS  0.210 
1 • /SNR  IS  0.100 


PIPELINE  CORRELATION  RESULTS) 

DIGITIZING  TO  4  BITS 
FOR  100  LINES 

RUNNING  AVE  LENGTH  8 
WINDOW 

9X  9  !  -0.050  0.105  -0.358  0.040  I  107 

FROM  RJESCED  1  Job  cannot  run  until  you  ralaaia  CAISOFTU' 

> 


. 05 ' 0 . ' ,2.54546' 

C  THE  SCENES  WERE  GENERATED  WITH  0.050  PROBABILITY  OF  EDGE 
RELATIVE  SHIFT  IS  0.0  PIXELS 
HEAN  INTENSITY  IS  0.448 
C  STD  DEV  IS  0.212 
1 , /SNR  IS  0.200 

( 

PIPELINE  CORRELATION  RESULTS) 

*.  .  DIGITIZING  TO  4  BITS 
FOR  100  LINES 

L  RUNNING  AVE  LENGTH  8 
WINDOW 

9X  9  !  -0.002  0.048  0.134  0. 01/  I  107 


( 

< 

c 

c 

c 

c 

c 

< 

c 

c 

c 

( 

( 

c 

( 

( 

( 

( • 


.05,1.,. 2. 1557, 

THE  SCENES  WERE  GENERATED  WITH  0.050  PROBABILITY  OF  EDGE 
RELATIVE  SHIFT  IS  1.000  PIXELS 
HEAN  INTENSITY  IS  0.528 
STD  DEV  IS  0.185 
1 ./SNR  IS  0.200 


PIPELINE  CORRELATION  RESULTS  1 
DIGITIZING  TO  4  BITS 
FOR  100  LINES 

RUNNING  AVE  LENGTH  8 
WINDOW 

9X  9  I  -0.091  0.120  -0.299  0.055  !  107 

> 


•05,0. ,.4,43, 

THE  SCENES  WERE  GENERATED  WITH  0.050  PROBABILITY  OF  EDGE 
RELATIVE  SHIFT  IS  0.0  PIXELS 
HEAN  INTENSITY  IS  0.483 
STD  DEV  IS  0.207 
1 ./SNR  IS  0.400 


PIPELINE  CORRELATION  RESULTS! 
DIGITIZING  TO  4  BITS 
FOR  100  LINES 


RUNNING  AVE  LENGTH  8 
WINDOW 

9X  9  !  0.015  0.100  -0.234 

> 


0.048  I  107 


.05, 1 . , .4,453, 

THE  SCENES  WERE  GENERATED  WITH  0.050  PROBABILITY  OF  EDGE 
RELATIVE  SHIFT  IS  I. 000  PIXELS 
HEAN  INTENSITY  IS  0.489 
STD  DEV  IS  0.232 
1 , /SNR  IS  0.400 


»  PIPELINE  CORRELATION  RESULTS! 

DIGITIZING  TO  4  BITS 
FOR  100  LINES 

, 

RUNNING  AVE  LENGTH  8 
WINDOW 

9X  9  t  -0,192  0.223  -0.429  0.045  !  107 

> 


•  05 1 •  8  f . 4  »  45 • 

THE. SCENES  WERE  GENERATED  WITH  0.050  PROBABILITY  OP  EDGE 
RELATIVE  SHIFT  IS  0.800  PIXELS 
MEAN  INTENSITY  IS  0.483 
STD  DEV  IS  0.171 
1 . /SNR  IS  0.400 


PIPELINE  CORRELATION  RESULTS t 
DIGITIZING  TO  4  BITS 
FOR  100  LINES 

RUNNING  AVE  LENGTH  8 
UINDOU 

9X  9  I  -0. 176  0.197  -0.352  0.059  !  107 

> 


•  05  f .If . 4  f  6543  f 


THE  SCENES  MERE  GENERATED  WITH  0.050  PROBABILITY  OF  EDGE 

RELATIVE  SHIFT  IS  0.400  PIXELS 

MEAN  INTENSITY  IS  0.543 

STD  DEV  IS  0.227 

l./SNR  IS  0.400  • 


PIPELINE  CORRELATION  RESULTS t 
DIGITIZING  TO  4  BITS 
FOR  100  LINES 

RUNNING  AVE  LENGTH  8 
UINDOU 

9X  9  !  -0.073  0.115  -0.308  0.055  I  107 

> 


t 


. 05 f  *4f  »  4  f 4574  f 


THE  SCENES  WERE  GENERATED  UITH  0.050  PROBABILITY  OF  EDGE 

RELATIVE  SHIFT  IS  0.400  PIXELS 

MEAN  INTENSITY  IS  0.424 

STD  DEV  IS  0.241 

l./SNR  IS  0.400 


PIPELINE  CORRELATION  RESULTS) 
DIGITIZING  TO  4  BITS 
FOR  100  LINES 


RUNNING  AVE  LENGTH  8 
WINDOW 

9X  9  !  -0.071  0.124  -0.239  0.049  !  107 


•  OS • • 2  > .  4 . 654  f 

THE  SCENES  WERE  GENERATED  WITH  O.OSO  PROBABILITY  OF  EDGE 

RELATIVE  SHIFT  IS  0.200  PIXELS 

MEAN  INTENSITY  IS  0.485 

STD  DEV  IS  0.179 

l./SNR  IS  0.400 


PIPELINE  CORRELATION  RESULTS’. 

DIGITIZING  TO  4  BITS 
FOR  100  LINES 

RUNNING  AVE  LENGTH  8 
WINDOW 

9X  9  I  -0.044  0.108  -0.234  0.048  t  107 

> 

FROM  RJESCED  S  Job  cannot  run  until  uou  release  CAISOFTW's  disks. 
18.11.57  >  . 


STAGE  V 


(  execution: 

>•05,0. , . 1  ,3435, 

r  THE  SCENES  HERE  GENERATED  WITH  0.050  PROBABILITY  OF  EDGE 
'  RELATIVE  SHIFT  IS  0.0  PIXELS 

MEAN  INTENSITY  IS  0.409 

STD  DEV  IS  0.25? 

*  I • /SNR  IS  0.100 


( 

PIPELINE  CORRELATION  RESULTS! 

DIGITIZING  TO  4  BITS 
(  FOR  100  LINES 

RUNNING  AVE  LENGTH  8 
(  WINDOW 

9X  9  1  -0.002  0.029  -0.082  0.029  I  107  I  -0.001  0.027  -0.064  0.027  I  107 

> 

( 


( 


c 


c 

.05.. 8. .1,4551. 

C  THE  SCENES  WERE  GENERATED  WITH  0.050  PROBABILITY  OF  EDGE 
RELATIVE  SHIFT  IS  0.800  PIXELS 
,  MEAN  INTENSITY  IS  0.S25 

V  STD  DEV  IS  0.198 

l./SNR  IS  0.100 

( 

PIPELINE  CORRELATION  RESULTS! 

C  DIGITIZING  TO  4  BITS 
FOR  100  LINES 

C  RUNNING  AVE  LENGTH  8 
WINDOW 

9X  9  I  -0.024  0.043  0.084  0.034  !  107  I  -0.021  0.052  -0.110  0.047  I  107 

(  >.05,1. , .1,415, 


( 


c 

THE  SCENES  WERE  GENERATED  WITH  0.050  PROBABILITY  OF  EDGE 
'  RELATIVE  SHIFT  IS  1.000  PIXELS 
MEAN  INTENSITY  IS  0.4B4 
STD  DEV  18  .a-ZLO 
C  l./SNR  IS  0.100  \ 


c 

PIPELINE  CORRELATION  RESULTS! 

.  DIGITIZING  TO  4  BITS 
V  FOR  100  LINES 

RUNNING  AVE  LENGTH  8 
C  WINDOW 

9X  9  I  -0.050  0.105  -0.358  0.092  I  107  !  -0.054  0.128  -0.520  0.116  I  107 

> 

(. 


c 

.05,  .8,  .4,45, 

THE  SCENES  WERE  GENERATED  WITH  0.050  PROBABILITY  OF  EDGE 
(  RELATIVE  SHIFT  IS  0.800  PIXELS 
MEAN  INTENSITY  IS  0.483 
STD  DEV  IS  0.171 
C  l./SNR  IS  0.400 


l 

PIPELINE  CORRELATION  RESULTS! 

DIGITIZING  TO  4  BITS 
l  FOR  100  LINES 

RUNNING  AVE  LENGTH  ■ 
l  WINDOW 

9X  9  I  -0.174  0.197  -0.352  0.0S8  I  107  !  -0.144  0.187  -0.339  0.091  I  107 

> 

l 


C 


c 

.03. .0. .4.4343. 

THE  SCENES  WERE  GENERATED  WITH  0.030  PROBABILITY  OF  EDGE 
(  RELATIVE  SHIFT  IS  0.400  PIXELS 
MEAN  INTENSITY  IS  0.543 
STO  DEV  IS  0.227 
(  t./SNR  IS  0.400 


PIPELINE  CORRELATION  RESULTS: 

DIGITIZING  TO  4  BITS 
(  FOR  100  LINES 

RUNNING  AVE  LENGTH  G 
(  WINDOW 

9X  9  I  -0.073  0.115  -0.308  0.088  I  107  !  -0.071  0.112  -0.311  0.084  i  107 

> 

c 


t 


,  /t. 03, 1, . ,4.433. 

V.  THE  SCENES  WERE  GENERATED  WITH  0.030  PROBABILITY  OF  EDGE 
RELATIVE  SHIFT  IS-4...0.00'  PIXELS 
,  HEAN  INTENSITY  IS  0.489 

V.  STD  DEV  IS  0.232 

l./SNR  IS  0.400 

C  J 

PIPELINE  CORRELATION  RESULTS! 

'  DIGITIZING  TO  4  BITS 
FOR  100  LINES 

(  RUNNING  AVE  LENGTH  8 

WINDOW 

9X  9  f  -0.192  0.223  -0.429  0.114  I  107  !  -0.183  0.213  -0.429  0.113  I  107 

(  > 


( 


.05. .4. .4,4574. 

(  THE  SCENES  WERE  GENERATED  WITH  0.030  PROBABILITY  OF  EDGE 
RELATIVE  SHIFT  IS  O.V>°  PIXELS 
MEAN  INTENSITY  IS  0.424 
'  STD  DEV  IS  0*241 
l./SNR  IS  0.400 


PIPELINE  CORRELATION  RESULTS! 

(  DIGITIZING  TO  4  BITS 
FOR  100  LINES 

(  RUNNING  AVE  LENGTH  8 

WINDOW 

9X  9  I  -0.071  0.124  -0.239  0.104  1  107  !  -0.049  0.118  -0.237  0.093  I  104 

<  > 


.05, 0 .  v . 4, 43 ,  /- 

THE  SCENES  WERE  GENERATED  WITH  0.050  PROBABILITY  OF  EDGE 

RELATIVE  SHIFT  IS  0.0  PIXELS 

HEAN  INTENSITY  IS  0.483 

STD  DEV  IS  0.207 

l./SNR  IS  0.400 


PIPELINE  CORRELATION  RESULTS! 

DIGITIZING  TO  4  BITS 
FOR  100  LINES 

RUNNING  AVE  LENGTH  8 
WINDOW 

9X  9  !  0.015  0.100  -0.234  0.099  I  107  !  0.017  0.102  -0.234  0.100  I  107 

> 


STAGE  I 


c  program:  pipeline  correlation 

C  CAI  DIU1S0N  RECON  OPTICAL 

C  0.  DEFOE 

C  DATE!  10  FEBRUART  1982 

C 

INTEGER  SEED 

real  ns 

DIMENSION  PK7.1I).  P2*7» 1 1 ) »  ILN( 10.3.2. A. 5) i 
>IASN*3.2.4.5> .  IASD*3.2.4.3>.  URI8.3.2.4.5).  L<4>. 
>EA*J.2.4.5.4>.  E2(3.2.4<5.4>»  EHX<3.2.4.5.4>.  K0<3>2>» 
>KD*4).  KF ( 4  > »  SF(4)>  ILD( 10.3.2.4.5) .  FRI3.2.4.3.4) > 
>SDF(4).PL1( IQ)  .PL2* 10)  »RKF<4) »KE(3.2.4.5) 

BATA  KD/10.8.4.4.2.1/.  KF/1.2.4.B/ 

NH-7)  NL-100*  MU-81  NAL>3I  NAD*2I  NO-41  NF-4 
JUIN-41  JMIX-10 
CM*. 23*  M*.02l  CMN-4. 

NA-JUIX*  NN-JUIXM)  NMI-JU1X-JNINM 
RH-5QRT* .5) 

KS1-NH4NA-1I  KS2-KS14NU-1 

NMH-NH-1*  NMP-NH.ll  NLM-ML4KS2 

NNM-NN-1)  NNP-NN4I*  NAP-NA41 

NMN-N8-I*  NMP-NM+ll  NMIN-NMI-1*  NUIP-NMI+1 

JUINM-JMIN-1 

NNN-CNHFM*  WN-WNH+WNHI  8NP-1.-8N 

DO  8  IN-l.NNI 

BO  7  ID-l.ND 

DO  4  IAD-l.NAD 

BO  3  IA-1 »NAL 

KECIA.IAD.I0.I8I-0 

3  CONTINUE 

4  CONTINUE 

7  CONTINUE 

8  CONTINUE 

DO  9  IF-l.NF 

9  RKF(IF)-1./TCF*IF) 

1  SP-SP2-0.I  NP-NE-0 

READ(5.5000)  PE.D1 -NS, SEED 
5000  FORMATUFIO.  4.110) 

IF(PE.EQ.O)  00  TO  9999 

D2-1-D1  I  C*(.5-ABS(.5-D1))*CM  I  NP-NH 

NPM-MP-1 

MA-HA  I  MU-MU 

DO  499  M-l.MLH 

MPN-NP  I  MP*MP4t 
IFIHP.EO.NMP)  HP-1 
IF(H.LT.NH)  00  TO  10 
HA-MA41 

IF (MA.EQ.NAP)  HA-1 
IF(N.LT.KSl)  80  TO  10 
MU-HH41 

ifinm.eo.nupi  mw-1 

ME-NE41 

10  DO  199  NR-lfNM 

N-NNP-NR  I  U-RAN<SEED)  . 

IF(U.BE.PE)  80  TO  20 

IU-S 

80  TO  30 

20  U-RAN*SEED1  »  IU-48U41  I  L<l*-0 

30  IFCIU.8E.3)  L(l)-t 

LC2)-LCI1-LC4)-L<1* 

IFdU.OE.4)  80  TO  100 
U-RAN(SEED) 

IF(IU.E8.2>  80  TO  SO 
JI1-48U41 

IFIIO.NE.il  80  TO  40 
LCJUfl 
00  TO  100 
40  LUU1-0 
80  TO  100 
30  JO-488 

IF* JU.0C.3)  80  TO  40 
LCIl-1  I  LTJ0421-1 
80  TO  100 
40  UJU-D-l 

IF* JU.EO.S)  JO-2 
LIJOl-l 

100  T1-T2-T3-T4-.3 

IF(N.LT.NN)  T1-P1INP.N4C ) 

IFIN.E8.il  80  TO  110 
Tl-PltHPN»N) 

IFCN.LT.NN1  T2-P1CHPH.N41) 

IFCN.8T.il  T4-P1*hPN.N-1) 

110  NA-CT18LCI )IT38LC3}4(T28LC214T48LC4))8RM)/CLC114L(319 
>CLC219LC4118RCC) 

1FC IU.E0.3)  80  TO  120 
8-1080  I  8-8A088  rsEED.IMfO)  1  80  TO  150 
120  U-RAN<8CE81 

IFCMA.LT.ON)  80  TO  130 
1FCMA.8T.NNP1  80  TO  t40 
0-8NP88 


130  8-(0A4UNH)8(l.-U)4U 

GO  TO  150 
140  0-<UA-UNH)8U 
ISO  P1(MP,N)»Q  I  SP-SP40 
SP2-SP24Q80  )  NP-NPil 
IF(N.EQ.l)  GO  TO  199 
Tl-Pl(HPM.N) 

T3-PKMP.N) 

T2-T1-T3 

8-TI-D18T2 

P2(HPNiN)-Q  I  SP-SP40 
8P2-SP2+08Q  1  NP-NP41 

199  CONTINUE 

4000  F0RNAT(211X>10(FS.3.1X>)1 

200  IFCN.LT.NM)  00  TO  499 

SIP-SORT*  *SP2-SP8SP/NP)/NP)  1  SIN-NS8SIP 
IFCH.GT iNHl  GO  TO  210 
DO  209  HM-l.NHM 
DO  208  N-liNNH 

Tl-Pt*HH.N>t8AUSS<SEED.0..SIN> 
T2-P2*HH.NH8AU8S(SEED.O.  .SIN) 

Tt-AHAXl (AHlNl(Tl.l.)iO.) 

T2-AHAX1  (AHINt  (T2. 1 . )  rO.  1 
PI (HH. NJ-PL1 (N)-Tl 
P2(MH.N)-PL2(N)-T2 

208  CONTINUE 

CALL  IMAGE* 10. PL1.PL2) 

209  CONTINUE 
GO  TO  220 

210  DO  219  N-l.NNN 

Tl-PKMPM .  N)  tGAUSS  ( SEED .  0 . .  SIN  > 

T 2-P2  ( MPM » N > * G AUS S ( SEE D . 0 .  .SIN) 

Tt-AHAXlt  AHINHT1 . 1 .  )  »0. 1 
T2-AMAX1  (AMIN1  (T2. 1 .1*0.) 

PI *MPM.N)-PL1 (N)-Tl 
P2(MPM.N)-PL2(N)-T2 
219  CONTINUE 

CALL  IMAGE*  10. PL1.PL2) 

220  DO  498  IA-1.NAL 


80  497  IAD-1.NAD 


80  494  ID-l.ND 
XBH-KDIID)  I  KS-288KSH 
IF(N.LE.KSl)  80  TO  230 
NAM-MA4JU1HM 
IF(HAN.OT.NA)  MAM-MAM-NA 
DO  229  18-1. NNI 
NAN-MAM- 1 

IF(MAM.LE.O)  MAN-MAM4NA 

IABN*IA.IAD.IB.IN)-IASN*IA.IAD.ID.IU)-ILN(MAM.IA.IAD.ID.IU) 
IA88* IA> IA8. ID. IMI-IASD* IA.1AD. ID. IMl-ILDCMAN.IA. IA8fI8.INI 
229  CONTINUE 
230  I81-I82-0 
.  DO  249  N-I.JMIN 


249  CONTINUE 

ILN*NA.IA.IAD. 18.11-181 
IL8CI1A.IA.IA8.  ID.  11-182 
TFINNI.EO.l)  80  TO  240 
80  239  IMH-I.NNIN 
N-JUIN4INN  I  I8-I8M41 


ILN*MA«tA.IA8.tD.IU)-I81 

ILB(MA.IA.IA8.ID.IH)-IS2 

239  CONTINUE 

240  IF(M.LT.N81)  80  TO  494 
IF*N.8T.M81)  80  TO  200 
80  279  IMH-t.NMI 
ltt-132-0 
IU-NNIP-INN 

80  249  MAN-INN. JOIN 

Ifl-Ill4ILNINAN.IA.IABtl8.I8l 

II2-l324tL8CHAN.tA.IAB.IB.INI 


1FCU.8T.MA1  8-848N 
80  TO  130 


24?  CONTINUE 

IASN(IA.IAB.ID.IH>*1S1 
IASD(IA.IAD,ID,IW)-IS2 
27?  CONTINUE 
GO  TO  300 

200  BO  20?  IH>1.NHI 

IASN(IA.IAD.ID,XU)*1ASN(IA.IAD»ID»IUHXLN(HA.IA»1AD»ID.IU) 

IASD(IA.IAD,I&,IU)*1ASD(IA,IAD,ID,IU)41L&(MA,1A,IAD,1D,XU) 

28?  CONTINUE 

2?0  IF(H.LE.KS2>  00  TO  300 
80  2??  1U*1 ,NNX 
BO  2??  IF.l.NF 
NH-HU-KF(IF>»  RF-AKF(IF) 

IF(HH.LE.O)  HK»MH4NU 

FR(IA.IAD.ID.iy.IF>*FR<IA,IAD,IB,IU,IF>-yR<HH,IA,IAD.ID.IU>BRF 
2??  CONTINUE 

300  BO  4?3  IUM.NHI 

T1*IASN(XA,1AD,XD,IU)  I  T2*IASD(IA,IAD,ID,iy> 

T3*A»S(T2> 

IF(T3.BT..3BABS(T1).AND.T3.BE.S.)  00  TO  301 
GO  TO  302 

301  KE(IA,IA0.ID,iy)*KE(IA,IAD.ID,IU)41 

302  UR(HU.IA.XAD.ID.IU) 

IHN.l-T.K82)  60  TO  4?5 
11*0)  S1»0.)  HN-NUP 

BO  4?4  IF*I,NF 
RF*RKF<IF) 

IFfN.GT.KS2)  GO  TO  330 
I2*KF<IF>  I  13*12-11  I  11*12 
00  30?  1*1.13 

NN*HN-1  I  S1*S14UR(NH> XA.IAB.IOiIU) 

30?  CONTINUE 
S1*S14RF 

FR(XA.IAD.IB.IUiXF)*Sl 
SO  TO  340 

330  S1*FR(IA.IAB.IB.1U.1F)4UR(NU.IA.IA0.IB.IU)*RF 
FRCXA. XAD.XOf IU.IF)*81 
340  01*81-01 

EA<IA.IA0.10.XU,IF)*EA(IA.IA0,I0.IU.IF)4S1 
E2<IA,IAB.ID.IU.IF)*E2<IA.IAD. 10,iy,IF14SlSSl 
82*EHX( XA. XAB* IB* IV, IF)  )  S2-ABS<S2) 

XFfABOf SI ) .0T.S2)  ENK(IA.IAB.IB.IH.IF)*S1 
4?4  CONTINUE 
4?3  CONTINUE 
4?0  CONTINUE 
4?7  CONTINUE 
4?B  CONTINUE 
4??  CONTINUE 

SF-SF/NP  • 

UR1TE(4.4001>  FE.01.SP.SIF.N8 

4001  FORNAT ( '  THE  ABOVE  SCENES  HERE  GENERATES  UITN  '.FS.3. 

>'  PROBABILITY  OF  EB8E  '/'  RELATIVE  SHIFT  IS  ’.F3.3.'  PIXELS'/ 
>'  HE AN  INTENStrr  IS  '.FS.3./'  STB  BEV  IS  '.FS.3>/'  l./SNR  1$  ' 
>.FS*3) 

RNE«1./NE 
BO  3??  10*1 .NO 
00  3??  IA*1.NAL 
00  3??  IAB-l.NAB 

HRITEI4.4002)  K0<I0)>IA.IAD.(XF<IF)>IF*1iNF) 

4002  FORHATf///'  PIPELINE  CORRELATION  RESULTS!  '/'  DIOITIZINO  TO  ' 
>.I2. '  BITS'/'  AL80RITHN  ',11.'  TYPE  '.II//'  RUNNING  AVE  '. 

>' LENGTH  ',!(!!. 2?X>, 11/'  HINDOO') 

IKJHIN-l 
00  3?f  IN-l.NHI 
11*1101 

00  SOS  IF'l.NF 
SKCAdA.IAO.IB.IV.IF) 

SBF  f IF)*SORT  f (E2( IA, IAB, IB. IU, IF)-S14S1FRNC)BRNE) 

3?t  CONTINUE 

URITE(4.4003)  I1,I1,{EA(IA,IA0.I0,IH,IF),E2(IA,IA0,I0,IU.IF), 
>(NX(IA.IA0.I0.IU.IF).SBF<IF),IF*1.NF),KE(IA,IA0,I0,IU> 

4003  F0RHAT(I3.'X'.I2,4f'  I',4F7.3).'  !'.I4) 

3??  CONTINUE 

00  TO  t 
???•  STOP 
ENO 


13.21.3?  >printf  Hull  fortran 
C  B. DEFOE  10/22/80 

C  ROUTINE  TO  FORM  GAUSSIAN  DISTRIBUTION  OF  RANDON  VARIABLES 

FUNCTION  GAUSS(ISEED.AiS) 

GAUSS‘0. 

DO  1  1*1.12 

8AUSS*SAUSS4RAN( ISEEO) 

1  CONTINUE 

8AUBS*<6AUSS-4.)A8tA 

RETURN 

ENO 


13.21.58  >printf  iaaaa  fortran 

C  0.  DEFOE  9/23/80  PR06RAH  TO  PRODUCE  IHA8E  ON  TELETYPE  FRON  ARRAY  P 
C  VALUES  OF  P  BETWEEN  0.  AND  l.t  'IS'  IS  SCALE  FACTOR 

C  INAGE  SEGNENTED  TO  FIT  ON  132  CHAR.  LINE  IF 

C  NECESSARY 

SUBROUTINE  IHAGEfNX.Pl >P2) 

DIMENSION  P1(10),P2(10) 

L0GICALB1  NCR,NFF,NLF.N<4.24).NN( 4,4) >NL1( 10.4) >NL2t 10.4) 

BATA  NCR/Z1S/.NLF/Z0B/.HFF/Z0C/ 

BATA  N/ZC8 .ZAB.ZBD.Z7B. ZC9 , ZAO , 

>  ZC8.ZAB.ZB0.Z7B.ZC9.Z40. 

>  ZC8. ZAD.ZBB.Z7B.Z7A. ZAO, 

>  ZCS.ZAD.ZBD.Z4D.Z5D.Z40. 

>  ZC8.ZD4iZF0.Z40.Z5D.Z5Cr 

>  ZC1.ZD4.ZE5.Z7A.Z7B.Z40. 

>  ZCl>ZD4rZES,Z7E,Z4E,Z40. 

>  ZC1.ZD4.ZE5.Z7E.Z40.Z7A, 

>  ZC1.Z04.ZE5.Z7E.Z40.Z40. 

>  ZC8.ZF0.Z7A.Z40.Z40.Z40. 

>  ZC8,ZFO,Z4E,Z40,Z40,Z40, 

>  ZC8.ZD4.Z7A.Z40.Z40.Z40. 

>  ZC8.Z04.Z4E.Z40.Z40.Z40, 

>  ZB4.Z40.Z7A.Z40.Z40.Z40, 

>  ZD4,Z4E,Z40,Z40.Z40.Z40, 

>  ZB4,Z7E,Z40,Z40,Z40,Z40, 

>  Z4C.Z4E.Z7E.Z40.Z40.Z40. 

>  Z4C.Z4E.Z40.Z40.Z40.Z40. 

>  ZAC, Z4E. ZAO, ZAO, ZAO. ZAO. 

>  Z60.Z7A.Z40.Z40.Z40.Z40/ 

BATA  NN/Z4E,Z40.Z40,Z40,Z40,Z40, 

>  Z7E, ZAO, ZAO, ZAO, ZAO, ZAO. 

>  Z4A, ZAO, ZAO. ZAO, ZAO. ZAO, 

>  ZOO,  ZAO,  ZAO,  ZAO,  ZAO,  ZAO. 

>  Z70, ZAO, ZAO, ZAO, ZAO, ZAO, 

>  ZAO,  ZAO,  ZAO,  ZAO.  ZAO,  ZAO/ 

EQUIVALENCE  (N(1,21),NN(1,1>) 

CALL  CSSUERR.'TERHDEF  ','B  ','XL  ','ON  ') 

00  10  IX*1,10 

L 1*23 .API (IX >41. 3 I  L2*23*tP2f IX)41«S 
DO  10  LX*1,4 

NL1(IX.LX)*N(LX.L1)I  NL2(IX,LX)«N(LX,L2> 

10  CONTINUE 
DO  12  LX*1,4 

12  URITEf4,4000)(NLlfIX,LX),lX*l,10),(NL2(IX,LX),IX*l,10) 

4000  F0RHAT('4',I0A1,IX,10A1) 

UNITE  (4,4000)  NCR 
URITE(4,4000)  NLF 

CALL  CSSUERR.'TERNOCF  ','B  '.'XL  '.’OFF  ') 

RETURN 

ran 


STAGE  II 


V -.L-  mo-,  ir*T** 


C  PROGRAM 1  PIPELINE  CORRELATION 

C  CAI  DIVISOR  RECON  OPTICAL 

C  D.  DEFOE 

C  DATE!  10  FEDRUART  1982 

C 

INTEGER  SEED 
REAL  NS 

DIMENSION  Pl(7>9).  P2I7»9> ,  ILN(0> 1 >2.3.3) » 

>IASNt 1,2,3,3) *  IASD(1,2,3,3),  NR<8.1 .2i3>3) ,  L(l)> 
>EAtl,2,3,3,l) ,  E2I1,2,3,3,D,  EHXI1,2,3,3,D,  K0I1.2), 

>KD(3>.  KFIl),  SFI1),  ILD<8.1>2>3>3)>  FRI1,2,3.3,1), 
>SDFID,PLlt9),PL2I9)«RKFID,KEIl,2,3,3) 

DATA  KD/A,4,2/  KF/8/ 

NN-71  NL-I001  NH«8I  NAL-1)  NAD-21  ND-31  NF*1 
JUIH-AI  JUIX-8 
CH-.25I  U-.02I  CNH«A. 

NA-JUIXI  NN-JWIXtl)  NUI-JUIX-JUIN+l 
RM-SORT(.S) 

KS1-NM4NA-1I  KS2-KSDNU-1 

NHN-NH-II  NMP-NM+U  NLH-NL4KS2 

NNM«NN-1 1  NNP-NNM)  NAP-NAM 

HUH-NV-11  NUP-NUH1  NIIIN»NNI-1I  NUIP-NUIM 

JUINN-JNIN-1 

NNH-CHH4U1  UN-UNH4UNNI  UNP-l.-UN 
DO  I  IF-1,NF 

1  RKF(IF>«t./KF(IF> 

2  SP-SP2-0.I  NP-NE-0 
DO  8  ItIMfNUI 

DO  7  ID-l.ND 
DO  A  1  AD-1,  NAD 
DO  S  IA-l.NAL 
KE(IA,IAD,lD,IN)-0 
DO  4  IF-I.NF 

EAIIA,IAD,ID,IN.IF)»E2llA.IAD,ID,IN.IF)«£NXIIA.IAD.I0.1U.IF)sO. 

4  CONTINUE 

5  CONTINUE 
A  CONTINUE 

7  CONTINUE 

8  CONTINUE 

READ(5,S000)  PE, Dl, NS, SEED 
5000  FORMAT! 3F10. 4, 110) 

IF(PE.EQ.O)  GO  TO  9999 

D2-1-D1  )  C*I.3-ADS1.3-DD)*CM  I  MP-NM 

MPM-MP-1 

MA-NA  I  HU-NH 

DO  499  HM.NLN 

HPM-HP  »  MP-HP+1 
IF(NP.EO.NMP)  HP-1 
IFIH.LT.NH)  GO  TO  10 
NA-NAA1 

IF I HA. ED .NAP)  HA-1 
IF(H.LT.KSl)  60  TO  10 

MU-HU41  . 

IFIHM.ED.NUP)  MU-1 

NE-NE41 

10  DO  199  NR-l.NN 

N-NNP-NR  I  U-RAN(SEED) 

IF(U.GE.PE)  80  TO  20 

IU-S 

DO  TO  30 

20  U-RAN(SEED)  I  1U-48U41  I  LID-0 

30  IFIIU.GE.3)  LID-1 

LI2)-LI3)-LI4)-LI1) 

IFIIU.8E.4)  80  TO  100 
U-NANISEED) 

IF1IU.E8.2)  80  TO  SO 
JU-48091 

IFIIU.NE.D  80  TO  40 
L(JV»1 
80  TO  100 
40  L(Jtt>*0 
80  TO  100 
30  JU-A8U 

IFIJU.8E.3)  80  TO  AO 
Llll-1  I  LI JU42)-t 
00  TO  100 
AO  LIJU-D-1 

IFIJU.EO.S)  JU-2 
LI JU)-1  * 

100  T1-T2-T3-T4-.S 

IFIN.LT.NN)  Tl-Pl(HP,Ntl> 

IFIN.EO.l)  80  TO  110 
T3-PKHPH,N) 

1FIN.LT.NN)  T2*P1IHPN,N41 ) 

1FIN.8T.D  T4-P1IHPN,N*D 

110  NA-IT19LID4T38LI3)4IT29LI2)4T4tLI4>)8RH)/(LID4LI3)4 

>ILI2)HI4))tRM> 

IFIIU.E8.S)  80  TO  120 

0*1000  I  O-OAOSOISEED.NA.O)  )  00  TO  ISO 


120  U-RANISEED) 

IFIUA.LT .UN)  60  TO  130 
1FIUA.GT . UNP)  80  TO  140 
Q-NNP9U 

IFIU.GT.UA)  O-O+UN 
60  TO  ISO 

130  0- 1 VA4UNH ) > 1 1 • -U  >  4U 
GO  TO  ISO 
140  0-(UA-UNH)*U 
150  PI  IMP, N)-Q  I  SP-SP40 
8P2-SP24080  I  NP-NP+1 
IFIN.EO.l)  GO  TO  199 
T1-P1IMPH*N) 

T3-P1(HP,N> 

T2-T1-T3 

0-T1-D18T2 

P2IMPM,N)-Q  I  SP-SP+Q 
SP2-SP24QIQ  •  NP-NP+1 

199  CONTINUE 

AOOO  F0RMATI2I IX, 10IF5. 3, IX) ) ) 

200  IFIH.LT.NH)  60  TO  499 
8IP-S0RTI!SP2-SP*SP/NP)/NP)  1  SIN-NStSIP 
IFtN.OT *NH)  60  TO  210 

DO  209  HN-1,NHH 
DO  208  N-DNNM 

T1  -PI  I  MM ,  N )  +GAUSS I  SEED ,  0 .  ,  SIN  > 

T2-P2 1  HN,N  >  +6AUSS I  SEED  ,0 . ,  SIN  > 

T1-AMAX1 IAHIN1 IT1 , 1 , ) ,0. ) 

T2-AMAX1 I AMIN1 I T2, 1 . ) ,0. ) 

P1IHN,N)-PL1IN)-T1 

P2INH,N)-PL2IN)-T2 

208  CONTINUE 

CALL  IHAGEI8,PL1,PL2) 

209  CONTINUE 
80  TO  220 

210  DO  219  N-1,NNN 
TI-P1(HPH,N)46AUSS(SEED,0. ,SIN) 
T2-P2(HPH,NHGAUSSISEED,0.  ,SIN) 

T1 -AMAX1 I ANINl I T1 ,1 . ) ,0. ) 
T2-AHAXHAMIN1IT2, 1 .  ),0.  ) 

PI IHPH,N)*PL1IN)-T1 
P2(MPM,N)-PL2IN)-T2 
219  CONTINUE 

CALL  INAGEI8,PL1,PL2> 

220  DO  498  IA-l.NAL 


DO  497  I AD-1, NAD 


DO  49A  ID-1, ND 
KSN-KDIID)  I  KS-2MKSH 
IFIM.LE.KSD  GO  TO  230 
HAN-HA4JUINH 
IFIHAH.6T.NA)  MAN-HAN-NA 
DO  229  I»>1,NVI 
NAN-NAN- 1 

I  IFINAH.LE.O)  HAH-HAH+NA 

IA8NIIA,IAD,ID,IH)-IASNI1A,IAD,1D,IU)-ILNIHAM,IA,IAD,1D,IU) 
IADD(IA,IAD,ID,IH>-IASD(IA,IAD,ID,IU)-tLD(HAM,IA,IAD,ID,IN) 
229  CONTINUE 
230  ISI-I82-0 

DO  249  N-DJNIN 


IS1-ISDI4 


IS2-IS24TT 

249  CONTINUE 

IUIHA,IA,IAD,ID,D>IS1 
ILDINA,IA,IAD,ID,1)-I82 
1FINNI.E0.D  60  TO  260 
DO  239  INN-1, RUIN 
N-JUIN4INN  I  IU-1NM41 


I81-IS14I4 


IS2-IS2II4 

1LNIHA,IA,IAD,ID,IU)-IS1 
ILDINA,  IA,IAD,ID,IU)-1S2 


23?  CONTINUE 

240  IFIH.LT.KS1)  SO  TO  404 
IFIH.6T.KS11  60  TO  280 
DO  270  IWN-l.NWI 
IS1>IS2«0 
IW-NNIP-IWN 
DO  240  NAH-IWN.JWIX 
ISl-IS14ILN(HAHr IA. IAD. ID. IU) 

IS2*IS2+ILD<HAH.IA.IAD.ID>IW> 

240  CONTINUE 

lASNIIA.IAO.IO.IWl-lSl 
IASD<IA»IAD.ID.IW>*IS2 
270  CONTINUE 
60  TO  300 

260  DO  280  IU-1.NWI 

lASNIIArlAD.IO.IWl-IASNdAf  1AD.  ID.1WHILNIHA.IA.IAD.ID.IU) 
IASD( IAi IADi I0.IW)-IASDdAfIADiIDfIU)4ILD(MA.IA.IAD.IDfIU> 

280  CONTINUE 

200  IFIH.LE.KS21  80  TO  300 
DO  200  IU-1.NNI 
DO  200  IF-l.NF 
HH-HW-KFIlFll  RF-RKF(IF) 

IF(NN.LE.O)  NN-HH4NW 

FRIIA.IAD.ID.1W. IFl-FRdA.IAD.lD.IW. IFl-WRIRN.lA.lADf ID.IW10RF 
200  CONTINUE 

300  DO  403  IW-l.NNI 

Tl-IASNdA.lAB.ID.INI  I  T2-IASDdA.IAD.ID.INI 
T3*A»SCT2> 

IF(T3.0T..S8ABS!T1>.ANB.T3.8E.5.>  00  TO  301 
SO  TO  302 

301  KEdA.IAD.ID.INl-KEdA.lAD.ID.IWlOl 

302  WRINW.IA.1AD.ID.IW) 

IF<H.LT.KS2)  60  TO  403 
11-01  Sl-O. I  NH«NWP 

DO  404  IF-l.NF 
RF-RNFdF) 

IF(H.6T.K82)  60  TO  330 
I2«KFdF)  I  I3-I2-I1  I  11-12 
DO  300  1-1.13 

HN-NN-1  I  Sl*SltWR(NN,IA»IAD»IJ.IW) 

300  CONTINUE 
SI-SURF 

FRdA.IAD.ID.IW.IFl-Sl 
80  TO  340 

330  SI-fRIIA.ZAS.IDflUrlFlOWRINV.IA.IAO.ZD.IVlFRF 
FRdA.IAD>ID>IN>IF)-Sl 
340  81-81-81 

EAdA.IAD.IB.IN.IFl-EAIIA.lAB.ID.IW.IFlOSl 
E2(IA.IAD>ID.IW>IF)-E2dA.lADr  ID.  IN)  IFHS14S1 
82-ENXIIA. IAD.  IB.  IN. IF 1  I  S2-A8SIS2) 

IF(AD8<81).6T.82)  ENX<IA.IAD>ID»IU.IF)-S1 
404  CONTINUE 

403  CONTINUE 

404  CONTINUE 

407  CONTINUE  . 

408  CONTINUE 
400  CONTINUE 

SP-3P/NF 

NRITE(4>4001)  PE.D1.SP.SIP.NS 

4001  FORNATI'  THE  ABOVE  SCENES  WERE  GENERATED  WITH  '.F5.3. 

>'  PROBABILITY  OF  EDGE  '/'  RELATIVE  SHIFT  IS  '.F5.3.'  PIXELS'/ 
>'  HERN  INTENSITY  18  '.FS.3./'  STD  DEV  18  'rFS.3»/'  l./SNR  IS  ' 
>»F3.3) 

RNC-l./NE  , 

BO  300  ID«1>NB 
BO  300  IA-1.NAL 
BO  300  IAB-1.NAB 

WRITE (4.40021  KDdB)  >IA>  IAD.  (KFI1F)  >  IF-1  .NF) 

4002  FORNATI///'  PIPELINE  CORRELATION  RESULTS!  '/'  DIOITIZINS  TO  ' 
>.I2. '  BITS'/'  ALOORITHH  '.II, '  TYPE  '>11//'  RUNNIH8  AVE  '> 
>'LENOTH  '.II/'  WINDOW') 

Il-JWIN-1 
BO  300  1W-1.NWI 
11*1101 

BO  300  IF-l.NF 
11-CAdA.IAB.IB.IW.IF) 

BBFIIFl-SORTI (E2IIA.1AB. tB>IW)IF>-S18818RNC)8RNC) 

SOS  CONTINUE 

WSITEI4.4003)  Il.Il.IEAdA.IAB.ID.IW,IF).£2dA.IAD.ID.IW,IF)» 
>ENX(IA.IAB.IB.IW.IF),8BF(IF),IF-1.NF),RE(!A.IA8.ID.IW) 

4003  FSftHATIIl.'X'. 12.11'  !'.4F7,3).'  !'.I4) 

300  CONTINUE 

80  TO  2 
0000  STOP 
END 


15.30.37  > 


STAGE  III 


c  program:  pipeline  correlation 

c  CAI  D1VISON  RECON  OPTICAL 

C  D.  DEFOE 

C  DATE:  10  FEBRUARY  1P82 

C 

INTEGER  SEED 
REAL  NS 

DIMENSION  PK  4 *  9  > »  P2< 4.9) .  ILNI8.3.3). 

>1 ASN< 3 v 3 > •  IASD(3* 3f2) .  UR(8<3.3>.  L<4>. 
>EA ( 3  >  3  >  »  E2(3.3>.  EMXI3.3). 

>KD  <  3 ) t  ILDCS.3.3.2) .  FR(I.I). 

>KE(3>3) 

DATA  KD/4.4.2/ 

NM-4I  NL-1001  NU-81  ND-3 
JUIN-41  JUIX-8 
CM-.2SI  U-.04I  CBH-4. 

NA- JBIX  •  NN-JUIX41)  NUI-JUIX-JWIN41 
RH-SORTI .5) 

KS1-NN4NA-1 I  KS2-KS14NU-1 

NMM-NM-1)  NMP-NM41)  NLM-NL4KS2-1 

NNM-NN-ll  NNP-NN4K  NAP-NA41 

NUM-NU-ll  NUP-NU41I  NBIH-NMI-ll  NUIP-NUI41 

JMlNM-JBIN-1 

BNH-CNHSUt  UN-UNH4UNHI  UNP-l.-WN 
KF-81  RF-l./KF 
1  SP-SP2-0.I  NP-NE-0 
DO  9  IU-1.NBI 
DO  8  ID-l.ND 
KECID.IM>-0 

EA(ID.IU)-E2<ID.IU)-EMX( ID. IW>-0. 

8  CONTINUE 

9  CONTINUE 

READC3.3000)  PE. D1 .NS. SEED 
SOOO  FORMAT I 3F 10*4.110) 

IF(PE.EB.O)  GO  TO  9999 
D2-1-D1  »  C-I . 5-ABS( .5-D1 ) )*CM 
HP-NHf  HPH-HP-1) 

HA-NA  »  HU-NM 
DO  499  H-l.NLM 


HP-HP41  MPf'sMP 

IF  (  HP  •  EQ  » NMP  >  MP-1 
IF(N.LT.NN)  BO  TO  10 
HA-MA41 

IF IMA .EQ . NAP)  MA-1 
IF(M.LT.KSl)  80  TO  10 
HM-HM41 

IFIHB.EO.NUP)  MU-1 
IF<M.LT.KS2)  GO  TO  10 
NE-NE41 

10  DO  199  HR-l.NN 

N-NNP-NR  »  U-RANISEED) 

IFIU.OE.PE>  GO  TO  20 

1U-S  • 

GO  TO  30 

20  U-RAN(SEED)  I  IU-4SU41  I  L(l)-0 

30  IF<IU.GE.3>  L( 1 >■! 

L<2>-L<3)-L(4)«L<1> 

IF < IU.QE. 4)  BO  TO  100 
U-RAN(SEED) 

IF( IU.E0.2)  GO  TO  SO 
JU-4SU41 

.  IF(lU.NE.l)  GO  TO  40 
L1JU>-1 
BO  TO  100 
40  L<J U>-0 

80  TO  100 
SO  JU-48U 

IFUU.BE. 3)  BO  TO  40 
LID-1  I  LCJU42>-1 
BO  TO  100 
40  LIJU-D-1 

IFIJU.EB.S)  JU-2 
LIJU>-1 

100  T1-T2-T3-T4-.S 

IFCN.LT.NN)  Tl-PKMP.N+l) 

1FIM.EB.1)  BO  TO  110 
T3-P1 IHPH.N) 

IF(N.LT.NM)  M-P1CHPN.M41) 

If  1N.BT.D  T4-P1 CMPN.N-1 > 

110  NA— IT18L1 1 >4T3SL13)4IT2SL(2)4T48L14> )SRH)/1L1 1 >4L (3)4 
>(L(2>4LI4>>MH> 

IF  1 1U.CB.S)  80  TO  120 

0-IB8M  I  B- GAUSS  (SEED.  BA.  0 )  I  BO  TO  ISO 
120  S-S4SI  StCB  > 

IFIMA.LT. IM>  BO  TO  130 
IFIBA.BT.NNP)  SO  TO  140 

AaMMARU 

IF1B.BT.MA>  8-84 BN 
BB  TB  ISO 


130  0-(UA4MNH)*(l.-U>4U 

GO  TO  ISO 
140  Q-IBA-BNH) SU 
ISO  Pl(MP.N)«a  >  SP-SP40 
SP2-SP24Q4Q  >  NP-NP41 
IF(H.EO.l)  GO  TO  199 
Tl-Pl(MPM.N) 

T3-PKMP.N) 

T2-T1-T3 

0-T1-D1ST2 

P2(HPH»N)-Q  I  SP-SP4G 
8P2-SP24QS0  I  NP-NP41 

199  CONTINUE 

4000  F0RMAT(2(1X.10(F5.3.1X)>) 

200  IFIN.LT.NM)  GO  TO  499 

8IP-SQRT ( <SP2-SPSSP/NP)/NP>  »  SIN-NS4SIP 
IF(M.GT.NM)  GO  TO  210 
DO  209  MH-l.NMM 
DO  208  N-l.NNM 

T1-P1(MH.N>4GAUS8(SEED.0« .SIN) 
T2-P2<MH.N)4GAUSS(SEED.O.  .SIN) 

T 1 - AMAX1 (AMINl(Tl.l.).O.) 

T2-AMAXI (AMIN1 (T2, 1 . ) .0. ) 

P1(MM.N)-T1 

P2(HM.N)-T2 

208  CONTINUE 

209  CONTINUE 
00  TO  220 

210  DO  219  N-l.NNM 

Tl-Pl < HPH.N>4BAUSS( SEED .0. .SIN) 
T2-P2(MPM.N)4GAUSS(SEED.O. .SIN) 

T1-AMAX1  (AMINKT1.1,  >,0.  ) 

T2-AMAX1 ( AHIN1 (T2. 1 . > .0. ) 

PI (MPM.N)-Tl 
P2(HPMrN)*T2 
219  CONTINUE 

220  DO  494  ID-1 .ND 
KS-4SSKDIID) 

IF (M.LE.KS1 >  GO  TO  230 
NAM-HA4NBI 

IF (MAN.GT . NA)  HAM-HAN-NA 
DO  229  IB-l.NUI 
NAN-NAN- 1 

IFIHAM.LE.O)  HAM-HAN4NA 

IASNI  ID.  IB  )  -  IASNI  ZD.  IU)-ILN(NANr  ID.  IB) 

DO  229  JD-1.2 

IA8D(ID.IU.JD)-IASDIID. IB. JD)-ILD(NAH. ID.IM.JD) 
229  CONTINUE 

230  I81-IS21-IS22-0 
DO  249  N-l.JBIN 


'  IS1-IS141S 


IS21-IS214IS 


IS22-IS224IS 

249  CONTINUE 

ILNIHA.ID.D-ISl 

ILDCNA.ID.lt  D-IS211  ILDIHA.  ID.l  >2)-IS22 
lIFINMI.EO.l)  GO  TO  240 
DO  2S9  IBH-l.NMIN 
N-JBIN4IMH  >  IB-IBN41 


I81-IS14I8 


IS21-IS214IS 


IS22-IS224IS 

ILN(NA.IS.IB)-I81 

ILD(NA.ID.tB.l)-IS21>  ILDIHA. ID. IB. 2>-IB22 
2S9  CONTINUE 

240  IF(M.LT.KSl)  80  TO  494 
IF1M.0T.K81)  60  TO  280 
DO  279  IBN-l.NMI 
181-1 82 1-1822-0 
IB-NB1P-IBN 
DO  249  HAM-IMK. JMZX 
I81-I814ILN1MAH. ID. IB) 

IS21 -1S2141LDCNAM. ID. IB.l ) )  IS23-IS224ILDIHAH.1D.IM. 
249  CONTINUE 

IASNI ID. IB) >181 

lASDIID.IH. 1 >-1821 >  IASD1 ID. IB. 2 >-1822 
279  CONTINUE 
00  TO  300 


280  DO  289  IW1.NUI 

:asncid,iw>-iasn(id.iw>+iln(ma>id»iw> 

DO  289  JD*1 *2 

IASDUD.IU, JD>*IASD<ID.IU. JD > +ILD < MA . ID , I U . JD) 

289  CONTINUE 

290  IF(M*LE.KS2)  60  TO  300 
DO  299  IU*1.NWI 
NM*NH-KF 

IF(HN.LE.O)  MM*HH+NM 
FR<ID.IW)*FR<ID.iy)-tlR<HM.IB.IH>«RF 
299  CONTINUE 

300  DO  495  IH-l.NHI 
T1>IASN(ID>IU) 

IFITl.OE.O.)  00  TO  301 
T2«IASDC1D.IU.1)I  00  TO  302 

301  T2*IASD<ID.IH.2> 

302  T3»ABS(T2) 

IF<T3.GT..5*ABS<T1),AND.T3.GE.3.>  00  TO  303 
00  TO  304 

303  KE<ID.iy>-KE<ID.IU)+l 

304  UR(MU.ID.IU) 

IF<M.LT.KS2>  00  TO  495 
S1*0.I  HH*NUP 
IF<N.GT.KS2>  00  TO  330 
DO  309  1*1. KF 

MM*MH-1  I  S1*S1+UR<MM.ID.IU> 

309  CONTINUE 
S1*S1BRF 
FR<ID.IU)*S1 
00  TO  340 

330  S1*FR<  ID.  IU) FUR  <MU.  ID.  IU  >  BRF 

FR< ID. IN ) *01 
340  S1-S1-D1 

EA  <  ID.  Ill)  *EA  (ID  .110481 
E2( ID. IU)*E2( ID. IU) 481*81 
S2*EMX<1D.IU>  I  S2*ABS(S2) 

IF(ABS(S1> .GT.82)  EMX< ID. IU)*S1 
495  CONTINUE 
490  CONTINUE 
499  CONTINUE 
8P-8P/NP 

URITE(6* 6001 )  PE.D1.SP.SIP.NS 

4001  FORMAT ( '  THE  SCENES  MERE  GENERATED  WITH  '.FS.3. 

>'  PROBABILITY  OF  EDOE  '/'  RELATIVE  SHIFT  IS  '.F5.3.'  PIXELS'/ 
>'  MEAN  INTENSITY  18  '.F5.3./'  8TD  DEV  IS  '.F5.3./'  l./SNR  IS  ' 
>.F5.3) 

RNE*1./NE 
DO  599  IDM.ND 
NRITE(4.4002>  KD<ID).NE.KF 

4002  FORMAT!///'  PIPELINE  CORRELATION  RESULTS!  '/'  DIGITIZING  TO  ' 
>*I2» '  BIT8'/'  FOR'. 15.'  LINES'//'  RUNNING  AVE  '. 

>'L£NQTH  '.II/'  UINDOU') 

Il-JHIN-1 
DO  599  IW1.NUI 
11*1141 
81*EA< ID. IN) 

S2*E2( ID. IN) 

SDF-SORT ( ( S2-S1 tSIFRNE > BRNE > 

S1-S1BRNE*  S2*SQRT(S2*RNE> 

WRITE  <  4  >  4003  >  II. II. 81. 82. 

>EMX (ID.IU).SDFfKE(ID.IU) 

6003  FORMAT <I3.'X'.I2.1('  t'»4F7.3>,'  I', 14) 

599  CONTINUE 
00  TO  1 
9999  STOP 
END 


17.32.06  > 
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lofftrrintf  pip«c  fortran 

PROGRAM:  PIPELINE  CORRELATION 
CAI  DIVISON  RECON  OPTICAL 
D.  DEFOE 

DATE!  10  FEBRUARY  1982 


STAGE  IV 


INTEGER  SEED 
REAL  NS 

DIMENSION  PI <  6 . 10 ) r  PZI4.10).  ILNI9).  ILDI9.2) 

>UR  <  8 ) »  L ( 4  > 

KD-41  KS-256 
NM-61  NL-100)  NU-8 
JUIN-JUIX-9 
CM-. 251  W-.04I  CUH»4. 

NA-JUIXI  NN-JUIX+1)  NUI-JUIX- JUIN+ 1 
RH-80RT ( .5) 

KS1-NN+NA-1 »  KS2-KSI  +NW- 1 
NMM-NM-ll  NMP-NM+ 1  *  NLM-NL+KS2-1 
MNM-NN-ll  NNP-NN+1)  NAP-NA+1 
HUM-NU-ll  NUP-NU+ 1 I  NUIH-NUI-l)  NUIP-NUI+1 
JUINM-JUIN-1 

UNH-CNH4UI  UN-UNH+ UNHI  UNP-l.-WN 
KF-81  RF-l./KF 
1  SP-SP2-0. )  NP-NE-0 
KE-OI  EA-E2-EMX-0. 

READIStSOOO)  PE.D1 iNStSEED 
5000  FORMAT I 3F 10.4.110) 

IF (PE.EQ.O)  GO  T"  9999 
82-1-01  I  C-< .5-  >BS< .5-01) )*CM 
NP-NMt  MPM-MP-1 
NA-NA  J  NU-NU 

DO  499  M-l.NLM  > 

MPM-  MP 

HP-MP+1 

IF(MP.EQ.NMP)  HP-1 
IFIN.LT.NN)  GO  TO  10 
MA-MA+1 

IF ( MA.EQ.NAP )  MA-1 
IF(M.LT.KSl)  GO  TO  10 
NU-HM+1 

IF ( MU . EQ  » NUP )  MU-1 
IFIH.LT . KS2 )  GO  TO  10 
NE-NE+1 

10  DO  199  NR-l.NN 

N-NNP-NR  )  U-RAN ( SEED ) 

IF(U.GE.PE)  GO  TO  20 
IU-5 

GO  TO  30 

20  U-RAN(SEED>  *  IU-4FU+1  I  L(l)-0 
30  IF( IU. GE • 3 )  LI  1 ) - 1 

L<2)-L<3)«L(4)-L(1) 

IF  I  IU.GE.  4  >  GO  TO  100 
U-R ANISEED) 

IF < IU. EG. 2 )  GO  TO  50 
JU-44U+1 

IF(lU.NE.l)  80  SO  40 
L<JU»-1 
80  TO  100 
40  L  <  JU )  -0 

00  TO  100 
SO  JU-44U 

IF( JU.GE.3)  GO  TO  40 
L<1)-1  »  L <  JU42  > « 1 

GO  TO  100  . 

40  LIJU-D-l 

XFI JU.EQ.5)  JO-2 
LIJU>-1 

100  T1-T2-T3-T4-.5 

IFIN.LT.NN)  Tl-PUHP.N+l) 

IFIM. EQ.l)  80  TO  110 
T3-P1INPH.N) 

IFIN. LT.NN)  T2-P1INPH.N41) 

IFIN.OT.D  T4-P1INPN.N-1) 

110  UA-IT1«H1)+T3*H3)4IT2»L(2)4T4*H4))*RH>/ILID+LI3>4 

>ILI2>*L(4)>*RH) 

IF  I IU.EO.S)  GO  TO  120 
0-IU9U  )  0-8AUS8 ISEED.UA.Q)  1  GO  TO  150 
120  U-RAN18EED) 

IFIUA.LT. UN)  GO  TO  130 
IF1HA.BT.UNP)  80  TO  140 
Q-UNPtU 

IFIU.GT «UA)  0-Q4UN 
00  TO  ISO 

130  0-IUA90NH)*I1.-U)+U 

00  TO  ISO 
140  0-IUA-UNH>*U 
150  PlIHP'Nl-0  •  8P-SP40 
BP2-SP24090  •  NP-NP41 
IF IN.EQ. 1 >  00  TO  199 
T1-P1IMPM.N) 

T3-P1IMP.N) 

T2-T1-T3 

0-T1-D19T2 

P2IMPM.N).0  •  8P-SP+0 
092-SP24090  I  NP-NP91 
199  CONTINUE 

4000  F0MMT(2(lX*10(FS.3tlX)l)  ..  ^ 


200  IFIM.LT.NM)  00  TO  409 

SIP*SQRT! !  SP2-SP4SP/NP) /NP>  I  SIN*NS»SIP 
IF(N.GT.NH)  00  TO  210 
DO  209  Hh* 1 . NHM 
DO  208  N«I.NNH 

T 1 »P1 ( HM .N)4GAUSS ( SEED • 0  * .SIN) 

T2*P2! HH. N ) +GAUSS ( SEED » 0 . .SIN) 

T1-AHAX1 < AHIN1 < T1 » 1 . ) .0. ) 

T2*AHAX1 ( AMI N 1 <T2.1.).0»> 

PI !HH»N)*T1 
P2!HH.N)«T2 

208  CONTINUE 

209  CONTINUE 
GO  TO  220 

210  DO  219  N-l.NNH 

T1 *P1 (HPH.N>+ GAUSS ( SEED. 0. .SIN) 
T2*P2!MPH»N)+GAUSSISEED»0. .SIN) 

T1-AHAX1 !AH1N1!T1»1«).0* ) 

T2*AHAX1  (ANINKT2.1  •  )  .0.  ) 

P1<NPM.N)»T1 
P2( HPH  »N  >  *T2 
219  CONTINUE 

220  IF(H.LE.KSl)  GO  TO  230 
IASN*IASN-ILN!HA) 

IASD1-IASD1-ILD!HA.1>I  IASD2*IASD2-ILD!HA.2> 
230  IS1*IS21*IS22»0 
DO  249  N'l.JUIN 


IS1-IS1+IS 


IS21=IS21+IS 


IS22-IS22+IS 

'249  CONTINUE 

ILN ( NA  >  * IS1 

ILD!HA.1)*IS21I  ILD!HA»2)-IS22 
260  IF(H.LT.KSl)  GO  TO  499 
IF (M.GT .KS1 )  GO  TO  280 
IASN*IASD1*IASD2*0 
DO  269  HAH“1  •  JUII X 
IASN*IASN4ILN!HAH) 

IASBl*IASD14ILD!HAH.l)f  IASD2-IASD2+ILDIMAH.2) 

269  CONTINUE 
GO  TO  300 

280  IASN*IASN4ILN!HA> 

IASD1-IASD14ILDIHA.1))  I ASD2*I ASD2+ILD ! HA » 2 > 

290  IF!H.LE.KS2>  GO  TO  300 
HH-HW-KF 

IF(HM.LE.O)  HH-MM+NH 
FR*FR-WR!HM)»RF 

300  Tl-IASN 
IFITl.OE.O.)  GO  TO  301 
T2-IASD1!  GO  TO  302 

301  T2«IASD2 

302  T3-ABS<T2>  . 

IF!T3.GT..5*ABS!T1).AND.T3.GE.5.>  00  TO  303 

00  TO  304 

303  KE*KE41 

304  UR(HH) 

IF ( H • LT * KS2 )  00  TO  499 
HH-NMP 

IF!M.8T *KS2)  GO  TO  330 
FR-O. 

DO  309  1*1 . KF  , 

MH-HH-1  I  FR*FR4WR!MN> 

309  CONTINUE 
FR-FRBRF 
00  TO  340 

330  FR*FR4WR!MW)*RF 
340  S1*FR-D1 
EA-EA981 
E2-E24S1SS1 
S2-ABS(EHX) 

IF!ABB!S1).GT.S2)  EHX-S1 
499  CONTINUE 
SP-SP/NP 

UR1TEI6.6001 )  PE.D1 .SP.SIP.NS 

6001  FORMAT! '  THE  SCENES  WERE  GENERATED  WITH  '.F5.3. 

>'  PROBABILITY  OF  EDGE  '/'  RELATIVE  SHIFT  IS  '.FS.J.*  PIXELS'/ 
>'  MEAN  INTENSITY  IS  '.FS.3./'  STB  DEV  IS  '.F5.3./'  l./SNR  IS  ' 
>.FS.3> 

RNE*1./NE 

WRITE! 6.6002  >  KD.NE.KF 

6002  FORMAT!///'  PIPELINE  CORRELATION  RESULTSt  '/'  DIGITIZING  TO  ' 
>.I2. '  BITS'/'  FOR ' » IS# '  LINES'//'  RUNNING  AVE  '• 

> 'LENGTH  '.II/'  WINDOW') 

I1*JUIN-1 

11*1141 

SDF*S0RTI!S2-S1»S1»RNE)*RN£) 

EA-EAMNE)  £2*SGRT!E2SRNE> 

WRITE! 6. 6003)  1 1 . 1 1 .EA.E2.ENX. SDF.KE 

6003  FORMAT !I3.'X'.I2»1!'  t'.4F9.3>.’  »'.I4) 

SO  TO  1 

9*9*  STOP 
END 


STAGE  V 


16*20*24  >*rintf  riree  Fortr*n 

C  PROGRAM!  PIPELINE  CORRELATION 

C  CAI  DIVISON  RECON  OPTICAL 

C  D.  DEFOE 

C  DATE!  10  FEBRUARY  19B2 

C 

INTEGER  SEED 
REAL  NS 

DIMENSION  Pl<6»10>*  P2( 6  f 10 ) t  ILN(9>>  ILD(9»2>.  ILDZ(9)> 
>WR<8) r  URZ(S)  L ( 4 ) 

KD-4!  KS-256 

NM«6i  NL-1001  NU«8 

JUIN-JUIX-9 

CM*. 231  U«.04»  CWH-4. 

NA-JUIXf  NN» JUIX+ 1 !  NMI-JUIX-JUIN+1 
RH*SQRT(.5)I  R2*SQRT(2. ) 

K81-NM+NA-1I  KS2*KS14NU-l 

NMM-NM-lt  NMP-NM+ll  NLM*NL+KS2-1 

NNM-NN-ll  NNP»NN+1I  NAP-NA+1 

NNM*NH-1»  NWP-NU+ll  NUIM-NWI-1!  NWIP-NWI+1 

JUINH-JUIN-1 

UNH-CNH4U)  WN-UNH+UNHI  UNP-l.-UN 
KF"8t  RF*1./KF 
1  SP-SP2-0.I  NP-NE-0 
KE-Ot  EA"E2"EMX"0. 

KEZ«0!  EAZ"E2Z-£MXZ-0. 

READ ( 5  r  5000  >  PE >D1 >NS • SEED 
3000  FORMAT (3F10.4illO) 

IF(PE.EQ.O)  GO  TO  9999 
D1H-D16.5 

D2-1-D1  I  C-< *S-ABS< *5-Dl > >»CM 
02H-D28.3*  D3-D14D2 
MP-NMI  MPM-MP-1 
HA-NA  I  MU-NU 
DO  499  M"1 »NLM 

MPM  «MP 

HP-MP+1 

IF(MP.EO.NMP)  MP-1 
IF(M.LT.NM)  GO  TO  10 
HA-MA41 

IF ( MA . EG • NAP  >  MA-1 
IF (M.LT . KS1 )  GO  TO  10 
HU-MU41 

IF(MU.EQ.NWP)  MU" 1 
IF <M«LT *KS2)  GO  TO  10 
NE"NE+1 

10  DO  199  NR"1 *NN  ' 

N-NNP-NR  I  U-RAN(SEED) 

IF(U.GE.PE)  00  TO  20 

IU-3 

<u  TO  30 

20  U-J<AH<SEED>  I  IU-49U+1  »  L(l)-0 

30  IF(IU.GE.3>  L ( 1 > *1 

L<2J"L<3»"L<4>*L<1) 

IF<IU*GE*4)  GO  TO  100 
U"RAN(SEED> 

IF< IU*EQ* 2 >  00  TO  30 
JU«4*U41 

IF(IU*NE*t)  BO  TO  40  1 

L(JU>"1 
60  TO  100 
40  L(JU)-0 

GO  TO  100 
30  JU"6*U 

IF  (  JU  •  GE  .  3 )  00  TO  60 
L(1 >"1  I  L< JU42)"1 
GO  TO  100 
60  L<  JU-1 ) ■ l 

IFIJU.EO.S)  JU-2 
L ( JU ) * 1 

100  T1«T2*T3"T4*.3 

IF<N.LT.NN)  Tl«Pl IMP  rN41 ) 

IF(M.EG.l)  GO  TO  110 
T3>P1(MPM>N) 

IF(N.LT.NN)  T2"P1(MPM»N41> 

IF(N.GT.l)  T4"P1<MPM»M-1> 


110  MA-<Tl*L<l)4T3*L<3>4<T2*l<2>4T4*L(4>)*RH>/(l.<l>4L<3>4 

>(L<2)4L<4>>*RH> 

IF(IU.EO.S)  GO  TO  120 

Q>IU*U  )  Q*GAUSS(SEED>WA.Q>  I  GO  TO  ISO 
120  U"RAN<  SEED) 

IFIUA.LT.UN)  00  TO  130 
IFIUA.GT.UNP)  00  TO  140 
0-NNP4U 

IF  <  U • GT . UA )  0*0+UN 
00  TO  ISO 

130  0»<NA4UNH)*<1.-U>4U 

GO  TO  ISO 
140  Q*<HA-NNH)*U 
ISO  P1<HP»N>«0  I  SP»SP+Q 
SP2*SP24Q*0  I  NP-NP41 
IF(M.EQ.l)  00  TO  149 
Tl-PHMPN.N) 

T3*P1<NP.N> 

T2-T1-T3 

Q-T1-D19T2 

P2(MPM»N>*0  I  SP-SP4Q 
SP2*SP2+0*Q  t  NP-NP41 

199  CONTINUE 

6000  FORMAT  (2(1X.10(F5»3.1X>>> 

200  IF(H.LT.NN)  00  TO  499 
8IP»SQRT(<SP2-SP*SP/NP>/NP>  I  SIN=NS*SIP 
IF(M.QT.NH)  00  TO  210 

DO  209  HH*1  f  NMM 
DO  208  N*1 iNNN 

Tl=Pl(MM>N)+OAUSS(SEEDfO. >SIN) 

T2=P2 (MM.N>4GAUSS(SEED.O. >SIN) 

T1»AMAX1(AMIN1(T1.1. )>0.) 

T2=AMAX1 <AHINl(T2f 1 . >  >0. > 

P1(MN.N)«T1 

P2<HN.N)»T2 

208  CONTINUE 

209  CONTINUE 
00  TO  220 

210  DO  219  N*liNNH 
T1=P1(MPN.N>4GAUSS(SEED.0. .SIN) 

T2=P2< MPM»N)+GAUSS( SEED >0. .SIN) 

T1-AHAX1<ANIN1<T1»1.  )  ,0.  > 

T2*AHAX1 (AM1N1<T2. 1, ) .0. > 

PI (HPN.N  > »T1 
P2<MPM,N)»T2 
219  CONTINUE 

220  IF(M.LE.KSl)  00  TC  230 
lASN-IASN-ILN(MA) 

IASD1-IASD1-ILD(NA,1)>  I ASD2* I ASD2-I LB ( MA , 2 ) 
IASDZ»IASDZ-ILDZ<MA> 

230  IS1~IS21*IS22«0 
IS2Z-0 

DO  249  N>l'JufN 


IS1*IS14IS 


1821*1821418 


IS22-IS224IS 

249  CONTINUE 

ILN<HA)-IS1 

ILD<HAt 1 >*IS21 I  ILD<NA.2)-IS22 
ILDZ(MA»IS2Z 

260  IF (M. LT.KS1 )  00  TO  499 
IF(H.OT.KSl)  GO  TO  280 
IASN«IASD1*IASD2«0 
IA3DZ-0 

DO  269  HAN-1. JNIX 
IASN-I ASN4ILN  <  NAN  > 

IASD1*IASD14ILD<MAH.1> I  IASD2*I ASD24ILD ( MAN. 2 ) 
IASDZ- I ASDZ4ILDZC  HAM ) 

269  CONTINUE 
00  TO  300 


280  IASN"XASN4ILN(HA) 

r5ISi*i^iS1  +  ILD<HA’1>l  IASD2-IASD2+ILD<MAf2> 
IASDZ-IASDZ+ILDZ<HA> 

290  IF(h.LE.KS2>  GO  TO  300 
HH-HW-KF 

IFIHH.LE.O)  HH-HH+NW 

FR»FR-WR<MM)*RF 

FRZ-FRZ-WRZ<MM)*RF 

300  T1»T4*IASN 

IF  <  T1 . GE . 0 . >  GO  TO  301 
T2*IASDli  GO  TO  302 

301  T2«IASD2 

302  73«ABS<T2> 

IF<T3.GT..3*ABS<T1>.AND,T3.GE.S.>  GO  TO  303 

GO  TO  304 

303  KE*KE+1 

304  WR(HW) 

T2=IASDZI  T3*ABS<T2> 

IF(T3.GT.ABS(T4) .AND.T3.GE.5. >  GO  TO  305 
.  GO  TO  306 

305  KEZ*KEZ41 

306  URZ  <  HU ) 

IF  <  N .LT  *  KS2 )  GO  TO  499 
HH*NWP 

IF (H.GT . KS2)  GO  TO  330 
FR=0. 

FRZ=0. 

BO  309  1=1, KF 
HM*HH-1  i  FR=FR+WR<HM> 

FRZ=FRZ+URZ(HM> 

309  CONTINUE 
FR=FR»..- 
FRZ=FRZtRF 
GO  TO  340 

330  FR=FR+WR<MW)*RF 

FRZ=FRZ+WRZ<MU)*RF  f., 

340  Sf*rR-Dl 
FA-EA+S1 
E2«E2*S1*S1 
S2=A8S<EMX) 

lF<ABS<sn.GT.S2>  EHX=S1 
S'ZtFRZ-DU  EAZ=EAZ+S1Z»  E2Z=E2Z+S1Z*S1Z 
499  CONTINUE**2”  IF  (  ABS<  S1Z  >  •  GT  •  S2Z  )  EMXZ=S1Z 
SP-SP/NP 

WRITE ( 6 f 6001 )  PEiDIfSPfSIPfNS 

6001  FORHAT*'  THE  SCENES  WERE  GENERATED  WITH  ' fF5 , 3f 

>  PROBABILITY  OF  EDGE  '/'  RELATIVE  SHIFT  IS  'fF5.3f'  PIXELS'/ 
X.Fsf*?  INTENSITY  IS  ',F5'3""  STD  DEV  IS  ' fF5.3f/'  1. /SNR  IS  ' 
RNE=1 ,/NE 

WRITE ( 6 . 6002  >  KDfHE.KF 

6002  FORHATI///'  PIPELINE  CORRELATION  RESULTS!  '/'  DIGITIZING  TO  ' 

®ITS  f  FOR.ISr'  LINES'//'  RUNNING  AVE  ', 

>'LENGTH  ',11/'  WINDOW') 

I1»JWIN-1 

I1»I1+1 

EA-EA4RNEI  E2-E24PNEI  SDF=S0RT<E2-EA*EA> (  E2=SQRT(E2>  ' 
55?*~*E*R*E*  E2Z=E2Z*RNEt  SDFZ=SQRT(E2Z-EAZ»EAZ> ,  E2Z»SQRT(£2Z) 

GO  TO  1 
9999  STOP 
END 


16.22.11  >los 

328.81  ARU'S,  .85  CONNECT  HRS 
LOGGED  OFF  AT  16.22.33  ON  07MAR82T( 


